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Abstract 

This  thesis  presents  a  method  for  building  static  Radar  Cross  Section  (RCS)  models 
of  aircraft  based  on  static  data  taken  from  limited  aspect  angle  windows.  These  models 
statistically  characterize  static  RCS.  This  is  done  to  show  that  a  limited  number  of  samples 
can  be  used  to  effectively  characterize  static  aircraft  RCS. 

The  aspect  angles  that  the  data  are  taken  from  are  nose-on  (0°±5°),  broadside 
(90°±5°),  and  tail-on  (180°±5°).  These  aspects  are  chosen  because  they  have  the  greatest 
probability  of  being  presented  to  the  point  defense  radar.  However,  any  aspect  angle  with 
any  window  size  can  be  modeled  with  the  method  presented  here. 

The  optimum  models  are  determined  by  performing  both  a  Kolmogorov  and  a  Chi- 
Square  goodness-of-fit  test  comparing  the  static  RCS  data  with  a  variety  of  probability 
density  functions  (pdf)  that  are  known  to  be  effective  at  approximating  the  static  RCS  of 
aircraft.  The  optimum  parameter  estimator  is  also  determined  by  the  goodness-of-fit  tests 
if  there  is  a  difference  in  pdf  parameters  obtained  by  the  Maximum  Likelihood  Estimator 
(MLE)  and  the  Method  of  Moments  (MoM)  estimators.  If  solving  for  the  MLE  results  in 
transcendental  equations,  numerical  methods  are  used  to  obtain  the  parameter  estimate. 


RADAR  CROSS  SECTION  MODELS  FOR 
LIMITED  ASPECT  ANGLE  WINDOWS 


/.  Introduction 


1.1  Overview 

In  the  problem  of  detection  of  aircraft  by  a  radar,  statistical  models  of  the  radar 
cross  section  (RCS)  are  necessary  (Skolnick,  1980:46).  The  models  must  take  into  account 
the  high  rate  of  RCS  fluctuation  and  the  large  magnitude  of  the  fluctuations  that  complex 
shapes,  such  as  aircraft,  exhibit  as  a  function  of  aspect  angle.  The  fluctuating  characteristic 
of  RCS  can  be  treated  as  a  random  variable  and  therefore,  probabilty  density  fuctions  (pdf) 
can  be  used  to  model  aircraft  RCS  (Dowdy,  1991:164).  The  parameters  of  the  pdfs  can  be 
estimated  using  either  measured  RCS  data  from  the  aircraft  itself  or  theoretical  RCS  data 
generated  from  the  electromagnetic  reflections  of  combinations  of  simple  point  scatterers 
(Maffett,  1989:239). 

A  high  rate  of  amplitude  fluctuations  requires  a  correspondingly  high  sampling  rate 
to  model  lobe  pattern.  The  Radar  Target  Scattering  range  (RATSCAT,  Holloman  AFB, 
NM),  an  Air  Force  organization  responsible  for  measuring  the  RCS  of  aircraft,  and  the 
sponsor  of  this  research,  is  currently  required  to  sample  the  RCS  of  static  aircraft  every 
0.1°  to  meet  characterization  requirements.  RATSCAT  would  like  to  determine  if  the  total 
number  of  samples  can  be  reduced  without  affecting  the  characterization.  This  thesis  will 
show  that  good  RCS  models  can  be  built  from  limited  aspect  windows,  thus  drastically 
reducing  the  total  number  of  data  samples  required. 

1.2  Problem  Statement 

The  goals  of  this  thesis  are  to  derive  a  method  for  building  RCS  models  using  a 
minimum  number  of  data  samples  and  to  determine  if  there  is  a  consistently  superior 
method  of  parameter  estimation  for  the  models.  Minimizing  the  value  of  the  test  statistic 
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in  goodness-of-fit  tests,  for  both  parameter  estimators  of  pdfs  and  the  pdfs  themselves, 
should  result  in  acceptable  RCS  models. 

1.3  Assumptions 

1.  The  aircraft  RCS  data  used  is  from  a  scale  model  of  a  jet  fighter,  coated  with 
metallic  paint.  It  is  assumed  that  the  aircraft  is  completely  symmetrical,  so  only  180°  (one 
side)  of  the  aircraft  is  sampled.  The  missile  RCS  data  is  obtained  from  an  actual  AIM-9 
missile.  However,  measurement  errors  distorted  the  amplitude  of  one  side  of  the  missile 
relative  to  the  other  side. 

2.  The  sampling  intervals  are  assumed  to  be  adequate  for  RATSCAT  standards.  The 
aircraft  data  is  sampled  at  0.5°  intervals,  the  missile  is  sampled  at  0.25°  intervals. 

3.  The  Maximum  Likelihood  Estimator  (MLE)  and  the  Method  of  Moments  (MoM) 
are  the  two  most  effective  parameter  estimators. 

4.  The  aircraft  and  missile  data  are  representative  of  their  class  of  object.  Data  from 
more  than  one  type  of  jet  fighter  or  missile  would  be  helpful  in  validating  the  effectiveness 
of  the  results,  but  were  unavailable. 

1.4  Scope 

This  study  will  investigate  the  sampling  requirements  for  building  a  static  RCS 
model.  No  account  of  dynamic  properties  effecting  RCS  such  as  wing  flutter,  pitch  and 
roll,  or  atmospheric  conditions  will  be  taken  into  account. 

The  parameters  will  be  estimated  using  the  Maximum  Likelihood  Estimator  and  the 
Method  of  Moments. 

The  study  will  consider  two  types  of  targets:  a  jet  fighter  aircraft  and  an  air-to-air 
missile.  This  will  be  done  to  compare  the  statistical  characteristics  of  the  two  types  of 
targets. 
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1.5  Approach 

A  computer  program  will  form  the  pdfs  using  RCS  data.  The  program  will  compute 
the  MoM  and  MLE  parameters.  Theoretical  values  for  MoM  and  MLE  parameters  will  be 
determined  in  Chapter  III. 

The  program  will  also  determine  the  test  statistics  for  the  Kolmogorov  and  Chi- 
Square  goodness-of-fit  tests.  Tables  are  then  used  to  determine  the  Level  of  Significance 
for  the  parameter  estimates,  pdfs,  and  data  sets. 

Trends  in  success  of  types  of  pdfs  and  parameter  estimators  in  relation  to  sample 
statistics  will  be  analyzed  in  Chapters  IV  and  V. 

1.6  Sequence  of  Presentation 

Chapter  II  presents  background  material  related  to  the  thesis  topic.  The  chapter 
also  provides  support  for  some  of  the  limitations  stated  here. 

Chapter  III  outlines  the  methodology  used  and  specific  assumptions  implicit  in  com¬ 
putations  required  to  determine  best  fit. 

Chapter  IV  provides  analysis  of  RCS  models  built  from  limited  aspect  window  versus 
full-range  sample  sets. 

Chapter  V  provides  analysis  of  RCS  models  built  from  limited  aspect  window  versus 
limited  aspect  window  sample  sets. 

Chapter  VI  presents  conclusions  drawn  from  the  study  and  recommends  topics  for 
further  study. 
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II.  Background 


2. 1  Overview 

This  chapter  provides  a  background  of  radar  and  statistical  concepts  used  in  the 
topic  development  in  the  following  chapters. 

The  radar  range  equation  is  the  primary  means  of  evaluating  radar  performance 
based  on  the  power  of  the  transmitted  signal,  the  gain  of  the  antenna,  the  wavelength  of 
the  transmitted  signal,  the  power  of  the  received  signal,  and  the  radar  cross  section  (RCS) 
of  the  target.  Of  all  the  parameters  of  the  radar  range  equation,  the  RCS  is  the  most 
difficult  to  specify  because  it  is  constantly  fluctuating  in  a  random  manner.  RCS  depends 
upon  the  size  and  surface  of  the  target,  aspect  angle  of  the  target  to  the  radar,  motion  of 
the  target  relative  to  the  radar,  frequency  of  the  transmitted  signal,  and  so  on  (Skolnik, 
1980:41). 

A  radar  system  can  detect  and  classify  an  unknown  target  by  comparing  character¬ 
istics  of  the  received  signal  with  a  dictionary  of  optimized  statistical  models  of  various 
types  of  targets.  The  statistical  models  are  developed  from  RCS  data  of  static  models  of 
the  aircraft.  Statistical  characteristics  of  the  data  are  used  to  estimate  the  parameters  of 
probability  density  functions  (pdf). 

This  chapter  will  review  RCS  characteristics,  types  of  probability  density  functions 
used  to  make  RCS  models,  methods  of  parameter  estimation,  goodness-of-fit  tests,  and 
radar  detection  techniques. 

2.2  RCS  Characteristics 

When  an  electromagnetic  wave  travelling  through  one  medium  strikes  the  surface  of 
a  second  medium,  a  portion  of  the  energy  of  the  wave  is  refracted  into  the  second  medium, 
and  the  rest  of  the  energy  is  reflected  back  into  the  first  medium.  The  amount  of  energy 
that  is  enters  each  medium  depends  on  two  variables:  the  angle  of  incidence  the  wave  has 
with  the  tangent  line  of  the  point  of  incidence  on  the  second  medium  and  the  impedance  of 
the  two  mediums.  Snell’s  law  of  reflection  states  that  the  angle  of  reflection  from  the  line 
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perpendicular  to  the  tangent  of  the  point  of  incidence  equals  the  angle  of  incidence  to  the 
tangent.  However,  the  angle  of  reflection  refers  to  the  opposite  side  of  the  perpendicular  line 
from  the  angle  of  incidence.  Therefore,  the  only  time  energy  will  be  reflected  back  toward 
the  origin  of  the  wave  (the  radar)  will  be  if  it  strikes  the  second  medium  perpendicularly. 
This  reflected  energy  would  be  the  signal  a  monostatic  (one  antenna)  radar  would  receive. 

For  a  perpendicular  wave,  the  magnitude  of  the  reflected  electrical  field  is 


Er 


(Z?  —  Zi)Ej 
Z2  T  Z\  ’ 


(2.1) 


where  Er  is  the  reflected  electromagnetic  field,  £,  is  the  incident  electromagnetic  field,  Z2 
is  the  impedance  of  the  second  medium,  and  Z\  is  the  impedance  of  the  first  medium.  A 
negative  value  indicates  a  wave  travelling  toward  the  origin  of  the  incident  wave  (Kraus: 
1984,  Ch.  10  and  12).  For  a  radar  signal,  the  first  medium  is  air  and  the  second  medium 
is  the  target  (aircraft).  Most  aircraft  have  a  metallic  (low  impedance)  skin.  Air  has  a  high 
impedance.  Therefore,  most  of  the  electromagnetic  wave  is  reflected  away  from  the  aircraft. 
Link  states  that  ’’the  Radar  Cross  Section  (RCS)  of  a  target,  simply  put,  is  a  measure  of 
the  electromagnetic  energy  reflected  from  the  target  divided  by  the  electromagnetic  energy 
incident  upon  the  target”  (Link,  1983:1-1).  A  complex  target,  such  as  an  aircraft,  has  a 
variety  of  surface  types.  An  electromagnetic  wave  incident  upon  it  will  be  scattered  in 
many  directions.  Some  portion  of  the  scattering  surfaces  will  reflect  the  wave  back  toward 
the  radar.  The  variety  of  types  of  surfaces  and  distances  from  the  radar  of  the  individual 
scattering  surfaces  will  cause  the  reflected  waves  to  possess  different  phases  and  amplitudes. 
Skolnik  states 


The  relative  phases  and  amplitudes  of  the  echo  signals  from  the  individual 
scattering  objects  as  measured  at  the  radar  receiver  determine  the  total  cross 
section.  The  phases  and  amplitudes  of  the  individual  signals  might  add  to  give 
a  large  total  cross  section,  or  the  relationships  with  one  another  might  result 
in  total  cancellation  (Skolnik,  1980:38). 

The  interaction  of  the  multiple  scattering  surfaces  causes  the  phase  of  the  received 
signal  to  fluctuate  much  more  rapidly  than  the  magnitude,  as  aspect  angle  is  changed 
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(Knott,  et  al.,  1985:180).  For  this  reason,  the  phase  of  the  received  signal  is  usually 
modeled  as  a  random  variable  having  a  uniform  pdf  over  the  interval  [0,  27r]  (Crispin  and 
Seigel,  1968:395),  (Maffett,  1989:218),  though  this  is  not  always  the  case  (Sandhu  and 
Saylor,  1984:490-507).  Assuming  a  uniform  phase  distribution  and  independent  scattering 
surfaces,  each  surface’s  phase  may  be  averaged  over  individually.  RCS  fluctuates  at  a 
slower  rate,  but  it  is  still  practical  to  treat  it  as  a  random  variable  (Skolnik,  1980:46). 


2.3  Modeling  RCS  with  Probability  Density  Functions 

Certain  types  of  pdfs  have  been  found  that  do  a  good  job  of  modeling  the  RCS. 
The  most  popular  pdfs  are  Rayleigh  functions  Swerling  used  in  his  four  fluctuation  models 
(Swerling,  1960:268-308).  One  pdf  modeled  a  complex  target  consisting  of  many  indepen¬ 
dent  scatterers  of  approximately  equal  echoing  areas: 


p(*|a)  =  ^exp(^), 

x  >  0, 


(2.2) 

(2.3) 


where  x  is  the  RCS  and  a  is  a  parameter  that  will  be  estimated  from  the  RCS  statistics. 
The  other  pdf  modeled  one  large  reflector  together  with  other  small  reflectors  and  is  named 
the  One- Dominant-Plus- Rayleigh: 


9a:3  .  — 3x2 

Hx|a)  =  _eXp(— ), 
x  >  0. 


(2.4) 

(2.5) 


(DiFranco  and  Rubin,  1980:Ch.  11).  These  pdfs  are  modifications  of  the  chi-square  density 
function  of  degree  of  freedom  2  and  4  respectively. 

Other  pdfs  that  have  been  found  to  exhibit  good  fit  to  experimental  data  are  the 
Beta,  Lognormal,  Normal,  Rayleigh  Squared,  and  Weibull  (Maffet,  1989:255),  (Dowdy, 
1991:165)  and  are  listed  in  Table  2.1.  In  this  table,  x  is  RCS,  y  is  RCS  squared,  and  a,  /?, 
p,  and  o  are  parameters  which  are  estimated  from  the  RCS  statistics. 
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Name 

PDF 

Range 

Beta 

**!<>■/»)= 

0<x<l 

a,P  >  0 

Lognormal 

0  >  0 

Normal 

WiUUIItMiTBBMBfl 

-oo<x<oo 

Rayleigh  Squared 

0 

Al 

a 

Wei  bull 

x  >  0 
a,/?  >  0 

Table  2.1.  Probability  Density  Functions 


2.4  Parameter  Estimation 

To  optimize  the  fit  of  pdfs  to  the  data,  the  parameters  of  the  pdfs  must  be  optimized. 
The  job  of  a  parameter  estimator  is  to  estimate  0  (the  parameter  of  the  pdf)  by  determining 
0  (the  parameter  of  the  observed  data).  Two  of  the  most  common  methods  used  to  estimate 
parameters  are  the  Method  of  Moments  (MoM)  and  the  Maximum  Likelihood  Estimator 
(MLE).  The  MLE  method  declares  0  to  be  the  values  of  0  that  gives  the  likelihood  function 
p(x\0)  its  maximum  value.  This  is  done  by  taking  the  first  derivative  of  the  likelihood 
function  and  setting  it  equal  to  zero: 

dp{x\e ) 

de 

and  solving  for  0  (Cramer,  1946:499).  The  second  derivative  should  be  examined  to  deter¬ 
mine  whether  6  is  actually  a  maximum  or  a  minimum.  Cramer  states 


The  importance  of  the  method  is  clearly  shown  by  the  two  following  proposi¬ 
tions: 

If  an  efficient  estimate  0  of  6  exists,  the  likelihood  equation  will  have  a  unique 
solution  equal  to  0. 

If  a  sufficient  estimate  0  of  0  exists,  any  solution  of  the  likelihood  equation  will 
be  a  function  of  ^(Cramer,  1946:499). 
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An  efficient  estimate  is  an  unbiased  estimate  that  satisfies  the  following  equation: 


alnP(_l|l>  =  kis-i) 
de  y  ’ 


where  A:  is  a  non-zero  function  of  6.  A  sufficient  estimate  is  an  unbiased  estimate  that 


satisfies  the  following  equation: 


d\np(x\0)  _  Q 


Since  the  maximum  of  \np(z\6)  occurs  for  the  same  value  of  9  as  the  maximum  of  p(x\0), 
the  Omle  will  be  an  efficient  and  sufficient  estimator  of  6  if  9Mle  is  unbiased.  An  estimator 
is  unbiased  if  the  expected  value  of  6  is  0: 

E(S)  =  9.  (2.9) 

Law  and  Kelton  list  these  desirable  statistical  properties  of  MLEs: 


1.  For  most  of  the  common  distributions,  the  MLE  is  unique;  that  is,  L(/3) 
is  strictly  greater  than  L(/3 )  for  any  other  value  of  (3  ( L(/3 )  is  the  likelihood 
function  of  /3). 

2.  Although  MLEs  need  not  be  unbiased,  in  general,  the  asymptotic  distribu¬ 
tion  (asn-»  oo)  of  9  has  mean  equal  to  9  (see  property  4  below). 

3.  MLEs  are  invariant;  i.e.,  if  <f>  =  h{9)  for  some  function  h ,  then  the  MLE 
of  (f>  is  h{9).  (Unbiasedness  is  not  invariant).  For  example,  the  variance  of  an 
exp (/3)  random  variable  is  /32,  so  that  the  MLE  of  this  variance  is  [x(n)]2. 

4.  MLEs  are  asymptotically  normally  distributed;  i.e.,  y/n(9  —  9)  N( 0,  u(0)), 

where  v(0)  =  (the  expectation  is  with  respect  to  x,-,  assuming  that  x, 

has  the  hypothesized  distribution)  and  denotes  convergence  in  distribution. 
Furthermore  if  9  is  any  other  estimator  such  that  y/n(9  —  9)  N( 0,cr2),  then 
v(9 )  <  a2.  (Thus,  MLEs  are  called  best  asymptotically  normal.) 

5.  MLEs  are  strongly  consistent;  that  is,  limn_inf  6  —  9  (with  probability  of 
l)(Law  and  Kelton,  1982:191). 


The  Method  of  Moments  determines  9,  by  determining  the  first  i  moments  of  the 
density  function  (p<)  and  equating  each  of  those  moments  to  the  corresponding  sample 
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moments  ( m* )  of  the  data,  n  multiple  parameters  can  be  determined  by  solving  n  equations 
simultaneously.  The  MoM  is  widely  used  mainly  because  it  is  easy  to  impliment.  However, 
it  was  shown  by  Fisher  that  the  asymptotic  efficiency  of  the  MoM  is  not  very  good  (Fisher, 
1921:309).  This  would  indicate  that  the  MoM  does  not  generally  give  parameters  with 
the  smallest  possible  variance  from  the  true  parameters  (Cramer,  1946:490).  The  superior 
efficiency  of  the  MLE  would  tend  to  indicate  that  the  MLE  should  be  chosen  over  the 
MoM  as  an  estimator  if  possible.  However,  for  pdfs  with  multiple  parameters,  the  MLE 
often  results  in  transcendental  equations,  in  which  case  numerical  methods  must  be  used  to 
determine  the  parameters.  The  MoM  can  then  be  used  as  the  initial  value  of  an  iteration 
process  (Shenton  and  Bowman,  1977:162).  Also,  if  the  MLE  parameters  are  arrived  at 
numerically,  depending  on  the  procedure,  the  arrived  at  parameters  may  not  have  an 
efficiency  as  high  as  the  actual  MLE  estimate.  The  MoM  by  itself  may  then  give  a  better 
estimate.  Finally,  there  are  many  pdfs  for  which  the  MoM  and  the  MLE  give  an  identical 
value  of  d. 

2.5  Goodness- Of- Fit 

In  the  last  section,  parameter  variance  was  mentioned  as  a  tool  for  determining,  a 
priori,  what  method  of  parameter  estimation  to  use  to  fit  a  pdf  to  the  observed  data.  To 
determine  if  the  observed  data  really  could  have  been  obtained  by  sampling  from  the  fitted 
pdf,  a  goodness-of-fit  test  is  used.  This  test  shows  how  well  a  fitted  pdf  fits  the  observed 
data  (Lewis  and  Orav,  1989:140).  A  goodness-of-fit  test  works  by  performing  a  hypothesis 
test  on  the  pdf  with  a  null  hypothesis  of: 

H0:  The  x,’s  are  independent  identically  distributed  (IID)  random  variables 

with  distribution  function  F. 

The  test  then  finds  a  test  statistic  by  making  some  type  of  comparison  between  the  fitted 
distribution  and  the  observed  data.  The  test  statistic  is  compared  with  a  range  of  possible 
values  of  the  test  statistic  for  a  given  number  of  samples.  Where  the  actual  test  statistic 
falls  in  the  range  determines  the  level  of  significance  a,  which  is  the  maximum  probability 
of  rejecting  the  null  hypothesis.  Therefore,  a  higher  value  of  a  corresponds  to  a  better  cor- 
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relation  between  the  data  and  the  distribution.  There  are  no  set  rules  for  determining  what 
an  appropriate  Level  of  Significance  is,  but  a  commonly  accepted  figure  is  5%  (Otnes  and 
Enochson,  1978:61),  (Conover,  1980:Chapters  4  and  5),  (Law  and  Kelton,  1982:  Chapter 
5). 

The  two  most  commonly  used  goodness-of-fit  tests  are  the  Chi-Square  and  the  Kol¬ 
mogorov  tests.  The  Chi-Square  test  groups  the  ordered  data  into  frequency  histograms. 
The  test  statistic  is  given  by 

T  =  (2.10) 

where  O,  is  the  actual  number  of  observations  in  bin  i  and  Ei  is  the  expected  number  of 
observations.  Although  the  optimum  number  of  histogram  intervals  has  not  been  deter¬ 
mined,  Law  and  Kelton  recommend  that  the  probability  of  an  observation  occurring  in  an 
interval  should  be  as  close  to  identical  as  possible  for  all  intervals.  This  is  done  to  make 
the  test  unbiased.  They  also  recommend  that  the  expected  number  of  observations  Ei  be 
greater  than  or  equal  to  5,  so  that  the  agreement  between  the  true  distribution  of  T  and 
its  asymptotic  Chi-Square  distribution  will  not  be  too  small.  This  is  desirable  because,  if 
the  null  hypothesis  is  true,  T  converges  with  the  Chi-Square  distribution  as  n  -*  oo.  The 
Chi-Square  test’s  chief  asset  is  its  applicability  to  any  type  of  distribution  function  (Law 
and  Kelton,  1982:197). 

The  Kolmogorov  goodness-of-fit  test  is  the  appropriate  name  of  the  Kolmogorov- 
Smirnov  goodness-of-fit  test  when  it  is  applied  to  comparisons  between  distributions  and 
observed  data.  The  test  statistic  is  the  maximum  vertical  distance  between  the  observed 
data  and  the  hypothesized  distribution: 

T  =  supx|F(x)  -  S(x)|,  (2.11) 

where  F(x)  is  the  hypothesized  distribution,  S(x)  is  the  empirical  distribution  of  the 
observed  data,  and  supr  indicates  the  supremum  over  all  x  (Conover,  1980:296).  The  chief 
advantage  of  the  Kolmogorov  test  over  the  Chi-Square  test  is  that  it  gives  exact  answers 
regardless  of  the  sample  size,  while  the  Chi-Square  test  only  asymptotically  approaches 
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exactness  as  the  number  of  samples  approaches  infinity.  Therefore,  the  Kolmogorov  test 
is  particularly  useful  for  small  sample  sets  (Conover,  1980:296). 

Unfortunately,  Law  and  Kelton  state  that  the  Kolmogorov  goodness-of-fit  test  is  , 
in  general,  only  applicable  to  hypothesized  distributions  whose  parameters  are  determined 
independently  of  the  observed  data.  One  notable  exception  is  the  Weibull  distribution 
function  (Law  and  Kelton,  1982:202).  However,  Law  and  Kelton  further  state,  ”the  ef¬ 
fect  of  this  misapplication  (using  the  Kolmogorov  test  for  a  distribution  with  parameters 
estimated  from  data)  is  not  well  understood”  (Law  and  Kelton,  1982:199).  In  addition, 
Breiman  states 


The  effect  that  this  has  on  the  level  of  the  test  is  not  well  known.  The  evidence 
we  have  is  that  the  effect  is  not  very  important.  For  moderate  to  large  sample 
size,  it  is  probably  safe  to  ingore  the  fact  that  6  was  estimated  (Breiman, 
1973:213). 

There  are  also  many  examples  of  the  Kolmogorov  test  being  applied  to  distributions  with 
estimated  parameters  (Dowdy,  1991:164),  (Maffett,  1989:Ch.l2).  Because  the  test  is  valid 
for  the  Weibull  distribution  and  evidently  does  not  cause  big  problems  for  other  distribu¬ 
tions,  the  results  of  the  Kolmogorov  test  will  be  determined  along  with  the  Chi-Square 
test  in  this  study. 

2. 6  Radar  Detection 

This  section  demonstrates  how  RCS  models  using  pdfs  are  used  in  radar  detection 
determination.  Radar  detection  is  an  application  of  statistical  decision  theory.  DiFranco 
and  Rubin  state 


The  basic  elements  of  a  statistical  decision  problem  are:  (1)  a  set  of  hypotheses 
that  characterize  the  possible  true  states  of  nature;  (2)  a  test  in  which  data 
are  obtained  from  which  we  wish  to  infer  the  truth;  (3)  a  decision  rule  that 
operates  on  the  data  to  decide  in  an  optimal  fashion  which  hypothesis  in  fact 
best  describes  the  true  state  of  nature;  (4)  a  criterion  of  optimality  that  reflects 
the  cost  of  correct  and  incorrect  decisions. 
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p(y) 


0  yO 

y 

Figure  2.1.  Neyman-Pearson  Decision  Test 

(DiFranco  and  Rubin,  1980:253).  For  a  radar  detection  problem,  the  hypotheses  are  usually 

H0  :  Noise  only  present 

Hi  :  Signal  plus  noise  is  present 

There  could  be  further  hypotheses  if  a  distinction  is  to  be  made  over  the  type  of  signal 
present. 

The  test  used  to  determine  which  hypothesis  is  true  is  called  the  Likelihood  Ratio 
Test.  In  this  test,  a  likelihood  ratio  is  formed  by  the  pdfs  of  the  two  hypotheses  and  is 
compared  with  a  threshold.  If  the  ratio  is  larger  than  the  threshold,  the  hypothesis  whose 
pdf  formed  the  numerator  of  the  ratio  is  declared  to  be  true.  If  the  ratio  is  less  than  the 
threshold,  the  other  hypothesis  is  declared  to  be  true  (Melsa,  1978:25). 

There  are  a  variety  of  methods  to  determine  what  the  threshold  should  be.  For 
radar  applications,  the  Neyman-Pearson  Test  is  usually  used.  This  test  chooses  a  value 
that  maximizes  the  probability  of  detection  Pd,  for  a  given  probability  of  false  alarm  PJa 
(Maffett,  1989:275).  This  method  can  be  understood  by  examining  Figure  1.  Assume  pi(y) 


2-9 


is  the  pdf  of  a  signal  present  plus  noise.  Assume  po{y)  is  the  pdf  of  just  noise.  If  Pja  is 
given  a  value  of  a0,  then 

Po(y)dy  (2.12) 

where  y  is  the  magnitude  of  the  received  waveform.  a0  and  p0(y)  are  known,  so  y0  can  be 
solved  for.  The  value  of  the  threshold  A  is  then 


Pi(Po) 

Po(lfo)’ 


and  the  complete  Ney man- Pearson  test  is 


(2.13) 


A(y) 


Pi(y) 

Po(y) 


> 

H0  A, 


(2.14) 


where  A (y)  represents  the  likelihood  ratio  test  (LRT).  This  test  reads:  choose  Hx  if  the 
LRT  is  greater  than  A,  choose  H0  if  the  LRT  is  less  than  A.  Pd  can  be  found  by  integrating 
Pi(y)  over  the  region  of  P]a : 


Pi  =  /  Pi(y)<*y  (2.15) 

Jyo 

(Maffett,  1989:275-277). 

If  the  a  priori  probabilities  of  the  two  hypotheses  (P[/f0]  and  P[/fi])  are  known,  a 
method  called  the  Bayes  Risk  test  can  be  used  to  determine  the  threshold.  Besides  P[/f0] 
and  P[#i],  the  costs  of  each  of  the  four  possible  choices  from  Figure  1  must  be  known. 
The  variables  for  the  possible  costs  are  as  follows: 


Coo  :  Cost  of  correct  no  target 

C0i  :  Cost  of  miss 

Cio  :  Cost  of  false  alarm 

Cn  :  Cost  of  correct  target  present 
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Usually,  C0o  and  Cn  have  a  cost  of  0.  If  Coo  =  Cn  =  0  and  C01  =  C10,  the  test  is  called 
the  ideal  observer  test.  The  objective  of  the  Bayes  Risk  criterion  is  to  minimize  Bayes 
Cost  ( B ): 


B  =  Cw?[d0tHo\  +  Colnd0,Hl]  +  CloP[duHo]  +  CllP[dl,Hl}  (2.16) 


Assuming  Cn  =  C0o  =  0,  Melsa  and  Cohn  show  the  LRT  for  the  Bayes  Risk  is 


A  (y) 


p(y\Hi)  |  c10P[^0] 
p(y\H0)  °  CoxPfflx] 


(2.17) 


(Melsa  and  Cohn, 1978:45).  The  idea  behind  the  test  is  that  even  though  the  received 
signal’s  amplitude  would  indicate  one  hypothesis,  it  is  much  more  important  to  correctly 
identify  the  other  hypothesis,  so  the  threshold  is  skewed  to  compensate.  As  an  example, 
assume  the  job  of  the  LRT  is  to  perform  target  identification  and  the  two  hypotheses  are 


#0  :  Target  is  a  missile 
Hx  :  Target  is  a  fighter 

The  cost  variable  would  then  be 


C00  :  Cost  of  correctly  deciding  missile 

C0i  :  Cost  of  deciding  missile  present  when  fighter  present 
C\0  :  Cost  of  deciding  fighter  present  when  missile  present 
Cn  :  Cost  of  correctly  deciding  fighter 


Assume  the  cost  of  incorrectly  identifying  a  missile  (Cjo)  is  five  times  higher  than  the  cost 
of  incorrectly  identifying  a  fighter  (C0i),  and  that  Coo  =  Cn  =  0.  Also,  assume  the  a  priori 
probabilities  of  Hx  and  H0  are  equal.  Then  the  LRT  would  be 


A  (y) 


p(y\Hi) 

p(y\n0) 


H| 

> 

< 


H 05. 


(2.18) 
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Most  radar  receivers  are  designed  to  perform  target  detection  based  on  a  number  of 
observations  of  a  randomly  fluctuating  waveform.  The  receivers  are  usually  designed  to 
optimize  detection  of  waveforms  that  behave  according  to  one  of  the  four  Swerling  cases 
for  randomly  fluctuating  waveforms  (Skolnik,  1980:46): 

Swerling  I:  Amplitude  of  an  entire  pulse  train  is  constant  but  behaves  as  a  single 
random  variable  with  a  Rayleigh  pdf,  from  scan-to-scan.  The  phase  of  each  pulse  is  assumed 
to  be  a  statistically  independent  random  variable  with  a  uniform  pdf. 

Swerling  II:  Similar  to  Swerling  I  except  the  amplitude  of  each  pulse  in  the  pulse 
train  is  a  statistically  independent  random  variable  with  a  Rayleigh  pdf. 

Swerling  III:  Similar  to  Swerling  I  except  the  amplitude  of  the  pulse  train  is  a  random 
variable  with  a  one-dominant-plus- Rayleigh  pdf. 

Swerling  IV:  Similar  to  Swerling  III  except  the  amplitude  of  each  pulse  in  the  pulse 
train  is  a  statistically  independent  random  variable  with  a  one-dominant-plus-Rayleigh 
pdf. 

DiFranco  and  Rubin  show  that  the  optimum  detection  receiver  for  these  waveforms  is 
a  single  pulse  matched  filter  followed  by  an  envelope  detector,  a  sampler,  an  integrator,  and 
finally,  a  threshold  comparator  (DiFranco  and  Rubin,  1980:Ch.l0  and  11).  Kulp,  Maffett, 
and  Dowdy  each  show  that  the  Rayleigh  and  the  one-dominant-plus-Rayleigh  pdfs  are  not 
always  the  optimum  pdfs  for  the  Swerling  cases  (Kulp,  1984:165-183),  (Maffett,  1989:280- 
285),  (Dowdy,  1991:167). 

2. 7  Summary 

This  chapter  presented  a  description  of  the  nature  of  RCS  and  how  it  is  used  in  radar 
det  ction  systems.  In  doing  so,  the  importance  of  the  pdf  was  demonstrated,  along  with 
methods  of  optimizing  the  pdf  in  the  Likelihood  Ratio  Test. 
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III.  Methodology 


3.  1  Introduction 

This  chapter  introduces  a  method  to  determine  the  optimum  pdf  and  parameter 
estimator  for  any  range  of  aspect  angle  needed.  Assumptions  made  will  be  explained  as 
they  arise.  The  MoM  and  MLE  estimators  will  be  determined  for  all  pdfs  considered  in 
this  thesis.  The  RCS  model  generating  program  ’’ASPECT”  is  presented  along  with  the 
required  inputs. 

3.2  Parameter  Determination 

As  explained  in  Chapter  II,  parameters  will  be  estimated  for  the  Rayleigh,  One- 
Dominant- Plus- Rayleigh,  Rayleigh  Squared  (negative  exponential),  Lognormal,  Normal, 
Beta,  and  Weibull  pdfs. 

3.2.1  Rayleigh  The  Rayleigh  pdf  for  n  independent  identically  distributed  (IID) 
data  samples  is  defined  as  (DiFranco  and  Rubin,  1968:376) 

"  x  —x? 

p(x!,x2,...,xn|a)  =  I]^exP(^)  (3-1) 

for  x  >  0.  To  determine  the  MLE  estimate  of  a,  a  sufficient  statistic  is  found  by  taking 
the  natural  logarithm  of  both  sides  of  the  equation,  giving: 

n  1  n 

lnp(xi,x2,...,x»  =  ln(J^x,)-2nlna-  —  ^x<2.  (3.2) 

«=1  i=l 

Taking  the  first  derivative  with  respect  to  a  and  setting  it  equal  to  zero  gives: 

51np(x1,x2,...,x„|a)  -2n  1  ^  2 

- az -  =  — +  ^L*<  <3-3) 

•=i 

=  0. 
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Solving  for  a  gives: 


&mle  = 


The  MoM  estimate  of  a  is  found  by  first  determining  the  second  moment  of  the  single 
sample  pdf: 


f°°  x3  — x 2 


The  solution  to  this  definite  integral  is  found  using  the  CRC  integral  #662: 


/°°  n  ,  ,v.  r (k) 

I  x  expi—  ax  )dx  =  — — , 

Jo  pak 


where  n  >  — 1  ,p  >  0, a  >  0,  and  k  =  Thus, 


u  =  /  1  x  P(2) 

=  2a2. 


Setting  p2  equal  to  the  sampled  second  moment  and  solving  for  a  gives: 


ocmoM  =  (^£>.-2)*» 


which  is  identical  to  &mle- 


3.2.2  Rayleigh  Squared  This  pdf  is  found  by  transforming  the  random  variable  x 
to  the  random  variable  y  =  x2.  This  is  done  by  substituting  y  into  the  Rayleigh  pdf  for  x 
and  multiplying  the  result  by  the  Jacobian  of  the  transformation, 


=  nSexP(^§-)llrl 


A*a2  rv2a2/,0x' 
»=1 

=  ndrxp(i5) 

1=1 


(3.10) 


(3.11) 
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The  MLE  estimate  of  a  is  found  by  taking  the  natural  logarithm  of  both  sides,  which  is 


1  " 

lnp(ii,i2,-,*n|a)  =  -nln2  -  2nlna  -  —  £  &2. 


(3.12) 


1=1 


Taking  the  first  derivative  with  respect  to  a  and  setting  the  result  equal  to  zero  gives: 


glnp(j1,x2,...,xn|a) 

da 


(3.13) 


=  0. 


Solving  for  a  gives: 


(3.14) 

(3.15) 


which  is  identical  to  the  MLE  for  the  Rayleigh  pdf. 

The  procedure  for  finding  the  MoM  estimate  of  a  is  identical  to  the  Rayleigh  MoM 
procedure,  however  the  variables  in  Equation  24  are  n  =  5,  a  =  p  =  2,  and  k  =  3. 
Thus, 


A*2 


l«z;2(^)3 

8a4. 


(3.16) 

(3.17) 


Setting  p2  equal  to  the  sampled  second  moment  and  solving  for  a  gives: 

Omom  =  ■  (3.18) 

8nr=t 

Unlike  the  Rayleigh  pdf,  the  Rayleigh  Squared  pdf  does  not  have  identical  MLE  and  MoM 
estimates. 
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3.2.3  One- Dominant- Plus- Rayleigh  For  n  IID  data  samples,  this  pdf  is  defined  as 

(DiFranco  and  Rubin,  1968:407) 

"  gx3  -3x-2 

p(i1,x2,...,i»  =  n2a4exp(  2a2 

»=i 

(3.19) 

n  >  0. 

Taking  the  natural  logarithm  of  both  sides  of  the  equation  gives: 

n  g  n 

lnp(xi,x2,  ...,xn|o)  =  nln9  -  nln2  -  4nlno  4-3^1nx<  -  — 

i=i  *a  )=1 

(3.20) 

Taking  the  derivatives  of  both  sides  with  respect  to  a  gives: 

dlnp(x!,x2,...,xn|a)  _  -An  3  , 

da  '  a  a3  ‘  ' 

(3.21) 

Setting  the  result  equal  to  zero  and  solving  for  a  gives: 

&mle  = 

(3.22) 

As  with  the  Rayleigh  pdf,  aMoM  =  aMLE. 

3.2.4  Lognormal  For  n  IID  data  samples,  this  pdf  is  defined  as  (Law  and  Kelton, 

1982:164) 

,  1  ,  A  1  ,-(lnx,-  -  /i)2^ 

(3.23) 

H 

V 

O 
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which  is  the  transformation  y  =  lnx  applied  to  the  Normal  pdf.  Taking  the  natural 
logarithm  of  both  sides  gives: 


lnp(xi,x2,„.,xn|<7,p)  =  -nlna  -  ^ln2  -  ^ln*  -  ^lnx, 


t=i 


i—  1 


Taking  the  derivative  with  respect  to  \i  gives: 


dlnp(x1,x2,...,xn|cr,/i)  -1^„  , 

- ^ 


1  =  1 


Setting  the  derivative  equal  to  zero  and  solving  for  p  gives: 

1  n 

P^mle  —  —  y^lnx,. 

nU. 

Taking  the  derivative  with  respect  to  a  gives: 

dlnp(x1,x2,...,x„|(r,/i)  -n  1  _  .  ,2 

- Yo -  =  T  +  • 


*  *  r=  i 


Setting  the  derivative  equal  to  zero  and  solving  for  a  gives: 


(3.24) 

(3.25) 


(3.26) 


(3.27) 


(3.28) 


vmle  =  (“HOn*.  “A)2)*-  (3-29) 

n7^i 

The  MoM  estimates  are  found  by  transforming  the  lognormal  pdf  to  the  normal  pdf.  If 


Vi i  V2-)  •••?  Vn  —  lnii,lnx2,  ...,lnx„ 


(3.30) 


and  x1,x2,...,x„  are  lognormally  distributed  random  variables,  then  3/1, y2, ...,  yn  have  a 
normal  distribution.  By  definition, 


A MoM 


& 

1  n 

~X>*. 


(3.31) 

(3.32) 
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Pmle 

A2  m.M 

1  n 

=?>* 

»=1 

i]T(lnx,)2 
n  »= 1 

(3.33) 

(3.34) 

&MoM  =  (A2m.m  “  PMom)5 

(3.35) 

(3.36) 

=  &MLE- 


3.2.5  Normal  For  n  IID  data  samples,  the  Normal  pdf  is  defined  as  (Law  and 
Kelton,  1982:163) 

1  -1  " 

P(xi,x7,...,xn\a,n)  -  -^=exp(— ^(x,-/r)2),  (3.37) 

for  — oo  <  x  <  oo.  Since  the  same  parameters  were  determined  for  the  Lognormal  pdf, 
the  procedure  to  arrive  at  the  MLE  and  MoM  Normal  parameters  will  be  skipped  and  the 
results  presented: 


A  MLE 


°MLE 


PMoM 


®MoM 


n7^i 


(3.38) 

(3.39) 

(3.40) 

(3.41) 
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3.2.6  Beta  For  n  IID  data  samples,  the  Beta  pdf  is  defined  as  (Law  and  Kelton, 
1982:165) 


p(xi,x2,...,xn\a,0) 


Tn(a  + 13) 
Tn(a)Tn(f3) 


n 


n 


for  0<*i,  x2, ...,  x„<l.  Taking  the  natural  logarithm  of  both  sides  gives: 


(3.42) 


\np(xi,xi,  ...,xn\a, /3)  =  nlnr(c*  4-  (3)  —  nlnr(aQ3.43) 

n  n 

-n\nT((3)  +  (a  -  ^^lnx;  4-  (/3  -  l)^ln(l  -  x,).  (3-44) 

t=l  1=1 

Taking  the  derivative  with  respect  to  a  and  setting  the  result  equal  to  zero  gives: 


dlnp(zi,x2,  ...,x„\a,(3) 
da 


d  In  T(a  4-/3) 

”  da 


dlnT(a) 
U  da 


4-  J^lnxi 

«=i 


(3.45) 


=  0. 


Taking  the  derivative  with  respect  to  /3  and  setting  the  result  equal  to  zero  gives: 

(<>- 

d(3 


d\np(xi xn\a, (3)  =  n^In(^  +  P)  _  4-  ]Tln(l  -  xt)  (3.46) 


i= 1 


=  0. 


Beckman  and  Tietjen  show  that  equations  3.45  and  3.46  can  be  rearranged  to 

®,(a)-19?(d  4-/3)  =  lnGj 
(/3)  -  *(d  + /3)  =  lnG2, 

where  $  is  the  digamma  function  and 

g,  =  ntx,)* 

<=i 


(3.47) 

(3.48) 


(3.49) 
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(3.50) 


<5,  =  fl(l  -  *.)-■ 

1=1 

(Beckman  and  Tietjen,  1978:253-258).  They  proceed  to  show  that  equations  3.47  and  3.48 
can  be  solved  simultaneously,  resulting  in  one  equation  with  one  variable: 

¥(/3)  -  ^{(^(lnG,  -  lnG2  +  ¥(/3)))  + /3}  -  lnG2  =  0.  (3.51) 


They  give  a  Fortran  root-solving  program  which,  for  the  assumptions  they  make,  gives 
six  decimal  place  accuracy.  They  also  give  a  table  of  selected  Gx  and  G2  values  and 
corresponding  aMLE  and  Pmle  values.  Law  and  Kelton  expand  the  table  to  cover  the  full 
range  of  Gx  and  G2  values  (Law  and  Kelton,  1982:213-214).  Values  not  listed  in  the  tables 
can  be  arrived  at  by  interpolation.  Beckman  and  Tietjen  claim  two  decimal  place  accuracy 
for  this  method. 

The  MoM  estimates  for  a  and  /3  can  be  found  without  resorting  to  numerical  methods. 
The  mean  and  variance  of  the  Beta  pdf  are  given  as  (Law  and  Kelton,  1982:165): 


a 


2 


a  +  /3 


a/3 

(a  +  /3  +  l)(a  +  /3)2 
-»*• 


Solving  for  /z2  gives: 

a/3  _  a 

^  ~  (a +  /3+  l)(a  +  /3)2  +  a  + /3 
a(a  +1) 

(a  +  /3)(a  +  /3  + 1) 


(3.52) 

(3.53) 

(3.54) 

(3.55) 


(3.56) 

(3.57) 

(3.58) 
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Solving  equations  3.52  and  3.56  simultaneously  for  a  and  (3  gives: 


OtMoM 


PmoM 


&!*<)* 


(3.59) 

(3.60) 


3.2.7  Weibull  For  n  IID  data  samples,  the  Weibull  pdf  is  defined  as  (Law  and 
Kelton,  1982:163): 

n 

P(*i»*  2,...,x„|a,/?)  =  J^a/?"axf_1exp(-(^)"),  (3.61) 

<=i  P 

Xi  >  0.  (3.62) 


Taking  the  natural  logarithm  of  both  sides  gives: 


lnp(®1,X2,-,*n|o,/?) 


n  In  a  -  na  In  f3  +  (a  -  1)  ^3  lnx,  -  (3  a  Y'  xf 


i=l 


-.=1 


(3.63) 


Taking  the  derivative  with  respect  to  (3  gives: 


d\np(xuxi,  ...,x„\a/3) 

dp 


—na 

~r 


+  a^-(o+1)  53  xf . 
»  =  1 


Setting  the  result  equal  to  zero  and  solving  for  f)  gives: 


0«LE  =  (1 £>,*)*• 

71  ,  =  1 

Taking  the  derivative  with  respect  to  or  gives: 


(3.64) 


(3.65) 


d\np(xl,X2,...,xn\a,P) 

da 

-^““^iflnx;. 

i=l 


- n  In  (3  4-  ^3  In**  +  P  “ln/jV^x®  (3.66) 

°  i= 1  «=1 

(3.67) 
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Setting  the  result  equal  to  zero,  substituting  the  expression  for  (3  in  equation  3.65  into  (3 
in  equation  3.66  and  solving  for  a  gives: 


-  _  / Hi=i  1 

OtMLE  -  (— — d - • 

2^,=i  xi  «  i=1 

d  can  be  solved  for  numerically,  using  Newton’s  Method: 


.  =  x 
”  +  1  "  f’(*n)‘ 


The  resulting  equation  to  find  a  using  Newton’s  Method  is 


<**+1  —  Oik  + 


i  r-~.fi  ,  i  r"  x&k  in*. 

-  ,  Ini-  4-  — —  _!  , _ L 

n  2^,= 1  T  &k  yt-  Ii» 


1  —  t  * 


i  t.:~x 


-(Er=,*f‘inr')3 


(3.68) 


(3.69) 


(3.70) 


The  initial  value  chosen  for  d  can  be  the  MoM  value  which  will  be  determined  next. 
The  first  two  moments  of  the  Weibull  pdf  are  (Derman,  1973:386) 


Thus, 


M  =  13  r(i  +  — ) 

a 

M2  =  /?2r(i  +  -). 

a 


a 2  (r(l+^-))2 

—  —  v  v  QMoM;' 

n2  r(  i  +  £-£— ) 

'  a  MoM  ' 


MoM  ' 
2 


_  (iEkM 


I  F"  T? 
n  Z^>  =  1  Li 


(3.71) 

(3.72) 


(3.73) 

(3.74) 


Derman  gives  a  table  for  determining  values  of  £  given  (Derman,  1973:386).  Once 

d  is  found,  (3  is  given  by 


0M 


oM  — 


_  »  E"=l  Xi 


r(i  + 


)' 


(3.75) 
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3.3  Data  Inputs 

The  program  ASPECT  accepts  data  inputs  in  nxm  matrix  form.  This  is  meant  to 
allow  windows  of  data  to  be  studied.  Because  the  RCS  surface  map  is  3-dimensional,  care 
should  be  taken  when  determining  window  size  to  prevent  data  distortion  at  the  window 
edges.  A  6 j  by  4>\  size  data  set  can  be  compared  with  any  d2  by  <f>2  data  set,  but  both  data 
sets  must  result  in  rectangular  matrices.  Every  sample  interval  increment  along  the  0-axis 
should  be  a  new  row  and  every  increment  along  the  <£-axis  should  be  a  new  column. 

ASPECT  is  designed  to  build  pdf  models  from  static  data  which  are  assumed  to  be 
Independent  Identically  Distributed  (IID)  random  variables.  Dynamic  data  may  be  used  if 
no  significant  correlation  exists  between  the  time  samples.  If  significant  correlation  exists, 
the  time  spacing  can  be  lengthened  to  a  point  where  there  is  a  maximum  acceptable 
correlation  between  the  time  samples.  The  resulting  data  should  be  checked  to  insure 
enough  samples  have  been  used  to  represent  a  statistically  significant  data  base  (Dowdy, 
1991:165). 

Past  efforts  toward  computerized  RCS  modeling  have  focused  on  building  models 
assuming  all  aspect  angles  have  an  equal  probability  of  being  presented  to  a  radar.  In 
these  cases,  data  from  a  narrow  range  of  aspect  angles  were  compared  with  the  whole 
range  of  values  from  the  aircraft.  The  values  from  the  windows  tended  to  be  concentrated 
in  one  narrow  range  of  the  spectrum  of  the  entire  data  set.  This  concentration  would 
result  in  small  sample  variances  at  most  aspect  angles.  The  mean  of  the  window  would 
vary  relative  to  the  full  range  of  values;  its  size  depended  on  what  part  of  the  aircraft  was 
being  modeled.  In  contrast,  the  mean  of  the  window  relative  to  itself  would  stay  fairly 
close  to  the  median  of  the  window  and  the  variance  would  be  much  larger  than  the  variance 
of  the  window  compared  to  full-range  data.  Henceforth  in  this  study,  when  parameters 
are  estimated  by  comparing  a  localized  window  of  data  against  the  full-range  of  aircraft 
data,  the  comparison  will  be  referred  to  as  Window  vs.  Full-Range.  When  the  parameters 
are  estimated  by  comparing  a  localized  window  of  data  with  itself,  the  comparison  will 
be  referred  to  as  Window  vs.  Window.  The  stability  of  the  mean  and  large  size  of  the 
variance  for  the  Window  vs.  Window  data  are  caused  by  the  fact  that  the  normalization 
of  the  data  causes  the  windowed  data  to  cover  as  many  values  as  the  full-range  data.  A 
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consequence  of  the  different  statistical  natures  of  the  Window  vs.  Full-Range  and  the 
Window  vs.  Window  comparisons  is  that  there  are  differences  in  the  level  of  effectiveness 
pdfs  have  with  the  two  comparisons. 

3.4  ASPECT  Program  Operation 

This  section  gives  a  step  by  step  explanation  of  the  operation  of  the  ASPECT  pro¬ 
gram. 


3.4.I  Step  1:  Input  Prompts  The  user  is  prompted  to  input  the  name  of  the  data 
file  to  be  used,  the  number  of  separate  windows  the  full-range  of  data  is  to  be  taken  from, 
and  the  length  and  width  of  each  of  the  full-range  windows.  The  program  then  linearizes 
the  data,  turning  it  into  a  single  row  of  values.  The  data  is  linearized  one  row  at  a  time 
so  the  position  of  the  modeled  window  sampled  within  the  full  range  doesn’t  shift.  The 
user  must  determine  before  hand  the  beginning  and  ending  linearized  sample  numbers  of 
the  modeled  window.  This  is  easily  done  for  a  Window  vs.  Window  comparison,  since  the 
total  one-dimensional  length  of  the  full-range  is  output  as  the  variable  ’data’,  ’data’  would 
be  the  ending  sample  in  this  case. 

The  data  is  normalized  as  it  is  input.  Because  ASPECT  performs  pdf  selection  by 
analyzing  the  statistical  characteristics  of  a  single  set  of  input  data,  the  relative  amplitudes 
of  the  samples  within  the  set  are  the  only  data  needed  to  generate  the  pdfs.  For  this  reason, 
normalization  can  be  performed  without  fear  of  distorting  the  empirical  distribution  of  the 
data.  ASPECT  normalizes  the  data  because  all  pdfs  except  the  Normal  require  minimum 
sample  values  of  0,  and  the  Beta  pdf  requires  a  maximum  sample  value  of  at  most  1. 

After  normalization,  the  user  is  prompted  for  the  beginning  and  ending  samples  of 
the  modeled  window.  It  is  from  this  sample  set  that  all  the  statistics  used  in  the  pdfs  are 
determined. 

3-4-2  Step  2:  Generation  of  Statistics  The  mean  and  second  moment  are  the  only 
statistics  necessary  to  determine  the  parameters  of  the  pdfs. 
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3.4-3  Step  3:  Parameter  Determination  In  this  section,  the  Beta  and  Weibull  pa¬ 
rameters  are  determined.  aMoM  and  (5m <>m  are  determined  automatically.  The  subroutine 
‘interp’  is  called  to  determine  the  Weibull  parameters,  ‘interp’  contains  Table  3.1  from 
Derman,  which  gives  values  for  dl~M  from  values  of  ‘interp’  uses  linear  interpolation 
to  find  —  - — ,  then  finds  Qmom- 

Once  the  Beta  and  Weibull  MoM  parameters  are  determined,  ASPECT  asks  the  user 
if  the  MLE  paremters  are  desired  for  the  Beta  or  Weibull  pdfs  (MLE  parameters  equal 
MoM  parameters  for  all  other  pdfs).  If  Weibull  MLE  is  desired,  the  subroutine  ‘newton’ 
is  called,  ‘newton’  prompts  the  user  for  an  initial  value  for  olmle  (omom  is  recommended) 
and  performs  a  Newton’s  Method  iteration  process  on  equation  3.70  until  equation  3.68  is 
true.  If  an  initial  value  for  &mle  is  chosen  too  far  from  the  actual  &mle,  the  iterations 
could  run  away  to  infinity,  in  which  case  the  ASPECT  program  operation  will  pause  and 
inform  the  user  that  oo  has  been  encountered  and  that  rc’  must  be  entered  to  break  out  of 
the  program.  If  another  attempt  is  made,  a  good  starting  point  usually  turns  out  to  be  a 
value  closer  to  zero  than  aMoM- 

If  the  Beta  MLE  parameters  are  desired,  the  ‘bmle’  subroutine  is  called.  This  sub¬ 
routine  performs  a  double  linear  interpolation  to  find  aMLE  and  (3Mle  from  Table  5.12  in 
Law  and  Kelton,  given  Gx  and  G2  as  described  earlier  in  this  chapter. 

If  the  MLE  parameters  are  produced,  they  will  be  the  parameters  used  in  the  cdf 
and  pdf  plots  and  the  goodness-of-fit  tests. 

3.4.4  Step  4:  PDF  and  CDF  Plots  After  all  pdf  parameters  are  determined,  the 
pdfs  are  plotted  individually,  as  is  a  histogram  of  the  windowed  data.  The  program  can 
be  manipulated  to  allow  multiple  graph  plots  for  visual  comparisons,  but  since  there  are 
many  occasions  when  one  pdf  or  histogram  will  dwarf  all  the  other  functions,  multiple 
graph  plots  are  not  done  automatically.  The  points  on  the  plots  are  found  by  finding 
values  for  the  pdfs  for  every  ^  increment  from  0  to  1. 

The  empirical  cdf  (ecdf)  is  found  by  first  comparing  the  windowed  data  with  the 
same  set  of  incremented  values  used  to  find  the  pdfs.  If  a  windowed  sample  is  within 

of  the  incremented  value,  a  counter  that  is  initialized  at  0  is  increased  by  1.  Each 
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windowed  data  sample  is  compared  with  each  incremented  value.  The  counter  is  zeroed  for 
each  new  incremented  value.  The  number  the  counter  assigns  to  each  incremented  value 
is  the  magnitude  of  the  increase  of  the  ecdf  at  the  incremented  value. 

The  other  cdfs  are  found  by  increasing  the  value  of  the  cdf  by  the  magnitude  of  the 
pdf  at  each  of  the  incremented  values  used  in  the  pdf  formation.  The  ecdf  and  the  cdfs 
are  then  all  normalized  by  dividing  each  cdf  by  its  largest  (final)  value. 

Each  cdf  is  then  individually  plotted  against  the  ecdf  for  a  visual  goodness-of-fit 
assessment.  As  with  the  pdf  plots,  the  program  can  be  manipulated  to  compare  more  than 
one  cdf  at  a  time. 

3-4-5  Step  5:  Goodness-of-Fit  Kolmogorov  and  Chi-Square  test  statistics  are  pro¬ 
duced  for  each  cdf.  The  Kolmogorov  test  statistics  are  found  as  explained  in  Chapter  2. 
To  determine  the  Chi-Square  test  statistics,  equal  probability  bin  intervals  must  be  found. 
Determining  the  interval  points  requires  finding  the  inverse  values  of  the  distribution  func¬ 
tions  for  the  required  bin  probabilities.  For  many  distributions,  the  inverse  values  cannot 
be  found  in  closed  form.  However,  the  approximate  interval  points  can  be  found  relatively 
easily  using  numerical  methods  as  demonstrated  in  the  following  example: 

Suppose  the  number  of  histogram  bins  chosen  was  three.  Then,  for  equal  probability, 
every  bin  would  have  a  33%  chance  of  a  sample  occurring  in  its  interval.  For  a  sample  set 
of  normalized  data,  the  first  bin  would  have  a  low  interval  of  0.  The  high  interval  point 
of  the  first  bin  and  the  low  interval  point  of  the  second  bin  would  be  the  input  value  that 
produced  the  cdf  value  closest  to  The  high  interval  point  of  the  second  bin  and  the  low 
interval  point  of  the  third  bin  would  be  the  input  value  that  produced  the  cdf  value  closest 
to  |.  The  high  end  interval  of  the  third  bin  would  be  1.  Values  from  the  window  would 
then  be  placed  in  the  bin  who’s  intervals  they  fell  between.  This  method  results  in  two 
decimal  place  accuracy. 

The  test  statistics  for  each  cdf  are  printed  out  for  the  Kolmogorov  and  the  Chi-Square 

tests. 
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3.4.6  Step  6:  Optimum  Model  Determination  The  optimum  pdf  is  declared  to  be 
the  pdf  with  the  smallest  Kolmogorov  test  statistic.  The  name  of  the  pdf  as  well  as  the 
estimated  parameters  are  printed  out.  The  Kolmogorov  test  statistic  is  used  because,  on 
occasions  when  the  Kolmogorov  test  and  the  Chi-Square  test  disagree,  visual  assessment 
usually  agrees  with  Kolmogorov. 

The  user  must  refer  to  tables  to  determine  the  level  of  significance  the  pdf  produces. 
3.5  Summary 

This  chapter  explained  the  computations  required  to  impliment  the  ‘ASPECT’  RCS 
model  optimization  program.  The  MLE  and  MoM  parameters  for  all  pdfs  used  were 
determined.  The  required  inputs  and  format  were  explained. 
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IV.  Determination  of  Optimum  RCS  Models  for  Window  vs.  Full-Range 

Sample  Sets 

4-1  Introduction 

This  chapter  will  begin  with  an  overview  of  the  types  of  parameters  of  several  pdfs 
and  the  effects  changes  in  those  parameters  will  have  on  the  shape  of  the  pdfs.  The  rest  of 
the  chapter  will  analyze  how  the  pdfs  handled  the  empirical  distribution  of  the  data  from 
each  aspect  window  when  compared  with  the  full-range  of  data.  Window  vs.  Full-Range 
comparisons  are  the  standard  method  for  forming  RCS  models. 

Three  different  types  of  data  will  be  tested  at  each  of  the  three  modeled  aspect 
windows.  The  first  type  will  be  fighter  data  with  samples  taken  in  the  0  =  90°  plane.  The 
second  type  will  be  fighter  data  with  samples  taken  in  a  wider  0  window:  0  =  90°±5°.  The 
third  type  of  data  will  be  missile  data  in  the  0  =  90°  plane. 

4-2  Parameter  Types 

The  pdfs  considered  in  this  investigation  will  be  characterized  by  three  parameters, 
namely  location,  scale,  and  shape. 

Location  parameters  determine  where  the  centroid  of  the  probability  densities  is 
located  in  the  range  of  the  pdf  values.  Location  parameters  are  used  in  the  Normal  and 
Lognormal  pdfs. 

Scale  parameters  compress  or  expand  the  probability  density  a  pdf  will  have  over 
its  range.  For  the  pdfs  used  here,  a  decreasing  scale  parameter  value  will  compress  pdfs, 
an  increasing  scale  parameter  value  will  expand  pdfs.  A  compressed  pdf  will  exhibit  very 
small  probability  densities  over  most  of  its  range  and  very  large  probability  densities  over 
a  small  portion  of  its  range.  An  expanded  pdf  will  exhibit  increased  probability  densities 
over  the  majority  of  its  range,  but  the  peak  probability  densities  will  be  much  lower  than 
those  of  a  compressed  pdf.  Scale  parameters  are  used  in  the  Normal,  Lognormal,  and 
Weibull  pdfs. 
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Shape  parameters  determine  both  the  shape  and  location  of  the  pdf.  They  can  cause 
the  pdf  to  behave  like  a  negative  or  positive  exponential  function,  a  lognormal  function,  or 
a  normal  function.  Shape  parameters  are  used  in  Beta,  Weibull,  and  Rayleigh-class  pdfs. 

All  parameters  used  in  this  study  are  linear  combinations  of  the  mean  (m)  and  second 
moment  (m2).  The  mean  and  second  moment  are  two  of  the  three  components  that  can 
determine  a  pdfs  effectiveness.  The  third  component  is  the  skew  of  the  pdf.  Skew  is  the 
difference  in  the  slopes  on  either  side  of  the  peak  values  of  a  pdf.  If  the  right  slope  is 
steeper  than  the  left  slope,  the  pdf  is  skewed  to  the  right.  If  the  left  slope  is  steeper  than 
the  right  slope,  the  pdf  is  skewed  to  the  left. 

The  following  subsections  discuss  each  pdf’s  response  to  changes  in  its  parameters. 

4-2.1  Normal  The  location  parameter  n  will  center  the  pdf  at  the  mean  of  the 
sample  set.  The  scale  paremeter  a  will  expand  the  pdf  for  increasing  values  of  variance 
and  compress  the  pdf  for  decreasing  values  of  variance.  The  Normal  pdf  is  never  skewed. 
Figure  4.1  shows  the  Normal  pdf’s  response  to  changes  in  (j.  and  o. 


Figure  4.1.  Normal  pdf’s  response  to  changes  in  fi  and  0 
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4-2.2  Lognormal  Since  the  lognormal  pdf  acts  as  a  normal  pdf  applied  to  the  nat¬ 
ural  logarithm  of  the  sample  set,  the  location  parameter  \x  centers  the  pdf  around  the 
mean  of  the  natural  logarithm  of  the  sample  set.  This  will  cause  the  pdf  to  be  skewed  to 
the  left  of  the  true  mean  of  the  sample  set.  The  skewing  will  be  further  enhanced  by  the 
sample  variable  in  the  denominator  of  the  non- exponential  portion  of  the  pdf  which  is  the 
result  of  the  Jacobian  of  the  transformation  y  =  In  x.  This  variable  will  result  in  enhanced 
probability  densities  for  lower-range  sample  values.  The  scale  parameter  n  compresses  the 
pdf  for  decreasing  sample  variances  and  expands  the  pdf  for  increasing  sample  variances. 
Figure  4.2  shows  the  Lognormal  pdf’s  response  to  changes  in  fi  and  o. 

4-2.3  Weibull  The  shape  parameter  a  will  cause  the  pdf  to  approximate  an  expo¬ 
nential  pdf  at  low  values  (a  <  3)  and  is  the  exponential  pdf  for  a  =  1.  At  higher  values  of 
a,  the  pdf  will  become  more  normal  in  appearance  and  center  on  the  sample  mean.  The 
scale  parameter  (3  compresses  the  pdf  for  small  values  of  (3  and  expands  the  pdf  for  large 
values  of  (3.  The  Weibull  pdf  is  always  skewed  to  the  right  for  values  of  a>4.  Figure  4.3 
shows  the  Weibull  pdf’s  response  to  changes  in  a  and  fi. 

4-2-4  Beta  The  Beta  pdf  has  two  shape  parameters,  a  and  f3.  The  ratio  and 
magnitudes  of  a  and  (3  determine  the  shape  of  the  pdf.  If  a  >  fi,  the  pdf  will  be  skewed  to 
the  right  against  the  high  boundary  of  the  data  range.  If  a  <  (3,  the  pdf  will  be  skewed  to 
the  left  against  the  low  boundary  of  the  data  range.  If  a  =  fi,  the  pdf  will  be  normal  shaped 
and  will  be  centered  in  the  data  range.  If  a  =  fi  =  1,  the  pdf  is  uniform,  with  a  magnitude 
of  1.  If  a  or  fi  is  less  than  1,  the  pdf  approximates  an  exponential  function,  increasing 
toward  the  high  range  boundary  if  a  >  fi  or  the  low  range  boundary  if  a  <  fi.  If  a  and 
fi  are  both  less  than  one,  then  the  pdf  will  approximate  two  exponential  functions.  One 
exponential  will  increase  toward  the  high  range  boundary,  the  other  will  increase  toward 
the  low  range  boundary.  Figure  4.4  shows  the  Beta  pdf’s  response  to  changes  in  a  and  f3. 

The  Beta  pdf  is  the  most  flexible  pdf  used  in  this  study,  making  it  useful  for  sample 
sets  that  do  not  exhibit  a  Normal-type  distribution.  However,  it  is  much  less  responsive  to 
changes  in  variance  than  the  Lognormal,  Normal,  or  Weibull  pdfs  and  cannot  attain  the 
highly  compressed  forms  required  to  model  sample  sets  with  small  variances. 
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Figure  4.2.  Lognormal  pdf’s  response  to  changes  in  n  and  a 


Figure  4.3.  Weibull  pdf’s  response  to  changes  in  a  and  (3 


Figure  4.4.  Beta  pdf’s  response  to  changes  in  a  and  (3 


4-2.5  Rayleigh  Class  The  Rayleigh  class  shape  parameter  o  is  a  linear  transforma¬ 
tion  of  the  second  moment  of  the  sample  set.  The  maximum  values  of  the  Rayleigh  class 
pdfs  are  always  skewed  to  the  left  of  the  sample  mean.  Small  values  of  a  will  compress  the 
pdfs  and  force  the  maximum  probability  density  closer  to  the  sample  mean.  Large  values 
of  a  will  expand  the  pdfs  and  force  the  maximum  probability  density  farther  away  from 
the  sample  mean.  The  Rayleigh  (y  =  x2)  pdf  is  actually  a  negative  exponential  function, 
whose  highest  probability  density  is  always  at  zero.  Because  the  Rayleigh  class  pdfs  are 
dependent  only  on  the  sample  second  moment,  they  are  less  responsive  to  changes  in  vari¬ 
ance  than  any  of  the  other  pdfs  used  in  this  study.  Figure  4.5  shows  the  Rayleigh  pdf’s 
responses  to  changes  in  a. 

4-3  Test  Results  for  Window  vs.  Full-Range  RCS  Models 

This  section  will  present  an  analysis  of  the  behaviour  of  the  pdfs  to  different  aspect 
windows  studied.  The  Kolmogorov  Level  of  Significance,  and  the  Chi-Square  Level  of 
Significance  for  Window  vs.  Full-Range  RCS  Models  are  presented  in  Appendix  A. 
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Figure  4.5.  Rayleigh  pdf’s  response  to  changes  in  a 

Due  to  the  Rayleigh-class  pdfs’  inability  to  compress  for  low  variances,  the  Rayleigh- 
class  pdfs  (Rayleigh,  One-Dominant- Plus- Rayleigh,  and  Rayleigh  Squared  made  poor  Win¬ 
dow  vs.  Full-Range  models  for  all  sample  sets.  For  this  reason,  the  performance  of  the 
Rayleigh-class  pdfs  will  not  be  discussed  in  the  following  sections,  but  the  results  are  listed 
in  the  tables  of  Appendix  A  along  with  the  other  pdfs. 

In  the  following  sections,  the  pdfs  with  the  lowest  Kolmogorov  and  Chi-Square  test 
statistics  are  presented.  If  more  than  one  pdf  has  the  lowest  Chi-Square  test  statistic,  the 
low  Chi-Square  pdf  with  the  lowest  Kolmogorov  test  statistic  is  presented  as  the  best  Chi- 
Square  pdf.  If  that  pdf  happens  to  be  the  pdf  with  the  lowest  Kolmogorov  test  statistic  of 
all  pdfs,  two  Chi-Square  pdfs  are  presented  for  comparison. 

4-3.1  Nose-On  Fighter:  6  =  90°  The  sample  set  is  characterized  by  a  medium 
range  mean  (m=0.52),  a  relatively  small  variance  (var=7.06  x  10-4),  a  small  number  of 
samples  (n=ll),  and  is  skewed  to  the  right.  The  cdfs  and  pdfs  with  the  lowest  Kolmogorov 
and  Chi-Square  test  statistics  are  compared  with  the  sample  set  and  presented  in  Figures 
4.6-4.9. 
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The  Kolmogorov  Goodness-of-Fit  Test  (KGFT)  declares  the  Weibull  MLE  and  MoM, 
Lognormal,  and  Normal  pdfs  to  have  a  Level  of  Significance  (LoS)  greater  than  20%  (the 
highest  LoS  given  by  Conover  (Conover,  1980:  Table  14).  However,  the  Weibull  MLE  has 
the  smallest  test  statistic  and  is  chosen  as  the  optimum  pdf.  The  Weibull  MLE  is  the  most 
compressed  pdf  (highest  probability  density  over  narrowest  range).  As  can  be  seen  by  the 
CDF  figures,  the  Empirical  Cumulative  Distribution  Function  (ECDF)  is  skewed  to  the 
right.  The  Weibull  pdfs  are  the  only  pdfs  able  to  match  this  skewing  to  any  degree. 

Because  of  the  small  number  of  samples,  the  Chi-Square  Goodness-of-Fit  Test  (CGFT) 
is  limited  to  one  degree  of  freedom  by  the  restrictions  listed  in  section  2.5.  For  this  rea¬ 
son,  the  miminum  test  statistic  is  shared  by  the  Weibull  MLE  and  MoM,  Lognormal,  and 
Normal  pdfs.  The  CGFT  is  more  conservative  than  the  KGFT,  declaring  the  LoS  to  be 
between  10%  and  25%.  This  range  of  LoS  is  the  second  highest  range  listed  by  Conover 
(Conover,  1980:  Table  2). 

4-3.2  Nose-On  Fighter:  6  =  90°±5°  This  sample  set  is  characterized  by  a  medium 
range  mean  (m=0.45),  a  small  variance  (var=2.21xl0~4),  a  large  number  of  samples 
(n=121),  and  a  rightward  skew.  The  pdfs  and  cdfs  with  the  lowest  Kolmogorov  and  Chi- 
Square  test  statistics  are  presented  in  Figures  4.10-4.13.  The  KGFT  declares  the  Weibull 
MLE  and  MoM  pdfs  to  have  a  LoS  greater  than  20%,  with  the  Weibull  MoM  having  a 
slightly  lower  test  statistic  than  the  Weibull  MLE.  The  Weibull  MoM  matches  the  test 
statistic  with  the  highest  compression.  The  Weibull  pdfs  are  the  only  ones  able  to  match 
the  rightward  skew  of  the  sample  set. 

The  CGFT  declares  the  Normal  pdf  to  have  the  highest  LoS  (10%  >  a  >  5%).  The 
CGFT  declares  all  the  other  pdfs  to  have  a  LoS  less  than  0.1%.  This  is  a  poor  showing 
for  all  pdfs  and  is  probably  due  to  the  few  degrees  of  freedom  (2)  used  in  the  test.  Time 
limitations  prevented  building  a  CGFT  with  the  many  degrees  of  freedom  which  can  be 
used  with  large  sample  sets.  Visual  assessments  of  Figures  4.10-4.13  show  the  Weibull 
MoM  pdf  to  be  virtually  indistinguishable  from  the  sample  set,  while  the  Normal  pdf 
cannot  match  the  rightward  skew  of  the  sample  set. 
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4-3.3  Nose-On  Missile:  6  =  90°  This  sample  set  is  characterized  by  a  medium 
range  mean  (m=0.47),  a  small  variance  (var=2.23  x  10-4),  a  medium  number  of  samples 
(n=21),  and  no  visible  skew.  The  best  Kolmogorov  and  ChiSquare  pdfs  and  cdfs  are 
presented  in  Figures  4.14-4.17. 

The  KGFT  declares  the  Weibull  MLE  and  MoM,  Lognormal,  and  Normal  pdfs  to 
have  a  LoS  greater  than  20%.  The  Normal  pdf  has  the  smallest  test  statistic,  the  Lognor¬ 
mal’s  test  statistic  is  slightly  larger.  Unlike  the  other  Nose-On  sample  sets,  the  optimum 
pdf  does  not  have  the  highest  probability  distribution.  In  fact,  the  Weibull  MLE  and  MoM 
and  Lognormal  pdfs  all  have  higher  probability  densities  than  the  Normal.  However,  the 
Weibull  pdfs  are  skewed  to  the  right  and  the  Lognormal  pdf  is  skewed  to  the  left.  The 
Normal  pdf  is  optimum  because  of  the  lack  of  a  skewed  sample  set  and  the  closeness  of  the 
values  of  the  peak  probability  densities. 

The  CGFT  declares  the  Weibull  MLE  and  MoM,  Lognormal,  and  Normal  pdfs  all  to 
have  a  LoS  between  2.5%  and  5%.  Though  a  visual  assessment  of  Figures  4.14  and  4.16 
shows  the  Lognormal  and  Normal  pdfs  to  be  a  very  close  match  to  the  sample  set,  the 
CDFT  declares  all  the  pdfs  to  be  a  poor  fit. 
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Nose-On  RCS  vs.  Window-Full  Range  Weibull  MLE  CDF 

Figure  4.6.  Top  KGFT  and  CGFT  CDF  of  Fighter  Data  at  Nose-On  ±5°,  0.5°  sample 
interval,  6  =  90°,  6-0  polarization,  1  GHz  Full-Scale 


Nose-On  RCS  vs.  Window-Full  Range  Weibull  MLE  PDF 


Figure  4.7.  Top  KGFT  and  CGFT  PDF  of  Fighter  Data  at  Nose-On  ±5°,  0.5°  sample 
interval,  0  =  90°,  6-0  polarization,  1  GHz  Full-Scale 
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Figure  4.8.  Top  CGFT  CDF  of  Fighter  Data  at  Nose-On  ±5°,  0.5°  sample  interval,  6 
90°,  6  —  6  polarization,  1  GHz  Full-Scale 
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Figure  4.9.  Top  CGFT  PDF  of  Fighter  Data  at  Nose-On  ±5°,  0.5°  sample  interval,  6 
90°,  6  —  6  polarization,  1  GHz  Full-Scale 
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Nose-On  RCS/5  vs.  Window-Full  Range  Weibull  MoM  CDF 


Figure  4.10.  Top  KGFT  CDF  of  Fighter  Data  at  Nose-On  ±5°,  0.5°  sample  interval, 
9  =  90°  ±  5°,  9  -  9  polarization,  1  GHz  Full-Scale 


Nose-On  RCS/5  vs.  Window-Full  Range  Weibull  MoM  PDF 


Figure  4.11.  Top  KGFT  PDF  of  Fighter  Data  at  Nose-On  ±5°,  0.5°  sample  interval, 
9  -  90°  ±  5°,  9  -  0  polarization,  1  GHz  Full-Scale 
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Figure  4.12.  Top  CGFT  CDF  of  Fighter  Data  at  Nose-On  ±5°,  0.5°  sample  interval, 
9  =  90°  ±  5°,  9  -  8  polarization,  1  GHz  Full-Scale 


Figure  4.13.  Top  CGFT  PDF  of  Fighter  Data  at  Nose-On  ±5®,  0.5°  sample  interval, 
9  =  90®  ±  5®,  8  -  9  polarization,  1  GHz  Full-Scale 
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P(x) 


Figure  4.14. 


P(x) 


Figure  4.15. 


Top  KGFT  and  CGFT  CDF  of  Missile  Data  at  Nose-On  ±5*,  0.5°  sample 
interval,  6  =  90°,  6-6  polarization,  18  GHz  Full-Scale 
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Top  KGFT  and  CGFT  PDF  of  Missile  Data  at  Nose-On  ±5°,  0.5®  sample 
interval,  6  =  90°,  6  —  6  polarization,  18  GHz  Full-Scale 
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Figure  4.16.  Top  CGFT  CDF  of  Missile  Data  at  Nose-On  ±5°,  0.5®  sample  interval, 
0  =  90°,  0  —  9  polarization,  18  GHz  Full-Scale 
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Figure  4.17.  Top  CGFT  PDF  of  Missile  Data  at  Nose-On  ±5°,  0.5°  sample  interval, 
0  =  90°,  9  —  0  polarization,  18  GHz  Full-Scale 
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4-3-4  Broadside  Fighter:  0  =  90°  This  sample  set  is  characterized  by  a  high  range 
mean  (m=0.67),  a  large  variance  (var=1.44  x  10~ 2),  a  medium  number  of  samples  (n=21), 
and  a  slight  rightward  skew.  The  best  Kolmogorov  and  Chi-Square  pdfs  and  cdfs  are 
presented  in  Figures  4.18-4.21. 

The  KGFT  declares  the  Beta  MLE  and  MoM,  Weibull  MLE  and  MoM,  and  Normal 
pdfs  to  have  a  LoS  greater  than  20%.  The  Beta  MLE  has  the  lowest  test  statistic.  Since 
the  sample  mean  is  greater  than  0.5,  the  Beta  pdf  is  skewed  rightward,  along  with  the 
Weibull  pdf.  The  large  variance  greatly  expands  the  variance  sensitive  pdfs  (Lognormal, 
Normal,  and  Weibull),  leaving  the  Beta  MLE  pdf  with  the  highest  probability  density. 

The  CGFT  agrees  with  the  KGFT,  declaring  the  Beta  MLE  and  Weibull  MLE  pdfs 
to  have  a  LoS  greater  than  25%.  However,  the  CGFT  rates  the  Beta  MoM,  Weibull  MoM, 
Lognormal,  and  Normal  pdfs  relatively  poorly  (10%  >  LoS  >  2.5%). 

4-3.5  Broadside  Fighter:  6  =  90°  ±5°  This  sample  set  is  characterized  by  a  medium 
range  mean  (m=0.55),  a  medium  sized  variance  (var=4.3x  10~3),  a  large  number  of  samples 
(n=231),  and  a  slight  rightward  skew.  The  resulting  pdfs  and  cdfs  are  presented  in  Figures 
4.22  and  4.23. 

The  large  percentage  of  off-broadside  samples  in  this  sample  set  compared  with  the 
Broadside  Fighter:#  =  90°  sample  set  has  lowered  the  sample  mean.  The  medium  range 
sample  mean  prevents  the  Beta  pdf  from  skewing.  The  medium  size  variance  expands  the 
Normal  and  Lognormal  pdfs  more  severely  than  the  Weibull  pdfs.  The  KGFT  declares  the 
Weibull  MLE  and  MoM  pdfs  to  have  a  LoS  greater  than  20%.  The  Weibull  MoM  has  the 
lowest  test  statistic  and  the  highest  probability  density. 

The  CGFT  gives  the  Weibull  MoM  the  lowest  test  statistic,  which  results  in  the 
highest  LoS  (10%  >  a  >  5%),  a  relatively  poor  result. 

4-3.6  Broadside  Missile:  9  =  90°  This  sample  set  is  characterized  by  a  high  range 
mean  (m=0.81),  a  large  variance  (var=1.96  x  10~2),  a  small  number  of  samples  (n=41), 
and  a  strong  rightward  skewing.  The  resulting  pdfs  and  cdfs  are  presented  in  Figures 
4.24-4.27. 
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The  KGFT  declares  the  Beta  MLE  and  MoM,  the  Weibull  MLE  and  MoM,  and  the 
Normal  pdfs  to  have  a  LoS  greater  than  20%.  However,  the  combination  of  a  wide  variance, 
a  high  mean,  and  a  strong  rightward  skew  gives  the  Beta  pdfs  much  lower  test  statistics 
than  the  other  pdfs.  The  Beta  MLE  has  the  highest  probability  density,  but  this  level  of 
compression  forces  the  peak  probability  density  far  to  the  right  of  the  sample  mean.  With 
the  /3  parameter  approaching  1.00,  the  Beta  MLE  pdf  approximates  a  positive  exponential 
function,  the  Beta  MoM  does  not  compress  as  severely  as  the  Beta  MoM  with  the  large 
sample  mean  and  therefore  attains  the  lowest  test  statistic.  Although  the  Weibull  pdfs  are 
skewed  to  the  right  as  well,  the  wide  variance  expands  the  pdfs  too  severely  for  them  to 
be  as  effective  as  the  Beta  pdfs. 

The  CGFT  declares  the  Beta  MLE  and  MoM  pdfs  to  have  a  LoS  greater  than  25%. 
The  CGFT  rated  the  Weibull  MLE  and  MoM  relatively  highly  as  well,  giving  them  a  LoS 
between  10%  and  25%. 
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Figure  4.18.  Top  KGFT  and  CGFT  CDF  of  Fighter  Data  at  Broadside  ±5°,  0.5°  sample 
interval,  0  =  90°,  0  —  0  polarization,  1  GHz  Full-Scale 


Figure  4.19.  Top  KGFT  and  CGFT  PDF  of  Fighter  Data  at  Broadside  ±5°,  0.5°  sample 
interval,  0  =  90°,  0-0  polarization,  1  GHz  Full-Scale 


4-17 


Broadside  RCS  vs.  Window-Pull  Range  Weibull  MLE  CDF 


Figure  4.20.  Top  CGFT  CDF  of  Fighter  Data  at  Broadside  ±5°,  0.5°  sample  interval, 
0  =  90°,  0-0  polarization,  1  GHz  Full-Scale 
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Figure  4.21.  Top  CGFT  PDF  of  Fighter  Data  at  Broadside  ±5°,  0.5°  sample  interval, 
9  =  90°,  0-0  polarization,  1  GHz  Full-Scale 


4-18 


1 


Broadside  RCS/5  vs.  Window-Full  Range  Weibull  MoM  PDF 

Figure  4.22.  Top  KGFT  and  CGFT  CDF  of  Fighter  Data  at  Nose-On  ±5®,  0.5°  sample 
interval,  9  =  90°  ±  5°,  0  -  9  polarization,  1  GHz  Full-Scale 


Broadside  RCS/5  vs.  Window-Full  Range  Weibull  MoM  PDF 

Figure  4.23.  Top  KGFT  and  CGFT  PDF  of  Fighter  Data  at  Nose-On  ±5°,  0.5°  sample 
interval,  0  =  90®  ±  5°,  9  -  0  polarization,  1  GHz  Full-Scale 


4-19 


1 

0.8 

0.6 

POO 

0.4 

0.2 

0 

0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

Broadside  RCS  vs.  Window-Full  Range  Beta  MoM  CDF 

Figure  4.24.  Top  KGFT  and  CGFT  CDF  of  Missile  Data  at  Broadside  ±5°,  0.5°  sample 
interval,  0  =  90°,  0-0  polarization,  18  GHz  Full-Scale 
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Figure  4.25.  Top  KGFT  and  CGFT  PDF  of  Missile  Data  at  Broadside  ±5°,  0.5°  sample 
interval,  0  =  90°,  0  —  0  polarization,  18  GHz  Full-Scale 
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Figure  4.26.  Top  CGFT  CDF  of  Missile  Data  at  Broadside  ±5",  0.5*  sample  interval, 
9  —  90°,  0  —  6  polarization,  18  GHz  Full-Scale 
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Figure  4.27.  Top  CGFT  PDF  of  Missile  Data  at  Broadside  ±5°,  0.5°  sample  interval, 
9  =  90°,  9-9  polarization,  18  GHz  Full-Scale 
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4-3.7  Tail-On  Fighter:  0  =  90°  This  sample  set  is  characterized  by  a  low  range 
mean  (m=0.23),  a  medium  sized  variance  (var=2.8  x  10-3),  a  small  number  of  samples 
(n=ll),  and  a  slight  skew  to  the  right.  The  resulting  pdfs  and  cdfs  are  presented  in  Figures 
4.28-4.31.  The  low-range  mean  skews  the  Beta  pdfs  to  the  left,  approaching  the  Lognormal 
function  in  shape.  The  medium  size  variance  expands  the  Lognorma  nd  Normal  pdfs  to 
a  greater  extent  than  the  Weibull  pdfs.  The  KGFT  declares  the  Beta  MLE  and  MoM, 
Weibull  MLE  and  MoM,  Lognormal,  and  Normal  pdfs  to  have  a  LoS  greater  than  20%.  The 
Weibull  MLE  pdf  matches  the  rightward  skew  with  the  highest  compression  and  attains 
the  highest  probability  density  and  lowest  test  statistic. 

The  CGFT  declares  the  Beta  MLE  and  MoM,  Weibull  MLE  and  MoM,  and  Normal 
pdfs  to  have  a  LoS  greater  than  25%.  The  one  degree  of  freedom  prevents  picking  a  winner 
from  these  pdfs. 

4-3.8  Tail- On  Fighter:  6  =  90°  ±  5®  This  sample  set  is  nearly  identical  to  the 
Tail-On  Fighter:  0  =  90°  sample  set.  It  is  characterized  by  a  low-range  mean  (m=0.25),  a 
medium  size  variance  (var=2.5x  10-3),  a  large  number  of  samples  (n=121),  and  a  rightward 
skew.  The  resulting  pdfs  and  cdfs  are  presented  in  Figures  4.32-4.35. 

As  in  the  last  section,  the  Beta  pdfs  are  skewed  to  the  left  along  with  the  Lognormal 
pdf,  and  the  Lognormal  and  Normal  pdfs  are  expanded  to  a  greater  extent  than  the  Weibull 
pdfs.  The  Weibull  pdfs  match  the  rightward  skew  of  the  sample  set.  The  KGFT  declares 
the  Weibull  MLE  and  MoM  and  the  Normal  pdfs  to  have  a  LoS  greater  than  20%.  The 
Weibull  MLE  attains  the  lowest  test  statistic  and  highest  probability  density. 

The  CGFT  declares  the  Beta  MoM,  Weibull  MoM,  and  Normal  pdfs  to  have  a  LoS 
between  5%  and  2.5%.  This  is  another  poor  performance  for  a  large  sample  set.  The  Beta 
MoM  and  Normal  pdfs  have  the  lowest  test  statistic. 

4-3.9  Tail-On  Missile:  6  =  90°  This  sample  set  is  characterized  by  a  medium  range 
mean  (m=0.58),  a  small  variance  (var=2.82  x  10-4),  a  small  number  of  samples  (n=21), 
and  a  slight  leftward  skew.  The  pdfs  and  cdfs  are  presented  in  Figures  4.36-4.39. 
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This  is  the  only  sample  set  that  is  skewed  to  the  left.  With  the  middle  range  mean 
forcing  the  Beta  pdf  into  a  Normal- type  distribution,  the  Lognormal  pdf  is  the  only  pdf 
that  can  match  the  leftward  skew.  In  addition,  the  small  variance  gives  the  Lognormal  and 
Normal  pdfs  the  compression  to  match  that  of  the  Weibull  pdfs.  The  KGFT  declares  the 
Weibull  MLE  and  MoM,  Lognormal,  and  Normal  pdfs  to  have  a  LoS  greater  than  20%. 
The  Lognormal  pdf  matches  the  highest  probability  density  with  the  lowest  test  statistic. 

The  CGFT  declares  the  Weibull  MLE  and  MoM,  Lognormal,  and  Normal  pdfs  to 
have  a  LoS  between  10%  and  25%.  These  pdfs  have  identical  test  statistics. 

4-4  Summary  of  Window  vs.  Full-Range  RCS  Models 

Statistical  analysis  of  the  fighter  and  missile  data  reveals  that  even  at  the  largest 
Window  vs.  Full-Range  variances,  high  compression  is  required  to  effectively  model  the 
sample  sets.  The  most  compressed  pdf  attains  the  lowest  Kolmogorov  test  statistic  in  seven 
of  the  nine  sample  sets  modeled.  Because  of  the  importance  of  maximum  compression, 
the  pdfs’  response  to  changes  in  variance  is  the  most  important  factor  in  determining 
best  goodness-of-fit.  The  Beta  pdfs’  limited  response  to  variance  changes  make  the  pdfs 
ineffective  for  small  variance  sample  sets  but  most  effective  for  large  variance  sample  sets. 
The  reverse  is  true  for  the  Lognormal  and  Normal  pdfs.  The  Weibull  pdfs  compress  at 
small  variances  at  the  same  rate  as  the  Lognormal  and  Normal  pdfs  but  do  not  expand 
at  the  same  rate  as  the  Lognormal  and  Normal  pdfs  for  large  variances.  Therefore,  the 
Weibull  pdfs  are  able  to  perform  effectively  at  both  high  and  low  variance  sample  sets. 
The  Weibull’s  response  to  variance  changes  enables  the  pdf  to  be  the  only  function  that 
attains  a  Kolmogorov  LoS  greater  than  20%  for  every  sample  set. 

The  skew  of  the  sample  sets  is  also  important  in  determining  goodness-of-fit.  In  seven 
of  the  nine  sample  sets  modeled,  the  sample  set  is  skewed  to  the  right.  This  enhances  the 
Weibull  pdfs’  effectiveness  and  hurt  the  Lognormal  and  Normal  pdfs’  effectiveness.  It  helps 
the  Beta  pdfs’  effectiveness  for  high-mean  sample  sets,  but  hurt  the  Beta  pdfs’  effectiveness 
elsewhere.  The  consistent  skewing  of  the  data  indicates  that  sample  variances  will  not  get 
much  lower  than  the  low  values  of  the  sample  sets  observed.  At  the  lowest  variance  levels, 
skewing  is  insignificant.  At  this  level  the  Normal  pdf  becomes  the  most  effective  model. 
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The  mean  only  effects  the  Beta  pdfs’  performance,  skewing  the  pdf  as  described  in 
section  4.2.4. 

The  mean  and  variance  of  the  fighter  and  missile  sample  sets  are  similar  in  two  of 
the  three  aspect  windows  modeled.  Both  nose-on  sample  sets  are  characterized  by  a  mid¬ 
range  mean  and  a  small  variance.  Both  broadside  sample  sets  posess  high-range  means 
and  large  variances.  The  sample  sets  posess  different  characteristics  at  the  tail-on  aspect. 
The  fighter  sample  set  posesses  a  low  range  mean  and  medium  variance,  while  the  missile 
sample  set  posesses  a  mid-range  mean  and  small  variance. 

The  variances  of  the  sample  sets  are  indicators  of  the  number  of  scattering  surfaces 
at  each  aspect  angle,  with  a  small  variance  indicating  few  scatterers  and  a  large  variance 
indicating  many  scatterers.  This  being  the  case,  the  magnitude  of  the  variance  at  each 
aspect  window  of  the  missile  and  fighter  sample  sets  is  typical  of  what  would  be  expected 
for  each  class  of  object. 

The  mean  of  each  sample  set  reflects  the  degree  of  perpendicular  reflection  of  the 
scattering  surfaces  at  each  aspect  window.  Again,  the  experimental  data  matches  what 
would  be  expected  for  both  classes  of  objects. 

The  relatively  consistent  rightward  skewing  of  the  data  was  not  expected.  Every 
sample  set  of  fighter  data  is  rightward  skewed.  However,  each  missile  aspect  angle  is 
skewed  differently:  the  broadside  is  skewed  rightward,  the  tail-on  is  skewed  leftward  and 
the  nose-on  is  not  skewed.  If  the  rightward  skewing  is  a  consistent  trait  in  aircraft  data, 
the  Weibull  would  be  the  most  effective  pdf  for  Window  vs.  Full-Range  models. 

The  MLE  performs  better  for  small  sample  sets,  while  the  MoM  performs  better 
for  large  sample  sets.  There  is  no  apparent  correlation  between  mean  or  variance  and 
parameter  estimator  performance. 
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Figure  4.28.  Top  KGFT  and  CGFT  CDF  of  Fighter  Data  at  Tail-On  ±5°,  0.5°  sample 
interval,  0  =  90°,  6  —  0  polarization,  1  GHz  Full-Scale 


Figure  4.29.  Top  KGFT  and  CGFT  PDF  of  Fighter  Data  at  Tail-On  ±5°,  0.5°  sample 
interval,  9  =  90°,  6-9  polarization,  1  GHz  Full-Scale 
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Figure  4.30.  Top  CGFT  CDF  of  Fighter  Data  at  Tail-On  ±5®,  0.5®  sample  interval, 
9  =  90°,  9-9  polarization,  1  GHz  Full-Scale 


Figure  4.31.  Top  CGFT  PDF  of  Fighter  Data  at  Tail-On  ±5®,  0.5°  sample  interval,  9  = 
90°,  9  —  9  polarization,  1  GHz  Full-Scale 
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Figure  4.32.  Top  KGFT  CDF  of  Fighter  Data  at  Tail-On  ±5°,  0.5®  sample  interval, 
0  =  90°  ±  5°,  9  -  8  polarization,  1  GHz  Full-Scale 


Figure  4.33.  Top  KGFT  PDF  of  Fighter  Data  at  Tail-On  ±5°,  0.5°  sample  interval, 
6  —  90°  ±  5®,  9  -  6  polarization,  1  GHz  Full-Scale 
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Figure  4.34.  Top  CGFT  CDF  of  Fighter  Data  at  Tail-On  ±5°,  0.5°  sample  interval, 
0  =  90°  ±  5°,  9  -  0  polarization,  1  GHz  Full-Scale 


Figure  4.35.  Top  CGFT  PDF  of  Fighter  Data  at  Tail-On  ±5°,  0.5°  sample  interval,  0  = 
90°  ±  5®,  0  -  9  polarization,  1  GHz  Full-Scale 
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Figure  4.36.  Top  KGFT  and  CGFT  CDF  of  Missile  Data  at  Tail-On  ±5°,  0.5°  sample 
interval,  0  =  90°,  0-0  polarization,  18  GHz  Full-Scale 


Figure  4.37.  Top  KGFT  and  CGFT  PDF  of  Missile  Data  at  Tail-On  ±5°,  0.5°  sample 
interval,  0  =  90°,  0-0  polarization,  18  GHz  Full-Scale 
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Figure  4.38.  Top  CGFT  CDF  of  Missile  Data  at  Tail-On  ±5°,  0.5°  sample  interval,  6 
90°,  6  —  6  polarization,  18  GHz  Full-Scale 


Figure  4.39.  Top  CGFT  PDF  of  Missile  Data  at  Tail-On  ±5°,  0.5°  sample  interval,  6 
90°,  6-6  polarization,  18  GHz  Full-Scale 
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V.  Determination  of  Optimum  Model  for  Window  vs.  Window 

Comparisons 


5. 1  Introduction 

This  chapter  will  analyze  how  each  pdf  handles  the  empirical  distribution  of  the  data 
from  each  aspect  window  when  the  statistical  characteristics  are  derived  from  the  aspect 
window  solely.  The  aspect  windows  and  types  of  aircraft  studied  will  be  identical  to  that 
of  Chapter  IV. 

Because  the  statistical  characteristics  are  derived  solely  from  the  aspect  windows, 
the  normalized  range  of  each  sample  set  will  stretch  from  zero  to  one.  This  is  in  contrast 
to  the  Window  vs.  Full-Range  comparisons  where  the  sample  set  could  be  restricted  to  a 
small  section  of  the  full  range  and  usually  was. 

The  expanded  range  of  the  sample  sets  should  result  in  a  consistently  larger  sample 
variance. 

The  Kolmogorov  Level  of  Significance,  and  the  Chi-Square  Level  of  Significance  for 
Window  vs.  Full-Range  RCS  Models  are  presented  in  Appendix  A. 

5.2  Test  Results  for  Window  vs.  Window  RCS  Models 

5.2.1  Nose-On  Fighter.O  =  90°  This  data  set  is  characterized  by  a  mean  of  0.74,  a 
variance  of  9.07  x  10-2,  a  small  number  of  samples  ( n  =  11),  and  a  strong  rightward  skew. 
The  best  fit  KGFT  and  CGFT  models  are  presented  in  Figures  5. 1-5.4. 

This  variance  is  four  times  larger  than  the  largest  Window  vs.  Full- Range  variance. 
The  mean  is  also  larger  than  most  of  the  Window  vs.  Full-Range  means.  The  wide 
variance  expands  the  pdfs  so  they  can  cover  a  sample  set  than  is  beginning  to  appear 
uniform.  The  KGFT  declares  the  Beta  MLE  and  MoM  pdfs  to  have  a  LoS  greater  than 
20%.  The  Beta  MLE  has  the  lowest  Kolmogorov  test  statistic.  However,  for  this  sample 
set,  visual  assessment  agrees  with  the  CGFT,  which  declares  the  Beta  MoM  to  be  the  only 
pdf  with  a  LoS  greater  than  25%.  The  average  separation  between  the  Beta  MoM  and 
the  sample  set  (0.0544)  is  much  smaller  than  the  difference  between  the  Beta  MLE  and 
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the  sample  set  (0.1275).  However,  the  ra  =  10  sample  climbs  at  a  much  higher  rate  than 

the  Beta  MoM.  For  this  data  sample,  the  Beta  MoM  is  less  expanded  than  the  MLE.  The 

* 

parameters  for  the  Beta  function  are  less  than  one  for  the  MoM  and  the  MLE.  This  causes 
the  pdf  to  approximate  a  negative  exponential  function  against  the  low-range  boundary 
and  a  positive  exponential  function  against  the  high-range  boundary.  The  MoM  pdf  has 
unequal  parameters  which  skew  the  emphasis  of  the  two  exponential  functions  toward  the 
positive  exponential  function  against  the  high-range  boundary.  This  closely  matches  the 
data  histogram.  The  MLE  has  small,  equal  parameters,  which  press  the  two  exponential 
functions  against  their  respective  boundaries.  This  is  the  result  of  an  extreme  form  of 
expansion. 

No  other  pdf  could  approximate  the  sample  set  because  of  the  high  variance. 

5.2.2  Nose-On  Fighter:  6  =  90°  ±  5°  This  data  set  is  characterized  by  a  mean  of 
0.64,  a  variance  of  2.9  x  10“2,  a  large  number  of  samples  (n  =  121),  and  no  discernable 
skew.  The  best  fit  KGFT  and  CGFT  models  are  presented  in  Figures  5. 5-5.8. 

Though  a  non-skewed  sample  set  with  a  mean  above  0.5  can  mean  trouble  for  the 
Beta  pdfs,  the  wide  variance  expands  all  other  pdfs  to  a  greater  extent.  The  KGFT  declares 
the  Beta  MoM  to  have  the  highest  LoS,  though  it  is  only  between  10%  and  20%.  There  are 
no  other  pdfs  with  that  range  of  LoS.  The  lack  of  skew  with  a  mean  shifted  off  mid-range 
hurts  the  Beta  pdfs,  but  the  large  variance  prevents  the  Normal,  Lognormal,  or  Weibull 
pdfs  from  compressing  enough  to  match  or  improve  the  performance  of  the  Beta  pdfs.  A 
visual  assessment  shows  the  Weibull  MoM  and  Normal  pdfs  matching  the  Beta  MoM  in 
tracking  the  ECDF.  The  CGFT  declares  the  Normal  and  Weibull  pdfs  to  have  the  lowest 
test  statistic.  The  CGFT  gives  the  pdfs  a  LoS  between  5%  and  2.5%. 

5.2.3  Nose-On  Missile:  0  =  90°  This  data  set  is  characterized  by  a  mean  of  0.50, 
the  widest  variance  of  any  sample  set  (var=0.1083),  a  small  number  of  samples  (n  =  21), 
and  no  discernable  skew.  The  best  fit  KGFT  and  CGFT  models  are  presented  in  Figures 
5.9-5.10. 
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The  sample  set  appears  to  have  a  fairly  uniform  distribution.  The  Beta  pdfs  will 
present  a  uniform  distribution  if  a  =  =  1.  If  a  and  /?  are  equal  but  less  than  one,  the 

pdf  is  expanded  toward  the  boundaries  of  the  range.  The  Beta  MoM  is  the  least  expanded 
of  the  Beta  pdfs  and  attains  the  lowest  Kolmogorov  test  statistic.  The  Beta  MLE  and 
MoM,  Normal,  and  Rayleigh  pdfs  attain  a  Kolmogorov  LoS  greater  than  20%.  Because 
the  mean  of  the  sample  set  is  mid-range,  the  Beta  pdfs  do  not  skew  and  therefore  are  able 
to  maintain  a  fairly  constant  probability  density  throughout  the  range. 

The  CGFT  agrees  with  the  KGFT,  giving  the  Beta  MoM  a  LoS  greater  than  20%. 
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Figure  5.1.  Top  KGFT  CDF  of  Fighter  Data  at  Nose-On  ±5°,  0.5°  sample  interval,  9 
90°,  6  —  0  polarization,  1  GHz  Full-Scale 
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Nose-On  RCS  vs.  Window-Window  Beta  MLE  PDF 


Figure  5.2.  Top  KGFT  PDF  of  Fighter  Data  at  Nose-On  ±5°,  0.5°  sample  interval,  6 
90°,  9  —  9  polarization,  1  GHz  Full-Scale 
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Nose-On  RCS  vs.  Window-Window  Beta  MoM  CDF 


Figure  5.3. 
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Figure  5.4. 


Top  CGFT  CDF  of  Fighter  Data  at  Nose-On  ±5®,  0.5°  sample  interval,  6 
90°,  6  —  6  polarization,  1  GHz  Full-Scale 


Top  CGFT  PDF  of  Fighter  Data  at  Nose-On  ±5°,  0.5°  sample  interval,  9 
90°,  9  —  6  polarization,  1  GHz  Full-Scale 
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Nose-On  RCS  vs.  Window-Window  Beta  MoM  CDF 

Figure  5.5.  Top  KGFT  CDF  of  Fighter  Data  at  Nose-On  ±5®,  0.5°  sample  interval,  0 
90°  ±  5°,  0  —  0  polarization,  1  GHz  Full-Scale 


Figure  5.6.  Top  KGFT  PDF  of  Fighter  Data  at  Nose-On  ±5°,  0.5°  sample  interval,  8 
90°  ±5°,  0  —  9  polarization,  1  GHz  Full-Scale 


5-6 


Figure  5.7.  Top  CGFT  CDF  of  Fighter  Data  at  Nose-On  ±5°,  0.5°  sample  interval,  6 
90°  ±  5°,  0  -  6  polarization,  1  GHz  Full-Scale 


5 

4.5 
4 

3.5 
3 

p(x)  2.5 
2 

1.5 
1 

0.5 
0 

0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

Nose-On  RCS/5  vs.  Window-Window  Weibull  MoM  PDf 

Figure  5.8.  Top  CGFT  PDF  of  Fighter  Data  at  Nose-On  ±5°,  0.5°  sample  interval,  6 
90°  ±  5°,  9  -  9  polarization,  1  GHz  Full-Scale 
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Figure  5.9.  Top  KGFT  and  CDFT  CDF  of  Missile  Data  at  Nose-On  ±5°,  0.5°  sample 
interval,  9  -  90°,  6-6  polarization.  18  GHz  Full-Scale 


Nose-On  RCS  vs.  Window-Window  Beta  MoM  PDF 


Figure  5.10.  Top  KGFT  and  CGFT  PDF  of  Missile  Data  at  Nose-On  ±5°,  0.5°  sample 
interval,  9  =  90°,  6  —  6  polarization,  18  GHz  Full-Scale 
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5.2.4  Broadside  Fighter:  6  =  90°  This  data  set  is  characterized  by  a  mean  of  0.66, 
a  variance  of  3.9  x  10-2,  a  small  number  of  samples  (n  =  21),  and  no  discernable  skew. 
The  best  fit  KGFT  and  CGFT  models  are  presented  in  Figures  5.11-5.14. 

The  KGFT  declares  the  Beta  MoM,  Weibull  MoM,  and  Normal  pdfs  to  have  a  LoS 
greater  than  20%.  The  Weibull  MLE’s  a  parameter  is  expanded  below  a  value  of  three, 
so  the  pdf  shifts  far  to  the  left,  where  it  begins  to  approximate  the  negative  exponential 
function.  Both  Beta  MLE  parameters  are  expanded  below  one,  so  the  pdf  is  pressed 
against  both  range  boundaries.  However,  the  small  variance  (for  Window  vs.  Window 
comparisons)  requires  a  pdf  that  can  approximate  a  normal  function.  The  Weibull  MoM 
and  Normal  pdfs  have  normal  shaped  functions.  Because  the  mean  is  above  0.5,  the  Beta 
MoM  is  skewed  right,  but  the  Beta  MoM  is  able  to  maintain  a  higher  compression  than 
the  Weibull  MoM  and  Normal  pdfs,  and  attains  the  lowest  Kolmogorov  test  statistic. 

Because  their  densities  are  centered  correctly,  the  CGFT  declares  the  Weibull  MoM 
and  Normal  pdfs  to  have  the  highest  LoS  (25%  >  a  >  10%).  Visual  assessment  tends  to 
agree  with  the  CGFT  again. 

5.2.5  Broadside  Fighter:  6  =  90°  ±  5°  This  data  set  is  characterized  by  a  mean  of 
0.64,  a  variance  of  1.8  X  10-2,  a  large  number  of  samples  (n  —  231),  and  a  rightward  skew. 
The  best  fit  KGFT  and  CGFT  models  are  presented  in  Figures  5.15-5.16. 

The  sample  variance  is  the  smallest  of  the  Window  vs.  Window  comparisons.  This 
allows  the  Weibull  MoM  to  compress  and  become  rightward  skewed  to  the  extent  that 
the  Weibull  MoM  attains  the  lowest  Kolmogorov  test  statistic,  although  none  of  the  pdfs 
match  the  compression  of  the  sample  set.  The  Beta  MoM  and  Weibull  MoM  pdfs  attain 
the  highest  LoS  (20%  >  a  >  10%). 

The  CGFT  declares  the  Weibull  MoM  to  have  the  lowest  test  statistic,  though,  as 
with  all  other  large  sample  sets,  the  test  statistic  results  in  a  poor  LoS. 

5.2.6  Broadside  Missile:  6  =  90°  This  data  set  is  characterized  by  a  mean  of  0.71, 
a  variance  of  4.46  x  10-2,  a  small  number  of  samples  (n  —  41),  and  a  rightward  skew.  The 
best  fit  KGFT  and  CGFT  models  are  presented  in  Figures  5.17-5.20. 
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The  Beta  MLE  is  expanded  into  a  postive  exponential  function  by  the  very  high 
mean.  The  Beta  MoM  does  not  skew  to  the  extent  of  the  Beta  MLE  and  approximates  a 
ramp  function.  The  Weibull  pdfs  approach  a  normal  function  and,  along  with  the  Normal 
pdf,  reach  their  maximum  probability  density  too  early.  The  KG  FT  declares  the  Beta 
MoM,  Weibull  MoM,  and  Normal  pdfs  to  have  a  LoS  greater  than  20%.  The  Beta  MoM 
has  the  smallest  test  statistic,  with  a  value  less  than  half  that  of  the  Weibull  MoM  or 
Normal  pdfs. 

The  CGFT  agrees  with  the  KGFT,  giving  the  Beta  MoM  a  LoS  greater  than  25%. 
Visual  assessment  clearly  shows  the  Beta  MoM  to  be  the  best-fit  pdf. 
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Broadside  R.CS  vs.  Window-Window  Beta  MoM  CDF 


Figure  5.11.  Top  KGFT  CDF  of  Fighter  Data  at  Broadside  ±5°,  0.5°  sample  interval, 
0  =  90°,  0-0  polarization,  1  GHz  Full-Scale 


Figure  5.12.  Top  KGFT  PDF  of  Fighter  Data  at  Broadside  ±5°,  0.5®  sample  interval, 
0  =  90°,  0-0  polarization,  1  GHz  Full-Scale 
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Broadside  RCS  vs.  Window- Window  Weibull  MoM  CDF 


Figure  5.13.  Top  CGFT  CDF  of  Fighter  Data  at  Broadside  ±5°,  0.5°  sample  interval, 
6  =  90°,  0  —  0  polarization,  1  GHz  Full-Scale 


Broadside  RCS  vs.  Window- Window  Weibull  MoM  PDF 


Figure  5.14.  Top  CGFT  PDF  of  Fighter  Data  at  Broadside  ±5°,  0.5°  sample  interval, 
0  =  90°,  0-0  polarization,  1  GHz  Full-Scale 
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Broadside  RCS  vs.  Window-Window  Weibull  MoM  CDF 

Figure  5.15.  Top  KGFT  and  CGFT  CDF  of  Fighter  Data  at  Broadside  ±5°,  0.5*  sample 
interval,  9  =  90°  ±5°,  9  —  0  polarization,  1  GHz  Full-Scale 


Broadside  RCS/5  vs.  Window-Window  Weibull  MoM  PDF 

Figure  5.16.  Top  KGFT  and  CGFT  PDF  of  Fighter  Data  at  Broadside  ±5°,  0.5°  sample 
interval,  9  =  90*  ±  5°,  9  —  9  polarization,  1  GHz  Full-Scale 
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Figure  5.17.  Top  KGFT  and  CGFT  CDF  of  Missile  Data  at  Broadside  ±5°,  0.5°  sample 
interval,  9  =  90°,  9  —  9  polarization,  18  GHz  Full-Scale 


Broadside  RCS  vs.  Window-Window  Beta  MoM  PDF 


Figure  5.18.  Top  KGFT  and  CGFT  PDF  of  Missile  Data  at  Broadside  ±5°,  0.5°  sample 
interval,  9  =  90°,  9  —  9  polarization,  18  GHz  Full-Scale 


5-14 


5.2. 7  Tail-On  Fighter:  9  —  90°  This  data  set  is  characterized  by  a  mean  of  0.66,  a 
variance  of  6.92  x  10-2,  a  small  number  of  samples  (n  =  11),  and  no  discernible  skew.  The 
best-fit  KGFT  and  CGFT  models  are  presented  in  Figures  5.21-5.24. 

Though  the  variance  is  not  the  largest  of  the  sample  sets,  the  samples  appear  almost 
uniformly  distributed.  There  are  many  pdfs  to  which  the  KGFT  gives  an  LoS  greater  than 
20%:  Beta  MLE  and  MoM,  Weibull  MoM,  Normal,  and  One-Dominant-Plus-Rayleigh. 
The  Beta  MLE  is  expanded  to  the  extent  that  the  exponential  functions  are  pressed  tight 
against  the  range  boundaries.  The  Beta  MoM  is  a  positive  exponential  function  with  a 
very  shallow  slope  until  very  close  to  the  high-range  boundary.  The  pdf  approximates  a 
uniform  distribution  to  a  certain  extent.  The  Weibull  MoM,  Normal,  and  One-Dominant- 
Plus-Rayleigh  pdfs  are  all  normal  functions  centered  near  the  sample  mean.  The  KGFT 
gives  the  Beta  MoM  the  lowest  test  statistic. 

The  CGFT  declares  the  Beta  MLE  and  MoM  pdfs  to  have  the  highest  LoS  (a  >  25%). 
Though  the  Beta  MLE  may  be  evenly  distributed  about  the  mean,  a  visual  assessment 
indicates  that  the  Beta  MLE  does  not  come  close  to  approximating  the  sample  set. 

5.2.8  Tail-On  Fighter:  9  =  90°  ±  5°  This  data  set  is  characterized  by  a  mean  of 
0.68,  a  variance  of  4.25  x  10-2,  a  large  number  of  samples  (n  =  121),  and  no  discernible 
skew.  The  best-fit  KGFT  and  CGFT  models  are  presented  in  Figures  5.25-5.28. 

The  variance  is  tighter  for  this  sample  set  than  the  Tail-On  Fighter:  9  =  90°  sample 
set.  The  KGFT  declares  the  Beta  MoM  pdf  to  have  a  LoS  greater  than  20%,  along  with 
the  lowest  test  statistic.  As  has  been  seen  in  many  other  sample  sets,  the  Beta  MoM  is  the 
only  pdf  able  to  maintain  high  compression  at  wide  variances  and  a  rightward  skew.  The 
Weibull  MoM  and  Normal  pdfs  produce  normal  functions  centered  on  the  sample  mean. 
As  variance  increases,  the  Weibull  pdf  becomes  normal  shaped,  then  begins  to  skew  left. 

The  CGFT  declares  the  Weibull  MoM  and  Beta  MoM  pdfs  to  have  the  highest  LoS, 
though,  as  with  all  the  large  sample  sets,  the  LoS  is  poor  (5%  >  a  >  2.5%).  The  Weibull 
MoM  has  the  lowest  test  statistic.  A  visual  assessment  shows  the  Beta  MoM  to  have  the 
best  fit. 
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5.2.9  Tail-On  Missile:  9  —  90°  This  data  set  is  characterized  by  a  mean  of  0.43, 
a  variance  of  0.103,  a  small  number  of  samples  (n  =  21),  and  slight  leftward  skew.  The 
best-fit  KGFT  and  CGFT  models  are  presented  in  Figures  5.29-5.30. 

With  the  second  highest  variance,  this  sample  set  is  nearly  uniform.  All  parameters 
for  the  Beta  pdfs  are  less  than  one,  so  they  are  expanded  to  opposing  exponential  functions. 
As  usual,  the  Beta  MoM  is  less  expanded  than  the  Beta  MLE,  so  the  Beta  MoM  has  a 
fairly  high  probability  density  across  the  sample  range.  The  Weibull  MoM  skews  leftward, 
but  skews  too  far.  The  high-range  probability  densities  are  not  large  enough.  The  Normal 
pdf  puts  too  much  emphasis  on  the  middle  range  of  the  sample  set.  The  KGFT  declares 
the  Beta  MoM  and  Weibull  MoM  to  have  a  LoS  greater  than  20%.  The  Beta  MoM  has 
the  lowest  test  statistic. 

The  CGFT  agrees  with  the  KGFT,  declaring  the  Beta  MoM  to  be  the  only  pdf  with 
a  LoS  greater  than  25%.  Visual  assessment  agrees  with  the  two  tests. 

5.3  Analysis  of  Window  vs.  Window  RCS  Models 

5.3.1  PDF  Performance  The  Window  vs.  Window  sample  sets  have  a  consistently 
high  mean:  Avg.  mean  =  0.68,  standard  deviation  =  9.4  x  10~2.  The  variance  ranges 
from  1.8  X  10~2  to  1.02  X  10_1.  This  range’s  low  end  is  the  high  end  of  the  Window  vs. 
Full-Range  variance  range,  and  the  range’s  high  end  is  a  level  of  magnitude  above  the  top 
of  the  Window  vs.  Full-Range  variance  range.  Consequently,  the  variance  sensitive  pdfs 
(Weibull,  Lognormal,  Normal)  are  expanded  much  more  than  the  Beta  pdfs.  The  Weibull 
MLE  shape  parameter  has  a  value  less  than  one  four  times,  at  these  values,  the  Weibull 
pdf  approximates  the  negative  exponential  function.  Because  only  one  of  the  sample  sets 
is  skewed  left,  the  extreme  expansion  of  the  Weibull  MLE  makes  the  pdf  ineffective  for 
Window  vs.  Window  modeling.  The  Weibull  MoM  is  not  as  variance  sensitive,  but  is  still 
skewed  left  when  the  sample  sets  are  not. 

The  scale  parameter  of  the  Lognormal  pdf  causes  the  pdf  to  approximate  a  negative 
exponential  function  when  the  parameter  is  greater  than  one.  This  is  the  case  in  eight  of 
the  nine  sample  sets.  In  the  ninth  case,  the  sample  set  has  the  smallest  variance.  The 
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pdf  is  still  centered  far  to  the  left  of  the  sample  set  values.  The  lognormal  pdf  is  more 
sensitive  to  variance  changes  than  any  other  pdf.  The  wide  variance  skews  the  pdf  to  the 
left  consistently,  making  the  pdf  ineffective. 

The  Normal  pdf  does  not  expand  as  fast  as  the  Lognormal.  When  the  Beta  pdf  is 
skewed  too  far  to  the  right  for  the  sample  mean,  the  Normal  is  able  to  outperform  the 
Beta  in  the  CGFT.  This  only  happens  at  the  lower  sample  variances. 

The  Weibull  MoM  performs  well  at  the  same  low  variances  that  the  Normal  performs 
well  on.  At  these  variances,  the  Weibull  MoM  outperforms  the  Beta  MoM  on  the  CGFT. 
At  the  lowest  Window  vs.  Window  sample  set,  the  Weibull  MoM  outperforms  the  Beta 
MoM  on  the  KGFT.  For  the  three  lowest  sample  variances,  a  visual  assessment  indicates 
the  Weibull  MoM  to  be  as  good  or  better  a  model  than  the  Beta  MoM. 

The  Rayleigh  and  One-Dominant-Plus-Rayleigh  pdfs  perform  slightly  better  than 
they  do  for  Window  vs.  Full-Range  models.  The  high  variances  expand  the  sample  sets  to 
the  extent  that  the  permanently  expanded  Rayleigh  class  pdfs  can  occasionally  model  the 
sample  sets  well.  This  permanent  expansion  is  due  to  consistent  sample  second  moments 
which  are  not  significantly  different  in  value  from  the  Window  vs.  Full-Range  sample  sets. 

The  Rayleigh  ( y  =  x 2)  pdf  is  consistently  ineffective,  as  it  is  for  the  Window  vs. 
Full-Range  models. 

The  Beta  pdf  is  made  for  the  Window  vs.  Window  sample  sets.  The  pdf  remains 
expansion  resistant  to  high  variances.  The  high  variances  result  in  a  multitude  of  sample 
distributions  which  also  favor  the  Beta  over  other  pdfs.  The  expansion  resistance  and 
versatility  of  the  pdf  give  it  as  high  a  Kolmogorov  LoS  as  any  pdf  in  every  sample  set.  The 
Beta  pdf  has  the  lowest  Kolmogorov  test  statistic  in  eight  of  the  nine  sample  sets.  The 
ninth  sample  set  posesses  the  smallest  variance  of  the  group. 

Perhaps  because  the  pdf  skews  at  the  routinely  high  means  of  the  Window  vs.  Win¬ 
dow  sample  sets,  the  Beta  pdf  does  not  perform  as  well  on  the  Chi-Square  tests.  This 
would  be  caused  by  the  uneven  distribution  about  the  sample  set  mean  that  the  Beta  pdf 
exhibits.  The  Weibull  and  Normal  pdfs  outperform  the  Beta  on  the  sample  sets  with  the 
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four  lowest  variances  when  the  CGFT  is  used.  Visual  assessment  agrees  with  the  CGFT 
once  out  of  the  three  cases  where  there  is  a  difference  between  the  KGFT  and  the  CGFT. 

5.3.2  Parameter  Estimation  The  MoM  outperforms  the  MLE  in  every  sample  set. 
In  the  one  sample  set  where  the  Beta  MLE  is  rated  better  than  the  MoM  by  the  KGFT, 
a  visual  assessment  clearly  indicates  the  Beta  MoM  to  have  the  better  fit.  Though  the 
consistently  wide  variances  may  be  the  cause  of  the  MoM’s  success  in  Window  vs.  Window 
models,  variance  plays  no  apparent  role  in  determining  the  optimum  parameter  estimator 
in  Window  vs.  Full-Range  models.  The  determiner  in  those  models  appears  to  be  the 
number  of  samples  in  the  sample  set.  In  the  Window  vs.  Window  models,  sample  size 
makes  no  contribution  toward  parameter  estimator  effectiveness. 

5.3.3  Sample  Set  Behaviour  The  wide  sample  variance  experienced  matches  expec¬ 
tations  for  this  type  of  model.  The  high  mean  of  the  Window  vs.  Window  models  matches 
the  rightward  skew  of  Window  vs.  Full-Range  models.  The  high  mean  of  the  Window  vs. 
Window  sample  sets  could  be  a  result  of  the  relatively  high  RCS  magnitudes  encountered 
at  the  modeled  windows.  The  relatively  high  magnitudes  would  cause  the  majority  of 
the  samples  to  be  large  along  with  a  few  small  samples.  This  hypothesis  is  backed  up  by 
the  fact  that  the  two  sample  sets  that  do  not  have  large  means  occur  in  the  Tail-On  and 
Nose-On  windows,  which  do  not  have  as  large  RCS  magnitudes  as  the  Broadside  windows. 

5-4  Summary 

This  chapter  studies  the  results  of  the  Window  vs.  Window  RCS  modeling  procedure. 
It  shows  that  the  variance  of  the  sample  sets  is  much  larger  than  the  variance  of  the  Window 
vs.  Full-Range  sample  sets.  The  mean  is  found  to  be  consistently  high  in  the  sample  range. 

It  is  shown  that,  for  the  large  variances  exhibited,  the  Beta  pdf  is  the  only  pdf  able 
to  handle  the  resulting  large  variety  of  sample  distributions. 

The  best  parameter  estimator  is  clearly  the  Method  of  Moments.  The  reason  for  this 
could  be  the  large  sample  variances,  but  this  is  uncertain. 
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Figure  5.19.  Top  KGFT  and  CGFT  CDF  of  Fighter  Data  at  Tail-On  ±5®,  0.5®  sample 
interval,  0  =  90°,  0-0  polarization,  1  GHz  Full-Scale 


Figure  5.20.  Top  KGFT  and  CGFT  PDF  of  Fighter  Data  at  Tail-On  ±5°,  0.5°  sample 
interval,  6  —  90°,  0-0  polarization,  1  GHz  Full-Scale 
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Figure  5.21.  Top  CGFT  CDF  of  Fighter  Data  at  Tail-On  ±5®,  0.5®  sample  interval, 
6  =  90°,  9-6  polarization,  1  GHz  Full-Scale 


Figure  5.22.  Top  CGFT  PDF  of  Fighter  Data  at  Tail-On  ±5®,  0.5®  sample  interval,  6  = 
90®,  9  -  6  polarization,  1  GHz  Full-Scale 
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Tail-On  RCS/5  vs.  Window- Window  Beta  MoM  CDF 

Figure  5.23.  Top  KGFT  CDF  of  Fighter  Data  at  Tail-On  ±5°,  0.5°  sample  interval, 
6  =  90°  ±  5°,  6  —  0  polarization,  1  GHz  Full-Scale 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

Tail-On  RCS/5  vs.  Window-Window  Beta  MoM  PDF 

Figure  5.24.  Top  KGFT  PDF  of  Fighter  Data  at  Tail-On  ±5°,  0.5°  sample  interval, 
6  =  90°  ±  5°,  0  -  9  polarization,  1  GHz  Full-Scale 
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Figure  5.25.  Top  CGFT  CDF  of  Fighter  Data  at  Tail-On  ±5°,  0.5°  sample  interval, 
9  =  90°  ±  5°,  0  -  6  polarization,  1  GHz  Full-Scale 


Figure  5.26.  Top  CGFT  PDF  of  Fighter  Data  at  Tail-On  ±5°,  0.5°  sample  interval,  9  = 
90°  ±  5°,  9  -  9  polarization,  1  GHz  Full-Scale 
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Figure  5.27. 


P(x) 


Figure  5.28. 


Top  KGFT  and  CGFT  CDF  of  Missile  Data  at  Tail-On  ±5°,  0.5°  sample 
interval,  6  =  90°,  0-0  polarization,  18  GHz  Full-Scale 


Top  KGFT  and  CGFT  PDF  of  Missile  Data  at  Tail-On  ±5°,  0.5°  sample 
interval,  6  =  90°,  6-0  polarization,  18  GHz  Full-Scale 
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VI.  Conclusions  and  Recommendations 


6.0.1  Conclusions  1.  The  statistical  characteristics  of  the  three  modeled  aspect 
windows  can  be  predicted.  These  characteristics  will  determine  the  optimum  pdf. 

2.  The  Weibull  pdf  will  consistently  be  the  optimum  pdf  for  narrow  variance  sample 
sets.  The  Beta  pdf  will  consistently  be  the  optimum  pdf  for  wide  variance  sample  sets. 

i 

3.  Though  never  quite  as  reliable  as  Window  vs.  Full-Range  models,  Window  vs. 
Window  models  perform  satisfactorily  for  the  Kolmogorov  Goodness-of-Fit  test  under  the 
specifications  of  Section  2.5.  The  decrease  in  reliability  is  probably  due  to  the  wider 
variance  of  the  Window  vs.  Window  sample  sets. 

4.  The  Method  of  Moments  will  be  the  superior  parameter  estimator  for  Window  vs. 
Window  models.  The  Maximum  Likelihood  Estimator  tends  to  be  the  superior  parameter 
estimator  for  small-sample-set  Window  vs.  Full-Range  models,  the  Method  of  Moments 
tends  to  be  the  superior  parameter  estimator  for  large-sample-set  Window  vs.  Full-Range 
models. 

5.  When  determining  goodness-of-fit,  a  visual  assessment  should  always  be  performed 
to  insure  the  pdf  with  the  lowest  test  statistic  is  the  closest  match  to  the  ECDF.  Both  the 
Kolmogorov  and  Chi-Square  tests  can  give  erroneous  results  under  certain  conditions.  One 
bad  sample  can  cause  a  close  fitting  pdf  to  fail  the  Kolmogorov  test,  while  skewed  sample 
sets  can  cause  problems  for  the  Chi-Square  test. 

6.0.2  Recommendations  1.  Build  signal  detection  and  signal  identification  models 
using  the  optimized  static  RCS  models  as  the  base  RCS  value.  Use  numerical  methods  to 
determine  the  received  detected  waveform  plus  noise. 

2.  Form  methods  for  measuring  data  skew.  This  would  provide  a  method  for  better 
explaining  why  a  particular  pdf  is  optimum.  Use  numerical  methods  to  determine  the  skew 
of  the  pdf  models. 

3.  Evaluate  more  data  sets  to  verify  the  statistical  characteristics  of  the  sample  sets 
noted  in  this  study. 
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4.  Determine  the  causes  of  the  varying  performances  of  the  MoM  and  MLE  parameter 
estimators. 
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Appendix  A.  Test  Statistics  and  Levels  of  Significance 

The  Kolmogorov  and  Chi-Square  results  are  presented  for  every  aspect  angle  and 
both  aircraft  types. 
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Table  A.l.  Fighter:  8  =  90°  Kolmogorov  Levels  of  Significance  for  Window  vs.  Full 
Range,  n=number  of  samples,  7\=Test  Statistic,  Wp=l-ot  Quantile 
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a  >20% 
a  >20% 
a  >20% 
a  >20% 
1%>  a 
1%>  a 
1%>  a 
2%> a >1% 


a  >20% 
a  >20% 
a  >20% 
a  >20% 
20%  >  a  >10% 
a  >20% 
1%>  a 
1%>  a 
1  %>  a 
1%>  a 


a  >20% 
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a  >20% 
a  >20% 
a  >20% 
a  >20% 
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1%>  a 
10%>  a  >5% 
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Table  A. 2.  Fighter:  0  =  90°  Chi-Square  Levels  of  Significance  for  Window  vs.  Full- 
Range,  n=number  of  samples,  DoF=  Degrees  of  Freedom,  Ti=Test  Statistic, 
Wp=l-a  Quantile 
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Table  A.3.  Missile:  0  =  90°  Kolmogorov  Levels  of  Significance  for  Window  vs.  Full- 
Range,  n=number  of  samples,  T2=  Test  Statistic,  Wp=l-a  Quantile 
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a  >20% 

Rayleigh 

0.5076 

0.255 

1%>  a 

Rayleigh-y  MLE 

0.7092 

0.255 

1%>  a 

Rayleigh-y  MoM 

HI 

0.255 

1%>  a 

Rayleigh-)- 1 

EB 

0.255 

1%>  a 

Beta  MLE 
Beta  MoM 

180° 

21 

0.4392 

0.344 

1%>  a 

Weibull  MLE 

0.1248 

0.226 

a  >20% 

Weibull  MoM 

0.1355 

0.226 

a  >20% 

Lognormal 

0.1018 

0.226 

a  >20% 

Normal 

0.1033 

0.226 

a  >20% 

Rayleigh 

0.6373 

0.344 

1%>  a 

Rayleigh-y  MLE 

0.8382 

0.344 

1%>  a 

Rayleigh-y  MoM 

0.8051 

0.344 

1%>  a 

Rayleigh-)- 1 

0.5624 

0.344 

1%>  a 

Table  A.4.  Missile:  6  =  90°  Chi-Square  Levels  of  Significance  for  Window  vs.  Full-Range, 
n=number  of  samples,  DoF  =  Degrees  of  Freedom,  Tx =Test  Statistic,  WP= 1-a 
Quantile 


pdf 

Window 

Center 

n 

DoF 

Tx 

wp 

Level  of  Significance 

Beta  MLE 
Beta  MoM 

0° 

21 

1 

28.7143 

10.83 

0.1%>  a 

Wei  bull  MLE 

1 

4.3333 

5.024 

5%>  a  >2.5% 

Weibull  MoM 

1 

4.3333 

5.024 

5%>  a  >2.5% 

Lognormal 

1 

4.3333 

5.024 

5%>  a  >2.5% 

Normal 

1 

4.3333 

5.024 

5%>  a  >2.5% 

Rayleigh 

2 

63 

13.82 

0.1%>  a 

Rayleigh-y  MLE 

2 

63 

13.82 

0.1%>  a 

Rayleigh-y  MoM 

2 

25.6667 

13.82 

0.1%>  a 

Rayleigh+1 

2 

63 

13.82 

0.1%>  a 

Beta  MLE 

90° 

41 

2 

0.3415 

2.773 

a  >25% 

Beta  MoM 

2 

1.3171 

2.773 

a  >25% 

Weibull  MLE 

2 

3.2683 

4.605 

a  >25% 

Weibull  MoM 

2 

3.26830 

2.773 

25%>  a  >10% 

Lognormal 

2 

17.5366 

13.82 

0.1%>  a 

Normal 

2 

5.4634 

5.991 

fig 

iO 

A 

8 

A 

o 

rH 

Rayleigh 

3 

38.6829 

16.27 

0.1%>  a 

Rayleigh-y 

3 

127.1207 

16.27 

0.1%>  a 

Rayleigh-y 

3 

127.1207 

16.27 

0.1%>  a 

Rayleigh+1 

3 

21.8537 

16.27 

0.1%>  a 

Beta  MLE 
Beta  MoM 

00 

o 

o 

21 

1 

25.6667 

10.83 

0.1%>  a 

Weibull  MLE 

1 

2.4286 

2.706 

25%>  a  >10% 

Weibull  MoM 

1 

2.4286 

2.706 

25%  >  a  >10% 

Lognormal 

1 

2.4286 

2.706 

25%  >  a  >10% 

Normal 

1 

2.5286 

2.706 

25%>  a  >10% 

Rayleigh 

2 

63 

13.82 

0.1%>  a 

Rayleigh-y 

2 

63 

13.82 

0.1%>  a 

Rayleigh-y 

2 

63 

13.82 

0.1%>  a 

Rayleigh+1 

2 

63 

13.82 

0.1%>  a 
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Table  A.5.  Fighter:  6  =  90*  ±  5°  Kolmogorov  Levels  of  Significance  for  Window  vs 
Full-Range,  n=number  of  samples,  Ti=Test  Statistic,  Wp  =  l-a  Quantile 


pdf 

Window 

Center 

n 

H 

wp 

Level  of  Significance 

Beta  MLE 

0° 

121 

0.3745 

0.148 

1%>  a 

Beta  MoM 

Weibull  MLE 

0.0733 

. 

a  >20% 

Wei  bull  MoM 

0.0553 

: 

a  >20% 

Lognormal 

0.1134 

0.124 

10%>  a  >5% 

Normal 

0.1059 

0.111 

20%  >  a  >10% 

Rayleigh 

0.5493 

0.148 

1%>  a 

Rayleigh-y  MLE 

0.7867 

0.148 

1%>  a 

Rayleigh-y  MoM 

0.6893 

0.148 

1%>  a 

Rayleigh -(-1 

0.4817 

0.148 

1%>  a 

Beta  MLE 

90° 

231 

0.1274 

0.107 

1%>  a 

Beta  MoM 

0.0877 

0.0895 

10% >  a  >5% 

Weibull  MLE 

0.0655 

0.070 

a  >20% 

Weibull  MoM 

0.0432 

0.070 

a  >20% 

Lognormal 

0.1248 

0.107 

1%>  a 

Normal 

0.0906 

0.100 

5%> a  >2% 

Rayleigh 

0.4475 

0.107 

1%>  a 

Rayleigh-y  MLE 

0.7008 

0.107 

1%>  a 

Rayleigh-y  MoM 

0.6517 

0.107 

1%>  a 

Ravleigh+1 

0.3334 

0.107 

1%>  a 

Beta  MLE 

O 

O 

00 

iH 

121 

0.1276 

5%> a  >2% 

Beta  MoM 

0.0986 

Q II 

20%  >  a  >10% 

Weibull  MLE 

0.0490 

| 

a  >20% 

Weibull  MoM 

0.0655 

0.097 

a  >20% 

Lognormal 

0.1364 

0.138 

5%>  a  >2% 

Normal 

0.0816 

0.097 

a  >20% 

Rayleigh 

0.3261 

0.148 

1%>  a 

Rayleigh-y  MLE 

0.5944 

0.148 

1%>  a 

Rayleigh-y  MoM 

0.3356 

0.148 

1%>  a 

Rayleigh+1 

0.2291 

0.148 

1%>  a 
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Table  A.6.  Fighter:  6  =  90°  ±  5°  Chi-Square  Levels  of  Significance  for  Window  vs.  Full- 
Range,  n=number  of  samples,  DoF=  Degrees  of  Freedom,  Ti=Test  Statistic, 
Wp= 1-a  Quantile 


pdf 

Window 

Center 

n 

DoF 

2\ 

Level  of  Significance 

Beta  MLE 
Beta  MoM 

0° 

121 

2 

2fi3.1736 

13.82 

0.1%>  a 

Wei  bull  MLE 

2 

19.2066 

13.82 

0.1%>  a 

Weibull  MoM 

2 

19.2066 

13.82 

0.1%>  o 

Lognormal 

2 

29.7438 

13.82 

0.1%>  a 

Normal 

2 

5.5702 

5.991 

10%>  a  >5% 

Rayleigh 

3 

409.2893 

16.27 

0.1%>  a 

Rayleigh-y  MLE 

3 

464.3306 

16.27 

0.1%>  a 

Rayleigh-y  MoM 

3 

484.0000 

16.27 

0.1%>  a 

Rayleigh +1 

3 

181.6860 

16.27 

0.1%>  a 

Beta  MLE 

90° 

231 

2 

30.5368 

13.82 

0.1%>  o 

Beta  MoM 

20.6667 

13.82 

0.1%>  a 

Weibull  MLE 

2 

11.8355 

13.82 

0.5%>  a  >0.1% 

Weibull  MoM 

2 

5.9502 

5.991 

10%>  ft  >5% 

Lognormal 

2 

81.1039 

13.82 

0.1%>  a 

Normal 

2 

20.6667 

13.82 

0.1%>  a 

Rayleigh 

3 

489.9307 

16.27 

0.1%>  a 

Rayleigh-y  MLE 

3 

638.4589 

16.27 

0.1%>  ft 

Rayleigh-y  MoM 

3 

369.9740 

16.27 

0.1%>  ft 

Rayleigh +1 

3 

262.7879 

16.27 

0.1%>  ft 

Beta  MLE 

00 

o 

o 

121 

2 

21.9339 

KE221 

0.1%>  ci 

Beta  MoM 

2 

6.1488 

7.378 

5%>  ft  >2.5% 

Weibull  MLE 

2 

9.8678 

10.6 

1%>  a  >0.5% 

Weibull  MoM 

2 

7.1405 

7.378 

5%>  o  >2.5% 

Lognormal 

2 

36.9504 

13.82 

0.1%>  o 

Normal 

2 

6.1488 

7.378 

5%>  ft  >2.5% 

Rayleigh 

3 

175.5702 

16.27 

0.1%>  a 

Rayleigh-y  MLE 

3 

223.4215 

16.27 

0.1%>  o 

Rayleigh-y  MoM 

3 

160.1157 

16.27 

0.1%>  ft 

Rayleigh -|-1 

3 

71.6860 

16.27 

.J 

0.1%>  ft 
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Table  A. 7.  Fighter:  9  =  90°  Kolmogorov  Levels  of  Significance  for  Window  vs.  Window, 
n=number  of  samples,  T^Test  Statistic,  lVr  =  l-a  Quantile 


pdf 

Window 

Center 

n 

H 

w 

v  p 

Level  of  Significance 

Beta  MLE 

0° 

11 

0.2283 

a  >  20% 

Beta  MoM 

0.2792 

B 

a  >  20% 

Weibull  MLE 

0.6355 

1%>  a 

Weibull  MoM 

0.4495 

0.468 

2%> a  >1% 

Lognormal 

0.7276 

0.468 

1%>  a 

Normal 

0.4274 

0.437 

5%>  a  >2.5% 

Rayleigh 

0.5165 

0.468 

1%>  a 

Rayleigh-y  MLE 

0.6349 

0.468 

1%>  a 

Rayleigh-y  MLE 

0.6453 

0.468 

1%>« 

Rayleigh+1 

0.4381 

0.468 

2%> a  >1% 

Beta  MLE 

90° 

21 

0.3018 

impii 

5%> a  >2% 

Beta  MoM 

0.1527 

a  >20% 

Weibull  MLE 

0.5058 

pffl 

1%>  a 

Weibull  MoM 

0.1546 

a  >20% 

Lognormal 

0.7743 

1%>  a 

Normal 

0.1578 

0.226 

a  >20% 

Rayleigh 

0.4013 

0.344 

1%>  a 

Rayleigh-y  MLE 

0.6492 

0.344 

1%>  a 

Rayleigh-y  MOM 

0.6450 

mm 

1%>  a 

Rayleigh+1 

0.2492 

IS 

20%>  a  >10% 

Beta  MLE 

00 

o 

o 

11 

0.2606 

USB 

a  >20% 

Beta  MoM 

0.1964 

B 

a  >20% 

Weibull  MLE 

0.5803 

Kb 

1%>  a 

Weibull  MoM 

0.2460 

0.308 

a  >20% 

Lognormal 

0.7475 

0.468 

1%>  a 

Normal 

0.2248 

0.308 

a  >20% 

Rayleigh 

0.3546 

0.391 

10%>  a  >5% 

Rayleigh-y  MLE 

0.5529 

0.468 

1%>  a 

Rayleigh-y  MoM 

0.5536 

0.468 

1%>  a 

Rayleigh+1 

0.2434 

0.308 

a  >20% 
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Table  A.8.  Fighter:  6  =  90°  Chi-Square  Levels  of  Significance  for  Window  vs.  Window, 
n=number  of  samples,  DoF=  Degrees  of  Freedom,  7\=Test  Statistic,  Wp= 1-a 
Quantile 


pdf 

Window 

Center 

n 

DoF 

m 

wp 

Level  of  Significance 

Beta  MLE 

0° 

11 

1 

2.2727 

2.706 

25%  >  a  >10% 

Beta  MoM 

1 

0.8182 

1.323 

a  >25% 

Weibull  MLE 

1 

4.8182 

5.024 

5%>  a  >2.5% 

Weibull  MoM 

1 

1.5455 

2.706 

25%  >  a  >10% 

Lognormal 

1 

7.3636 

7.879 

1%>  a  >0.5% 

Normal 

1 

2.2727 

2.706 

25%  >  a  >10% 

Rayleigh 

1 

4.4545 

5.024 

5%>  a  >2.5% 

Rayleigh-y  MLE 

1 

4.4545 

5.024 

5%>  a  >2.5% 

Rayleigh-y  MoM 

1 

4.4545 

5.024 

5%>  a  >2.5% 

Rayleigh+1 

1 

2.2727 

2.706 

25%  >  a  >10% 

Beta  MLE 

90° 

21 

1 

13.8571 

10.83 

0.1%>  a 

Beta  MoM 

1 

5.0952 

6.635 

2.5%>  a  >1% 

Weibull  MLE 

1 

17.6667 

10.83 

0.1%>  a 

Weibull  MoM 

1 

3.1905 

3.841 

10% >  a  >5% 

Lognormal 

1 

55.3810 

10.83 

0.1%>  a 

Normal 

1 

3.1905 

3.841 

10%>  a  >5% 

Rayleigh 

2 

14.6190 

13.82 

0.1%>  a 

Rayleigh-y  MLE 

2 

48.1429 

13.82 

0.1%>  a 

Rayleigh-y  MoM 

2 

41.6667 

13.82 

0.1%>  a 

Rayleigh+1 

2 

4.7143 

5.991 

10%>  a  >5% 

Beta  MLE 

t— * 
00 
o 

o 

11 

1 

0.0909 

1.323 

a  >25% 

Beta  MoM 

1 

0.0909 

1.323 

a  >25% 

Weibull  MLE 

1 

3.7273 

3.841 

10%>  a  >5% 

Weibull  MoM 

1 

1.1818 

1.323 

a  >25% 

Lognormal 

1 

7.3636 

7.879 

1%>  a  >0.5% 

Normal 

1 

0.8182 

1.323 

a  >25% 

Rayleigh 

1 

2.2727 

2.706 

25%  >  a  >10% 

Rayleigh-y  MLE 

1 

7.3636 

7.879 

1%>  a  >0.5% 

Rayleigh-y  MoM 

1 

7.3636 

7.879 

1%>  a  >0.5% 

Rayleigh+1 

1 

2.2727 

2.706 

25%  >  a  >10% 

Table  A.9.  Missile:  6  =  90*  Kolmogorov  Levels  of  Significance  for  Window  vs.  Window, 
n=number  of  samples,  T^Test  Statistic,  Wp=\-a  Quantile 


Window 

Center 


Beta  MLE 
Beta  MoM 
Wei  bull  MLE 
Wei  bull  MoM 
Lognormal 
Normal 
Rayleigh 
Rayleigh-y  MLE 
Rayleigh- y  MoM 
Rayleigh+1 


Beta  MLE 
Beta  MoM 
Weibull  MLE 
Wei  bull  MoM 
Lognormal 
Normal 
Rayleigh 
Rayleigh-y  MLE 
Rayleigh-y  MoM 
Rayleigh+1 


Beta  MLE 
Beta  MoM 
Weibull  MLE 
Weibull  MoM 
Lognormal 
Normal 
Rayleigh 
Rayleigh-y  MLE 
Rayleigh-y  MoM 
Rayleigh+1 


n 

T\ 

wp 

Level  of  Significance 

21 

0.1603 

EEia 

a  >20% 

0.1303 

0.226 

a  >20% 

0.2760 

0.287 

10% >  a  >5% 

0.2324 

0.259 

20%>  a  >10% 

0.5055 

0.344 

1%>  a 

0-1919 

0.226 

a  >20% 

0.2048 

0.226 

a  >20% 

0.3221 

0.344 

2%> a  >1% 

0.2998 

HI 

5%> a  >2% 

0.3545 

HH 

1%>  Q 

41 

0.2087 

0.212 

10%>  a  >5% 

0.0716 

a  >20% 

0.3189 

0.255 

1%>  a 

0.1475 

0.167 

a  >20% 

0.6765 

0.255 

1%>  a 

0.1557 

0.167 

a  >20% 

0.3405 

0.255 

1%>  a 

0.5720 

0.255 

1%>  a 

0.5796 

0.255 

1%>  a 

0.2365 

0.237 

5%>  a  >2% 

21 

0.2377 

20%  >  a  >10% 

0.1214 

0.226 

a  >20% 

0.2715 

0.287 

10% >  a  >5% 

0.1956 

0.226 

a  >20% 

0.3850 

0.344 

1%>  a 

0.2520 

20%  >  a  >10% 

0.2887 

0.321 

5%> a  >2% 

0.2947 

0.321 

2%> a  >1% 

0.2454 

0.259 

20%>  a  >10% 

0.3430 

0.344 

ss 

f—4 

A 

a 

A 

IN 

A-iU 


Table  A. 10.  Missile:  6  =  90°  Chi-Square  Levels  of  Significance  for  Window  vs.  Window, 
n=number  of  samples,  DoF=  Degrees  of  Freedom,  7\=Test  Statistic,  Wp= 1- 
a  Quantile 


pdf 

Window 

Center 

T, 

wp 

Level  of  Significance 

Beta  MLE 

0° 

21 

1 

4.3333 

5.024 

5%>  a  >2.5% 

Beta  MoM 

1 

0.5283 

1.323 

a  >25% 

Weibull  MLE 

1 

6.6190 

5.0246.635 

2.5%>  a  >1% 

Weibull  MoM 

1 

7.0000 

7.879 

1%>  a  >0.5% 

Lognormal 

1 

18.4286 

10.83 

0.1%>  a 

Normal 

1 

4.3333 

5.024 

5%>  a  >2.5% 

Rayleigh 

2 

8.1429 

9.210 

2.5%>  a  >1% 

Rayleigh-y  MLE 

2 

5-8571 

5.991 

10%>  a  >5% 

Rayleigh-y  MoM 

2 

6.2381 

7.378 

5%>  a  >2.5% 

Rayleigh+1 

2 

15.000 

13.82 

0.1%>  a 

Beta  MLE 

90° 

41 

2 

59.8780 

13.82 

0.1%>  a 

Beta  MoM 

2 

1.5610 

2.773 

a  >25% 

Weibull  MLE 

2 

16.4310 

4.605 

0.1%>  a 

Weibull  MoM 

2 

3.5122 

4.605 

25%>  a  >10% 

Lognormal 

2 

104.4878 

13.82 

0.1%>  a 

Normal 

2 

3.5122 

4.605 

25%>  a  >10% 

Rayleigh 

3 

21.5610 

16.27 

0.1%>  a 

Rayleigh-y 

3 

75.9512 

16.27 

0.1%>  a 

Rayleigh-y 

3 

82.7805 

16.27 

0.1%>  a 

Rayleigh+1 

3 

14.4878 

16.27 

0.5%>  a  >0.1% 

Beta  MLE 

180° 

21 

1 

5.4762 

6.635 

2.5%>  a  >1% 

Beta  MoM 

1 

1.6667 

2.706 

25%>  a  >10% 

Weibull  MLE 

1 

5.8571 

6.635 

2.5%>  a  >1% 

Weibull  MoM 

1 

6.2381 

6.635 

2.5%>  a  >1% 

Lognormal 

1 

12.7143 

10.83 

0.1%>  a 

Normal 

1 

6.2381 

6.635 

2.5%>  a  >1% 

Rayleigh 

2 

8.5238 

9.210 

2.5% >  a  >1% 

Rayleigh-y 

2 

12.3333 

13.82 

0.5%>  a  >0.1% 

Rayleigh-y 

2 

6.2381 

7.378 

5%>  a  >2.5% 

Rayleigh+1 

2 

8.9048 

9.210 

2.5%>  a  >1 

A- 11 


Table  A.  11.  Fighter:  0  —  90°  ±  5°  Kolmogorov  Levels  of  Significance  for  Window  vs. 
Window,  n=number  of  samples,  Ti=Test  Statistic,  Wp  =  l-a  Quantile 


pdf 

Window 

Center 

n 

■a 

wp 

Level  of  Significance 

Beta  MLE 

0° 

121 

0.1816 

1%>  a 

Beta  MoM 

0.1074 

20%>  a  >10% 

Weibull  MLE 

0.1943 

B|pl| 

1%>  a 

Weibull  MoM 

0.1142 

10%a  >5% 

Lognormal 

0.5295 

0.148 

1%>  a 

Normal 

0.1236 

0.124 

10%>  a  >5% 

Rayleigh 

0.3953 

0.148 

1%>  a 

Rayleigh-y  MLE 

0.5850 

0.148 

1%>  a 

Rayleigh-y  MoM 

0.5735 

0.148 

1%>  a 

Rayleigh+1 

0.2376 

0.148 

1%>  a 

Beta  MLE 

90° 

231 

0.1576 

jBTjf 

1%>  a 

Beta  MoM 

0.0780 

20%  >  a  >10% 

Weibull  MLE 

0.1244 

1%>  a 

Weibull  MoM 

0.0758 

20%  >  a  >10% 

Lognormal 

0.4818 

1%>  a 

Normal 

0.0898 

10%>  a  >5% 

Rayleigh 

0.3877 

fill 

1%>  a 

Rayleigh-y  MLE 

0.6199 

0.107 

1%>  a 

Rayleigh-y  MoM 

0.6046 

0.107 

1%>  a 

Rayleigh+1 

0.2611 

0.107 

1%>  a 

Beta  MLE 

180° 

121 

0.1296 

mmMSM 

5%>  a  >2% 

Beta  MoM 

0.0679 

a  >20% 

Weibull  MLE 

0.1862 

mm 

1%>  a 

Weibull  MoM 

0.1199 

0.124 

10%>  a  >5% 

Lognormal 

0.5354 

0.148 

a  >1% 

Normal 

0.1197 

0.124 

10%>  a  >5% 

Rayleigh 

0.3212 

0.148 

1%>  a 

Rayleigh-y  MLE 

0.5675 

0.148 

1%>  Q 

Rayleigh-y  MoM 

0.5674 

0.148 

1%>  a 

Rayleigh+1 

0.2004 

0.148 

1%>  a 
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Table  A. 12.  Fighter:  6  —  90°  ±5°  Chi-Square  Levels  of  Significance  for  Window  vs.  Win¬ 
dow,  n=number  of  samples,  DoF=  Degrees  of  Freedom,  7\=Test  Statistic, 
Wp=l-a  Quantile 


pdf 

Window 

Center 

n 

DoF 

Ti 

wp 

Level  of  Significance  \ 

Beta  MLE 

0° 

121 

2 

50.4463 

13.82 

0.1%>  a 

Beta  MoM 

2 

9.2066 

9.21 

2.5%>  a  >1% 

Weibull  MLE 

2 

21.6860 

13.82 

0.1%>  a 

Weibull  MoM 

2 

6.5620 

7.378 

5%>  a  >2.5% 

Lognormal 

2 

149.3884 

13.82 

0.1%>  a 

Normal 

2 

6.5620 

7.378 

5%>  a  >2.5% 

Rayleigh 

3 

78.2975 

16.27 

0.1%>  a 

Rayleigh-y  MLE 

3 

252.6777 

16.27 

0.1%>  a 

Rayleigh-y  MoM 

3 

235.0744 

16.27 

0.1%>  a 

Rayleigh-}- 1 

3 

22.0165 

16.27 

0.1%>  a 

Beta  MLE 

90° 

231 

2 

60.7965 

0.1%>  a 

Beta  MoM 

2 

15.7749 

13.82 

0.1%>  a 

Weibull  MLE 

2 

23.9567 

13.82 

0.1%>  a 

Weibull  MoM 

2 

15.6883 

13.82 

0.1%>  a 

Lognormal 

2 

278.7662 

13.82 

0.1%>  ft 

Normal 

2 

16.1645 

13.82 

0.1%>  a 

Rayleigh 

3 

234.5195 

16.27 

0.1%>  « 

Rayleigh-y  MLE 

3 

527.5065 

16.27 

0.1%>  a 

Rayleigh-y  MoM 

3 

484.2597 

16.27 

0.1%>  a 

Rayleigh-}- 1 

3 

118.2857 

16.27 

0.1%>  a 

Beta  MLE 

o 

O 

00 

rH 

121 

2 

24.3306 

13.82 

0.1%>  ft 

Beta  MoM 

2 

8.7107 

9.210 

2.5%>  a  >1% 

Weibull  MLE 

2 

16.2314 

13.82 

0.1%>  a 

Weibull  MoM 

2 

8.2975 

9.210 

2.5%>  «  >1% 

Lognormal 

2 

188.5950 

13.82 

0.1%>  a 

Normal 

2 

8.8760 

9.210 

2.5%>  ft  >1% 

Rayleigh 

3 

50.1983 

16.27 

0.1%>  a 

Rayleigh-y  MLE 

3 

211.5207 

16.27 

0.1%>  « 

Rayleigh-y  MoM 

3 

211.5207 

16.27 

0.1%>  ft 

Rayleigh-}- 1 

3 

20.6116 

16.27 

0.1%>  ft 
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Appendix  B.  ASPECT  RCS  Modeling  Program 


This  appendix  presents  the  APECT  RCS  modeling  program  along  with  the  ’new 
ton’,’bmleVchi’,  and  ’interp’  subroutines.  The  ASPECT  MATLAB  code  inputs  a  2 
dimensional  array  and  determines  the  optimum  pdf  model. 


^ASPECT  RCS  Modeling  Program 
c  =  0; 

while  c~  =  1 
clear 

s*sinput( 'Input  number  of  file  to  be  used.lscttm, 
2=miccilel80,  3=  wid,  4  =  mU«ile‘); 
if  z=  3=  1 
load  cttm 
aaittm; 
elceif  a  ==2 
load  miscilelfiO 
s=miccilel80; 
elceif  i==3 
load  wid 
z=wid; 
eUeif  z==4 
load  missile 
z= missile; 
end 
k=0; 

t  =  input(’Input  number  of  separate  section*  to 
be  included  in  fulbrange  data. 

You  will  be  prompted  to  input  the  beginning 
and  ending  «ample  of  each  section.’) 
for  os  1:1 

beginl=input(’Input  beginning  length  sample’) 
beginw=input('Input  beginning  width  cample*) 
endl=input(’Input  ending  length  cample*) 
endw=input('Input  ending  width  cample') 
for  i  =  beginl:endl 
for  j=beginw:endw 
kssk  +  l ; 
xl(k)=j(i,j); 
end 
end 
end 

mescage=’Data  ic  now  linearised.  Window  will 
be  a  cingle  row  table  with  length  equal  to 
aspect  length  multiplied  by  aspect  width.' 
list  l’; 

c**(sl-min(*l  ))/abc(max(sl)-min(zl)); 

dataak 

j*0; 

y=0; 

losinput(’Input  beginning  window  cample') 
h>sinput('Inpul  ending  window  cample’) 
magc=hi-lo+l 


for  ialrdata 
if  i;Io — iihi 
cO<»)=0; 
end 

if  iislokiishi 

i-i+is 

tO(i)-»(i); 

*Om(j  )=•(>); 
if  c0m(j)c=iO 
c0m(i)=10‘(-«); 
elceif  i0m(j)3sl 

c0m(j)sel-10‘(*6); 

end 

ncOm(j  )=  l*cOra(  j ); 

*0ml(j)=c0m(j)‘2; 
c0m4(  j  )s«Om(j  )  *  4 ; 
»0mgl(j)=»Om(j,-(l/m»g); 

)sniOm(j  )  ‘  ( 1  /m»g); 

end 

if  iilo— i*hi 
*Ou(i)=cO; 
end 

if  i^slolciiahi 

y=y+i; 

cOu(i)=.(i); 

cOum(y)=c(i); 

end 

end 

Gl  =  prod(*Omgl)%  for  beta  mle  alpha 
G2sprod(f0mg2)%  for  beta  mle  beta 
c01=cO+10*(-6); 

KMOMENT  GENERATION 

mamean(cOm)%»ample  mean 
m2amean(c0ml)9lcample  second  moment 
m4=mean(c0m4)%sample  4th  moment 
M2=(ctd(e0m))‘2fc  cample  variance 
var=m2'm '  2%  cample  variance 
d  =  m  ‘  2/ m2%  for  weibutl  mom  interp 
medsmedian(cOm)Kcample  median 

(or  i«lo:hiH  lognormal  parameter  ect 
lmean(i)z]og(c01(i)); 
end 

lmeanx«um(lmean)/mag  %  log  mean 
for  islorhi 

Ivar(i)=s(log(s01(i))-lmean)'2; 
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end 

|var=*um(lvar)/mag  %  log  variance 


pause; 
plot(x.W); 
WMLEP  =  [x’  W»]; 


94BETA,  WEIBULL,  AND  RAYLEIGH 
PARAMETER  DETERMINATION 

rl=(3*m2/4)‘(0.5)J&  Aayleigh-fl  parameter 

r=(m2/2)‘(0.5)%  Rayleigh  parameter 

rymle  =  r9i  Rayleigh-y=x*2  mle  parameter 

rymom=(m2/8)*(0.25)%  Ry  mom  parameter 

a=m*(m-m2)/(var)%a  beta  estimate  using  method  of  moments 

b=(l-m)*(m-m2)/(var)fcb  beta  mem 

interp  %  Determine  Weibull  MoM  Parameters 

k=input('If  Weibull  mle  is  desired  press  1,  if  mom  press  RETURN’); 
if  keel 

newton  %  Determine  Weibull  MLE  parameters 
end 

kl=input(’If  Beta  mle  is  desired  press  1,  if  mom  press  RETURN’); 
if  kl  =  =  l 

bmle  %  Determine  Beta  MLE  parameters 
val = (  gam  ma(  a-f  b  )  )  /  ( gamroa( a)* gamma(  b ) ) %  beta 
end 

%  PDF  FORMATION 
for  is  1:200 
x(i)=i/200; 

B(i)=»»l*(x(i)(»-l))-((l-«(i))-(b-l)); 
if  B(i  )==inf 
B(i)=B(i-l)*4; 
end 

W(»)=(al/be)*(x(i)/be)‘(al-l)*exp(-(x(i)/be)*al); 

R(  i  )r=  2*(  x(  i )/  m2  )*ex  p(-(  x(i )  ‘  2  )/(  m2)); 

Rymle(i)=(l/(2*(rymle*2)))*cxp(-(x(i)*2)/(2*(rymle  *  2))); 

Rymom(i)=(l/(2*(rymom‘2)))*exp(-(x(i)‘2)/(2*(rymom‘2))); 

Rl(i)=(9*x(i)-3/(2»(rl)*4))-’exp((.3*x(i)-2)/(2*(rl)*2)); 

L(i)=((»(ir*qr»(lvar)*(2*pi)-(l/2))-(-l)rexp(-(log(x(i)) 

-lmean)‘(2)/(2*Ivar)); 

N(i)=((l/(2*pi*var»*(l/2))*exp(.((x(i>m)*2)/(2*var)); 

end 

hold  off; 

%hold  on, 
i  =0.025:0.05:0.975; 

[Data,  uj=hist(s0am,i); 
hist(s0um,i ); 
for  i  =  l:200 
for  j  =  l:20 

if  (i/200)i=l0*j/200dr(i/200)il0*(i*l)/200 
Data2(i)=Data(j); 
end 

if  i  ==200 
Data2(i)=0; 
end 
end 
end 

plot(x,Data2); 

DataPsfx’  Data2’]; 

save  DataP.dat  DataP  /ascii; 


save  WMLEP.dat  WMLEP  /ascii; 
WMoMP  =  [x’  W’J; 
save  WMoMP.dat  WMoMP  /ascii; 
pause; 

%hold  on; 
plot(x,L); 

LognormalP=(x’  L'J; 

save  LognormalP.dat  LognormalP  /ascii; 
pause; 

9thold  on; 
plot(x,B); 

%title(’Beta  MLE  pdf:  0  deg’); 

0MoMP=[x’  B’J; 

save  BMoMP.dat  BMoMP  /ascii; 

BMLEP  =  [x’  B’J; 

save  BMLEP.dat  BMLEP  /ascii;pause; 

9th  old  on; 

plot(x.N); 

NormalP  =  (x’  N’J; 

save  NormalP.dat  Normal?  / ascii; 

pause; 

fthold  on; 

plot(x.R); 

RayleighP  =  [x’  R'J; 

save  RayleighP.dat  RayleighP  /ascii; 

pause; 

9£hold  on, 
plot(x.Rymle); 

%hold  on; 

Rayleighy P  =  [x’  Rymle’}; 

save  RayleighyP.dat  RayleighyP  /ascii; 

pause; 

plot(x,Rymom); 

%hoId  on; 

RayleighyoP=[x’  Rymom’J; 

save  RayleighyoP.dat  RayleighyoP  /ascii; 

pause; 

96ho!d  on; 
plot(x,Rl); 

RayleighlP=(x’  Rl’J; 

save  RayleighlP.dat  Rayleigh  IP  /  ascii; 

pause; 

[B  max ,  i  B \ = max(  B  ) 

[W  max  ,i  W]=max(  W) 

[Lmax,iL]=max(L) 

(Nmax,iN]=max(N) 

(Rmax,iR]=max(R) 
(Rymlemax,iRymle]=max(  Rymle) 

[Ry  mommax  ,i  Ry  mom)  =  max(Ry  mom) 
[Rlmax,iRl]=max(Rl) 

HCDF  Formation 

j=0; 

for  i  =  l:200 
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x(»)=i/200; 

sP(i)=0; 
for  j  =  1 :  mag 
if  *Oum(j  )==0<5di==l 
*P(i)SS«P(i)+l; 
end 

if  x(i)i(«0mji(j>l/(400))4x(i)i=(»0ttni(j)+l/(400)) 

end 

end 

end 

for  i=l:200 
if  i  =  =  l 
sPl(»)=sP(i); 

BPl(i)«B(i); 

WPl(i)=W(i); 

RPl(i)«R(i); 

RymlePl(i)=Rymle(i); 

RymomPl(i)«Rymom(i); 

RiPl(i)=Rl(i); 

LP 1  ( i )  =  L(  i ) ; 

NPl(i)ssN(i); 

else 

sPl(i)=sP(i)-f  *Pl(i-l); 

BPl(i)=B(i)+BPl(i-l); 

WPl(i)=W(i)+WPl(i-l); 

RPl(i)=R(i)-f  RPl(i'l); 
RymleFl(i)=Rymle(i)+RymlePl(i*l); 
RymomPl(i)=Ryraom(i)-f  RymomPl(i-l); 
RlPl(i)=Rl(i)+RlPl(i-l); 

LPl(i)=L(i)+LPl(i-l); 

NPl(i)=N(i)+NPl(i-l); 

end 

end 

sPl  =sPl/max(sPl); 

BPl=BPl/max(BPl); 

WPJ  =  WPl/max(WPl); 

RP 1  =  RP 1  /  max(RP  1 ); 

Ry  mleP  1 =Ry  mleP  1  /m»x(  Ry  mleP  1 ); 
RymomPlsRymomPl/max(RymomPl); 

RlPl  =  RlPl/max(RlPl); 

LPl  =  LPl/max(LPl); 

NPl  =  NPl/max(NPl); 

plot(x,sPl,x,RPl); 
title(‘R*y!eigh  cdf:  0  deg’); 

R*yleighC=lx’  RPl’V, 

save  RayleighC.dat  RayleighC  /ascii; 

Sprint; 

pause; 

plot(x,*Pl,X,RlPl); 

title(  ’  Rayleigh -fl  cdf:  0  deg*); 

%print; 

Rayleigh  1Csb(x*  RlPl’l; 

save  RayleighlC.dat  RayieighlC  /ascii; 

pause; 

p}ot(x,*Pl,x,Rym1ePl); 
title(’RayIeigh-y  mle  cdf:  0  deg*); 


Rprint; 

RayleighyCss[x’  RymlePl’); 

save  RayleighyC.dat  RayleighyC  /ascii; 

pause; 

plot(x,sPl,x,RymomPl); 
title(>Rayleigh-y  mom  cdf:  0  deg*); 
Kprint; 

Rayleigh  yoC=lx’  RyraomPl*]; 

save  RayleighyoC.dat  RayleighyoC  /ascii; 

pause; 

plot(x,sPl,x,LPl); 

ti tle( 'Lognormal  cdf:  0  deg’); 

Kprint; 

LognormalC=[x*  LPl’); 

save  LognormaiC.dat  LognormalC  / ascii; 

pause; 

plol(x,sPl,X,NPl); 
title( ’Normal  cdf:  0  deg*); 

%print; 

NormalC=[x*  NPl’J; 

save  NormalC.dal  NormalC  / ascii; 

pause; 

plot(x,sPl,x,BPl); 

title('Beta  MoM  cdf:  0  deg*); 

BMoMC  =  [x'  BPl’J; 

save  BMoMC.dat  BMoMC  /ascii; 

BMLEC=[x’  BPl’J; 

save  BMLEC.dat  BMLEC  /a*cii;ftprint; 
pause; 

plot(x,sPl,x,WPl); 

title(’Weibuil  MoM  cdf:  0  deg’); 

WMoMC=[x’  WPl’J; 

save  WMoMC.dal  WMoMC  /ascii; 

WMIEC  =  [x’  WPl’J; 

save  WMLEC.dat  WMLEC  /ascii; 

DaiaC=(x’  sPl’J; 

save  DataC-dat  DataC  /ascii; 

Rprint; 

pause; 

%Goodne»s-Of’Fit 

for  i=  1.200 

gfb(i)=abs(sPl(i)-BPl(i)); 

gfw(i)=abs(*Pl(i)-WPl(i)); 

gfr(i)=abs(sPl(i)-RPl(i)); 

gfn(i)=abs^Pl(i>RlPH»»; 

gfrymle(i)=abs(sPl(i)-RymlePl(i)); 

gfrymom(i)=abs(sPl(i)-RymomPl(i)); 

gfl(i)=abs(sPl(i)-LPl(i)); 

gfn(i)=abs(sPl(i)-NPl(*)); 

end 

%  Kolmogorov  Test  Statistic 

kgfbBmax(gfb) 
kgfw  =  max(gfw) 
kgfrsmax(gfr) 
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kgfrl  =max(gfrl ) 
kgfrymle=max(gfrymle) 
kgfrymom=max(gfrymom ) 
kgfl  =  max(gfl) 
kgfnsmax(gfn) 

pdf=min(kgfb,kgfw); 
pdf=min(pdf,kgfr); 
pdf=min(  pdf.kgfrl ); 
pdf=min(pdf,kgfrymle); 
pdfs  min(  pdf  ,kgfry  mom); 
pdf=min(pdf,kgfl); 
pdf=min(pdf,kgfn); 
if  pdf—  =  kgfb 

messages  'Bet*,  is  the  optimum  pdf,  alphas' 

a 

measages’beta=’ 

b 

elseif  pdf=  =  kgfw 

messages ’Weibull  is  the  optimum  pdf,  alphas’ 
al 

messages ’betas’ 
be 

elseif  pdf=  =  kgfr 

messages ’Rayleigh  is  the  optimum  pdf,  alphas’ 

(m2/2)‘0.5 
elseif  pdfsskgfrl 

messages’One-Dominant-Plus-Rayleigh  is  the  optimum  pdf, 
alphas’rl 

elseif  pdfsskgfrymle 

messages’ Rayleigh- y  mle  is  the  optimum  pdf,  alphas’ 
rymle 

elseif  pdfsskgfrymom 

message  s’ Ray  lei  gh-y  mom  is  the  optimum  pdf,  alpha=' 
rymom 

elseif  pdfsskgfl 

messages’Lognormal  is  the  optimum  pdf,  sigmas’ 
sqrt(lvar) 

messages ’mus’ 
lmean 

elseif  pdfsskgfn 

messages ’Nomnal  is  the  optimum  pdf,  sigmas' 
var 

messages  ’mu  = ' 
m 
end 

messages’Enter  RETURN  for  Chi-Square  GOF’ 
pause; 

%  Chi-Square  Goodness  of  Fit 
chi 

csinput(’Enter  i  to  quit,  enter  another  number  to  continue'); 
ifc  =  =[] 
csO; 
end 
end 


ftINTERP:  Weibull  MoM  Parameter  Formation 


load  tablewl.mat 
for  i  =  l:100 
diff(i)sd-  table  wl(i); 
end 
k=0; 

for  isl:100 
if  diff(i)»=0 

ksk  +  1; 

diffn(k)sdiff(i); 

end 

end 

(o,p]smin(abs(diifn)); 

y2=(p-l)/100; 

yl=(p)/100; 

jclstablewl(p); 

*2s  table  wl(p-f  1); 

ys(x2-d)*yl/(x2-xl)+(d-xl)*y2/(x2-xl ); 
almomsl/y 

bemomsm/gamma(lHhy) 

al=almom 

besbemom 


*NEWTON:  Weibull  MLE  Parameter  Formation 

j=0; 

almlesinput(' Input  beginning  iteration  (al  suggested)') 
for  islo:hi 

j=j+i; 

sn0m(j)=s(i)+10’(-6); 
ln*n01(j  )slog(sn0m(j)); 
sn01(j)ssn0m(j); 
sn02(j)sl-sn0m(j); 

end 

sum2=sum(lasn0l); 

suml  s  l ; 

»um3=l; 

jsO; 

while  (suml/sum3-l/almle)"=sum2/mag 
for  i= 1  .mag 

Insn0(i)s(an0m(i)'almle)*log(sn0m(i)); 
sn0c0(i)ssn0m(i)'  almle; 
Insn02(i)s(sn0m(i)‘almle)*(log(sn0m(i)))*2; 
end 

snmlssnro(lnsnO); 
sum3ssum(sn0c0); 
sum4  =  sum(lnsn02); 

almlesalmle-f  (sum2/mag-f  l/almle-(suml/sum3))/ 
((l/alm!e)‘2+(sum3*sum4-(suml)*2)/(*um3)*2) 
bemles(sum3/mag)*(l /almle) 

j*j+i; 

if  js  =  2Q 
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break; 

end 

end 

al  =  almle 
be  =  bemle 

%BMLE:  Beta  MLG  Parameter  Formation 

load  betat.mat; 
if  Gl  =  =  G2&Gl  =  =0 
a=0 
b=0 
end 

if  GUG2 
Gla=G2; 

G2a=Gl; 

elaeif  Gl;  =  G2&G2'=0 
Gla=Gl; 

G2a=G2; 

end 

if  G2"  =  0 
for  i=2:122 

if  Glai  =  betat(i,l)^Gla^betat(i*l,l) 
Glxl  =  betal(i-l,l ); 

Glx2=betat(i,l); 

end 

if  Gla==0 
Glxl=0-0); 

Glx2=005, 

end 

end 

u=0; 

k=0; 

for  i=2:122 

if  betat(i,l)==Glxl 
u=u  +  l; 

»(»)=i; 

if  G2ai=betat(ip2)&G2aj.betal(i«l>2) 
G2xl  l=betat(i-l,2); 
G2xl2=betat(i,2); 
al  l  =  betat(i- 1,3); 
al2=betat(i,3); 
bl  l  =  betat(i>l,4); 
b!2=betat(i,4); 
end 
end 

if  betat(i,l)=  =  Glx2 
k=k  +  l; 
h(k)=i; 

if  G2a;  =  belat( i ,2)4rG2a(i,  =  be t at(i- 1,2) 
G2x21  =  betat(i-1,2); 
G2x22=b«Ut(i,2); 
a21=betat(i- 1 ,3); 
a22=betat(i,3); 
b21  =  betat(i- 1 ,4); 
b22=betat(i,4); 


end 

end 

end 

if  betat(v(u),2)iG2a 

G2xll  =  betat(v(u)-1,2); 
G2xl2=betat(v(u),2); 
all=betat(v(n)-l,3); 
al2=b«tat(v(n).3); 
bll  =  betat(v(u)-l,4); 
bl2=betat(v(u),4); 
end 

if  betat(h(k),2);G2* 

G2x21=betat(h(k)-1,2); 

G2x22=betat(k(k),2); 

a21  =  betat(h(k>l,3); 

a22=betat(h(k),3); 

b21  =  betat(h(k>l,4); 

b22=betat(h(k),4); 

xl=G2xll; 

x2=G2xl2; 

x  =  G2a; 

yl=all; 

y2=al2; 

y=(x2-x)*yl/(x2-xl)+(x-xl  )*y2/(x2-xl); 

al=y 

yl=bll; 

y2=bl2; 

y=(x2-x)*yl/(x2-xl)+(x-xl)*y2/(x2-xl); 

bl=y 

xl  =  G2x21; 
x2=G2x22; 
yl=a21; 
y2=a22; 

y=(x-xl)*y2/(x2-xl)-(x-x2)*yl/(x2-xl); 

a2=y 

yl=b21; 

y2=b22; 

y=(x-xl)*y2/(x2-xl)-(x-x2)*yl/(x2-xl); 

b2=y 

xlsGlxl; 

k2=G1x2; 

x  =  Gla; 

yl=al; 

y2=a2; 

y=(x2-x)*yl/(x2-xl)+(x-xl  )*y2/(x2-xl); 

a=y 

yl=bl; 

y2=b2; 

y=(x2-x)*yl/(x2-xl)+(x*xl  )*y2/(x2-xl); 
b=y 

elieif  Gla==0 
xl=G2xll; 
x2=G2xl2; 
x  =  G2a; 
yl=al  1; 
y2=al2; 

y=(x2-x)*yl/(x2-xl)+(x-xl  )*y2/(x2-xl); 
al  =  y 
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ylabll; 

y2=bl2; 

y=(x2-x)*yl/(*2-xl)+(x-xl)*y2/(x2-xl); 

bl  =  y 

xl=G2x21; 
x2=G2x22; 
y  1=4.21; 
y2=a22; 

y  =  (x2-x)*y  l/(x2-xl)+(x-xl  )*y2/(x2-xl); 

a2=y 

yl  =  b21; 

y2=b22; 

y=(x2-x)*yl/(x2-xl)+(x-xl)*y2/(x2-xl); 

b2=y 

xl=Glxl; 

x2  =  Glx2; 

x  =  Gla; 

y  1=4.1; 

y2  =  a2; 

y  =  (x2-x)*yl/(x2-xl)-(xl-x)"y2/(x2-xl); 

*=y 

yl  =  bl; 
y2=b2; 

y=(x2-x)*yl  /(x2-xl)-(xl-x)*y2/(x2-xl); 
b=y 
else 

xl  =  G2xll; 
x2=G2xl2; 
x=G2a; 
yl=all ; 
y2=a!2; 

y=(*2-x)*yl/(x2-xi)+(x-xl)*y2/(x2-xl); 

*l=y 

yl=bll; 

y2=bl2; 

y=(*2-x)*yl/(x2-xl)+(x-xl  )*y2/(x2-xl); 
bl=y 

xl=G2x21; 

x2=G2x22; 

yl=*21; 

y2=»22; 

y=(x2-x)*y  l/(x2-xl)+(x-xl  )*y2/(x2-xl ); 

a2=y 

yl  =  b21; 

y2=b22; 

y=(x2-x)*yl/(x2-xl)-f(x-xl  ),*y2/(x2-xl); 

b2=y 

xl=Glxl; 

x2=Glx2; 

xsQli; 

yl=al; 

y2=»2; 

y=(x2-x)*yl/(x2-xl)+(x-xl)*y2/(x2-xl); 

*=y 

yl  =  bl; 

y2=b2; 

y=(x2-x)*y  l/(x2-xl)+(x-xl  )*y2/(x2-xl); 

b=y 


end 

if  GUG2 
bb  =  a; 
aa  =  b; 
b  =  bb 

4  =  44. 

end 

end 

%CHI:  Determine  Chi-Square  Goodness- of- Fit 

if  ma;i=10&magtl5 
for  i  =  1:200 

bc(i)  =  abs(  BP  l(i)-0.5); 
wc(i)=ab*(WPl(i)-0.$); 
rc(i)=abf<RPl(i)-0.5); 
rymlec(i)=ab«(RymlePl(i)-0.5); 
rymomc(i)=ab«(RymomPl(i)-0.5); 
rlc(»)=abs(RlPl(i)-0.5); 
lc(i)=abs(LPl(i)-0-5); 
nc(i)=abs(NPl(i)-0.5); 
end 

[bc,*amb]  =  min(bc); 

(wc,jimw]=min(wc); 

[rc,samr]  =  min(rc); 
[rymtec,samrymle]=min(rymlec); 
[rymomc,iamrymom]=min(  rymomc); 
[rlc,samrl]  =  min(rlc); 

[lc,samt]  =  min(lc); 

[nc,«amn}=min(nc); 

NBl=0; 

NB2=0; 

N  Wl  =0; 

N  W2=0; 

N  Rl  =0; 

NR2=0; 

NRymlel=0; 

NRymle2=0; 

NRymoml=0; 

NRymom2=0; 

NRl  1=0; 

NR12=0; 

NL1=0; 

NL2=0; 

NN1=0; 

NN2=0; 

jOuml  =tort(s0um); 
for  i  =  l:mag 
if  s0uml(i)i=samb/200 
NB1=NB1  +  1; 
end 

if  i0uml(i);samb/200 
NB2*NB2+1; 
end 

if  s0uml(i)i=«amb/200 
NWUNWl  +  1; 
end 
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if  •0uml(i);a&mw/200 
NW2=NW2+1; 
end 

if  a0uml(i);=a&mr/200 
NRl  =NRl  +  1; 
end 

if  *0uml(i)i,a*.mr/200 
NR2=NR2+1; 
end 

if  s0uml(i)j=*4mrymIe/2OO 
NRymlel  =  NRymIel  +  l; 
end 

if  «0uml(i)i«4mrymle/200 
NRymle2  =  NRymle2-f  1; 
end 

if  s0uml(i);=44mrymoin/200 
NRymoml  rNRymoml  +  1 ; 
end 

if  c0uml(i}(.aAmrymom/200 
NRymom2=NRymom2  +  l ; 
end 

if  •Oum](i)i=*arnrl /2O0 
NRllxNRll  +  1; 
end 

if  *Ouml(i )i,anmrl/200 
NR12=NR12+1; 
end 

if  a0uml(i)i=sanml/200 
NLlsNLl  +  1; 
end 

if  s0uml(i)i*4ml/200 
NL2=NL2+|; 
end 

if  s0um](i)i=a*mn/200 
NNl  =  NNl  +  l; 
end 

if  «0uml(i)t**mn/200 
NN2=NN2+1; 
end 
end 

TSB1  =s((  NB  1.0.5*mng) '  2)/(0.5*in*f  ); 
TSB2=((NB2-0.5*mng)  ‘  2)/(0.5*m4g); 

TSBsTSBl +TSB2 

TSW1=((N  Wl-0.5*m4g)*2)/(0.$*m4g); 

TS  W2=((  N  W2-0.S*m*g)*  2)/(0.5*mng); 

TSWssTS  W1+TSW2 

TSRl  =((  NRl-0.5"mog)‘2)/(0.S*m4g); 

TSR2=((NR2-0.5*m*g)'2)/(0.5*m»g); 

TSRsTSRl  +TSR2 

TSRymlel  =:((  NRymlel-0.5*m4g)‘  2)/(0.5*mng), 
TSRymle2=((NRymle2-0.5*m4g)'2)/(Q  5"m4g); 
TSRymle  =  TSRymlel+TSRymle2 
TSRymoml=((NRymoml-0.b*m4g)‘2)/(0.5*m4g); 
TSRymom2=((NRymom2-0.5*mAg)’2)/(0 
TSRymom  =  TSRymoml+TSRymom2 
TSRl  1=((  NRl  l-0.5*m4g)‘2)/(0.5*m4g); 
TSR12*((NR12-0.5»ra*g)*2)/(0  S*m*g); 

TSRl  =TSR1 1-+TSR12 


TSLl=((NLl-0.5*m4g)  *2)/(G.5*m4g); 
TSL2=((NL2'0.5*m4g)'2)/(0.5*m*g); 
TSL=TSLl+TSL2 

TSNl={(NNl*Q.b*m4g)*2)/(0.&*m4g); 
TSN2=((NN2-0.6"m*g)'2)/(0-5*m4g); 
TSN  =TSN 1 +TSN2 
end 

Xpnusc; 

if  15i  =  mnglem4gi20 
for  »=1:200 

bcl(i)=*b*(BPl(i)-l/3); 
bc2(i)=4ba(BPl(i)-2/3); 
wcl(i)  =  4ba(WPl(i>l/3), 
wc2(i)=4ba(WPl(i)-2/3); 
rcl{  i)=*b*(RPl(i)-l/3); 
rc2(i)=4ba(RPl(i)-2/3); 
rymomcl(i)=*ba(RymomPl(i)-l  /3); 
rymomc2(i)  =  4ba(RymomPl(i)*2/3), 
rymlec  l(i  )  =  *b«(  RymleP  1  ( i )- 1  /3); 
rymlec2(i)=4b«(Ry  mlePl(i)-2/3); 
rlcl(i)=4ba(RlPl(i)-l/3); 
rle2(i)=4b«<RlPl(i)-2/3)1 
lcl(i)=4ba(LPl(i)-l/3); 
lc2(i)=4b*(LPl(»>-2/3); 
ncl(i)=4ba(NPl(i)-l/3); 
nc2(»)=4b»<NPl(i)-2/3)t 
end 

[bcl,«4mbl]  =  min(bcl); 
(wcl,nmwl]=min(wcl ); 

[rc  1  ,**mrlj  =  min(  rc  1 ); 

[rymlec  1  ,«4mrymiel]=  min(  rymlec  1 ); 

[ry  momc  1  ,«mrymoral)=min(rjrmomc  1 ); 
[rlcl  ,t4mrll]=mim(rlcl); 

(lc  1  ,»4mll)  =  min(lcl ); 
[nc],*4mnl}=tnin(ncl); 

[bc2,*4mb2]=  min(bc2); 

(wc2,*4mw2]  =  min(wc2), 

[rc2,*»mr2]  =  min(  rc2); 

(ryreomc2,«4mrymom2]=min(rymomc2); 

(rymlec2,*4mrymle2]=min(rymlec2); 

{rlc2,<4inrl2]=min(rlc2); 

(Ic2,«4mt2]=min(lc2); 

[nc2,*4mn2]  =  min(nc2); 

NB1=0; 

NB2=0, 

NB3=0; 

N  W 1  =0; 

NW2=0, 

NW3*0; 

NR1=0; 

NR2=0; 

NR3=0; 

NRymlel=0; 

NRymle2=0; 

NRymle3=0; 

NRymoml»0; 

NRymom2«0; 

NRymom3=0; 
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NRllssO; 

NR12»0; 

NRl3=sO, 

NLlsO; 

N  L2~0; 

NL3=0; 

N  N  1=0; 

NN2  =  0. 

NN3  =  0; 

»Ouml  s=4ori(#0um); 

for  i  =  1  :mtg 

if  jOum l(i )4  =  *ambl/200 
NB 1  =  NB 1  +  1 ; 
end 

if  «0um l(i)i»*mb! /200&«0uml(i)i**amb2/200 
NB2=  NB2+ 1 ; 
end 

if  *0uml(i)^*amb2/200 
NB3=NB3+1; 
end 

if  »0um)(i)i  =  5Amwl/200 
NW1  =  NW1  +  1; 
end 

if  «0uml(i);famwl/200&«0uml(i)i=»amw2/200 
NW2=NW2+1, 
end 

if  *0um)(t);»amw2/200 
NW3  =  NW3+1; 
end 

if  *0uinl(i);s=**mrl/200 
NRl  —  NRl  +  1; 
end 

if  *0uml(i)i*amr1  /200<£#0uml(i)i**an»r2/200 
NR2=NR2+1; 
end 

if  *0uml(i )i«amr2/200 
NR3=NR3+1; 
end 

if  •Ouml(i)is«*imyml«  1/200 
NRymlel  xNRymlel  +  1 ; 
end 

if  fOuml(i)t,»*mrymle1/200&«Oaml(i)is*ARirymle2/200 
NRy2xNRy2+l; 
end 

if  »0uml(i)^»amrymle2/200 
NRymle3sjNRymle3+  1 ; 
end 

if  i0«ml(i)is»mr)rmoml/300 
NRymomlsNRymoml  +  l; 
end 

if  •0aml(i)4«*fnryfnom1/200<e*0uml(i)ise*inrymoin2/200 
N  Ry  mom]E  NRy  mom2+ 1 ; 
end 

if  •0aml(i)i»»mryfnom2/200 
NRymom3*NRymom3  +  l, 
end 

if  *Oaml(i)i««Binrl  1/200 

NR1I-NA11  +  1; 


end 

if  *Ouml(i)ti*mrl l/2Q0&*Guml{i)i=i*mrl2/200 
NR12  =  NR12+1; 
end 

if  *Ouml(i)i«*mrJ 2/200 
NR13=NR13+1; 
end 

if  •0uml(i)i=*amll/200 
NL1  =  NL1  +  1; 
end 

if  *0uml(i )i,*anr>Il/200Ari0uml(i)i=**ml2/200 
NL2=NL2+1; 
end 

if  »0uml(i)i*aml2/200 
NL3=NL3+1; 
end 

if  aOuml(i)i  =  «*mn  1/200 
NNlsNNl-f  1; 
end 

if  *0uml(i)i«*mnl/200&*0uml(i );=*amn2/200 
NN2=NN2+1, 
end 

if  «0uml(i)i«amn2/200 
NN3=NN3+1; 
end 
end 

TSBl=((NBl-(l/3)“mag)'2)/((  1/3)“  mng); 
TSB2=((NB2-(l/3)-m*g)*2)/((l/3)*mag); 
TSB3=((NB3-(l/3)*roag)'2)/((  1/3)*  mag); 

TSB=TSB1 +TSB2  -f  TSB3 

TSWl  x((N  Wl-(1  /3)*mag)'2)/((  1  /3)*mag); 

TSW2x((NW2-(l/3)*n»»g)‘2)/((  1/3)"  mag); 

TSW3s((N  W3-(l/3)*mag)‘  2)/((  1  /3)*mag); 

TSW=TSWl-fTSW2+TSW3 

T$Rlx((NRl-(  l/3)“m*g)'2)/((  l/3)*mag); 

TSR2*((NR2-(l/3)-mag)*2)/((1/3)-mag); 

TSR3=((NR3-(  l/3)"mag)‘2)/((  l/3)“mag); 

TSR*TSRl+TSR2+TSR3 

TSRymUl=((NRym\el-(l/3)*mag)'2)/((l/3)-mag); 
TSRymle2«((NRymle2-(  1  /3)*mag)'2)/((  l/3)“mag); 
TSRyrole3=((NRymle3-(  1  /3)“mag)*2)/((l  /3)*mag); 
TSRyrale=TSRymlel  +TSRymle2-fTSRymU3 
TSRymoml=((NRymoml-(l/3)mmag)‘2)/((l/3)*inag); 
TSRjrmom2=((NRymom2-(  1  /3)“mag)‘2)/((l  /3)“mag); 
TSRymom3»({  NRymom3-(  1  /3)*mag)‘  2)J  ((  1  /3)“mag); 
TSRymom=TSRymoml  -f  TSRymom2+TSRymom3 
TSRlla((NRH-(l/3)-mag)  2)/((l/3)-mag); 
TSRl2=((NRi2-(  1  /3)“mag)*2)/((  1  /3)*mag); 
TSRl3x((NRl3-(  I  /3)*mag)‘  2)/((  1/3)* mag); 

TSRl  »TSRl  1  -f  TSR12+TSR13 
TSLl*((  NLl-(  l/3)*mag)'  2)/((  1  /3)*mag); 
TSL2c((NL2-(  J/3)“mag)‘  2)/((  1  /3)“mag); 
TSL3*((NL3-(  l/3)“mag)‘2)/((  1/3)“ mag); 

TSL*TSLl +TSL2+TSL3 
TSNl=s((NMl-(l/3)*m*g)‘2)/((  1/3)*mag); 
TSN2at((NN2*(l/3)“mng)*  2)/((  l/3)*mng); 
TSN3*((NN3-(l/3)“wag)*  2)/((  1/3)"  mag); 

TSNxTSNl  +  TSN2-f  TSN3 
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end 

if  20,  =  m*g4fn.*g;2S 
(or  is  1.200 

bcl(i)=*o«(BPl(i)~0.25), 
bc2(i  )  =  *b»(  BP  l(i  )-0.5); 
bc?(i)=»b*(BPl(i)-0  75); 
wcl(i)=*b*(  WPl(i)-0  25); 
wc2(i)=*b*(  WPl(i)-0.5); 
wc3(i)**b«(  WPl(i)-0  75); 
rcl(i)=*b*(  RPl(i)-0  25); 
rc2(i)s=*b*(RP  l(i)>0.5); 
rc3(i)s=*b#(RPl(i)-0.75); 
rymlec  1(  i  )  =  *b«(  RymleP  1  (»)*0. 26); 
rymlec2(i)=*b<(RymleP  1  (i)-0  5); 
rymlec3(i)=ab*(  RymleP  l(»)-0. 75); 
rymomc  l(i  )  =  *bi(  RymomP|(i)-0  25); 
rymomt2(i)s*b*(  Ry  rnomP  l(i)-0.5), 
rymomc3(i)=»bi(  RymomPl(i)-0.75); 
rlcl(t)s*bft(fUPl(t)-0.25). 
rlc2(i)=*b*(RlPl(i).o  5); 
r!c3(i)=»b*(RlPl(i)-0  75); 
lcl(i)  =  4b*(LPl(i)-0  25); 
k2(i)=4b«(LPj(i)-0.5), 
lc3(i)=*b*(LPl(i)-0.75), 
nc  l(i)s*bs(NP  t(i)«0.25); 
nc2(i)sab«(NPl(i)*0  5); 
nc3(i)s=4b*(NPl(i)-0.75); 
end 

[bcl  ,j»mbl]=min(  be  1 ), 

[wcl  ,i4mwl)=min{  wc  1 ); 

[rcl  ,aamrl]=min(rcl ); 

(rymlec  1  .•amrym)e)]  =  min(rymlecl ); 
(rymomcl  ,umrymoml]  =  min(  rymomcl ); 
(r  1c  1  ,iamrl  1)  =  min(  r  1c  1  ), 

Del  ,#»mll)=min(lc  1 ), 

(ncl,umnl)  =  min(ncl), 

[bc2,i*mb2)  =  min(  be 2); 

[wc2,umw2)zmin(  wc2), 
|rc2,#*mr2)=min(rc2); 

(rymlec2,«amrymle2]smin(ryml«c2); 

(rymomc2,*»mrymom2]  =  min(  rymomc2); 

frlc2,*antrl2]acmin(rlc2); 

Ilc2,»*ml2)amin(lc2); 

[nc2,«*mn2)=min(  nc2), 
(bc3,**mb3]=min(bc3), 

[wc3,umw3]  =  miB(  wc3); 

[rc3,«*mr3]at  min(  rc3 ); 
(rymlec3,*amrymle3]smiTt(rytn)ec3); 
[ij,momc3,umrymom3|i=min(  ry  momc3 ), 

(r  lc3,aamr  13]»min(  rlc3), 
[Ie3,»4ml3)«min(lc3), 

[nc3,*4mn3]»min(  ne3), 

NB1-0; 

N02»O; 

NB3*0; 

HB4-0; 

NWUO; 

NW2«0; 


NW3=0; 

NW4=0; 

NRl=0; 

NR2=0; 

NR3=0; 

NR4=0; 

NRy  dale  1  =sO; 

NRymle2=sO; 

NRymle3=sO; 

NRymle4sO; 

NRy  mom  1=0; 

NRymom2=0; 

NRymom3=0, 

NRymom4=0; 

NRll=0; 

NR12=0; 

NR13=0; 

NR14=0; 

NLl=0; 

Nt2=0; 

NL3=0; 

NL4=0; 

NN  1=0; 

NN2=0; 

NN3=0; 

NN4  =  0; 

iOBml=iort(iOum); 
for  isl  m*| 

if  «0uml(i);s«4mb  1/200 
NBlsNBl-f  1; 
end 

if  *0uml(i)i.«4mbl/2004f»0uml(i)i=:»4mb2/200 
NB2=NB2+1; 
end 

if  *0um  l(i )i*4mb2/200&»0um](i ),=*amb3/200 
NB3*NB3-rl; 
end 

if  *0uml(i)itAmb3/2OO 
NB4  =  NB4  +  1, 
end 

if  iOuml(i)i=umw  1/200 
NW1  =  NW1  +  1; 
end 

if  «Owm l(i)^,«4mw  1/200 iri0um!fi);  —  «Amw2/200 
NW2*NW2+l; 
end 

if  «0u  m  l(i  )i«amw2/200 Ic  *0n  m  1(  i  );**amw3/?00 
NW3=NW3+], 
end 

if  *0nm](i)ic*mw3/200 
NW4.NW4+1; 
end 

if  «0uml(i);mft*mt!/200 
NRl*NRl  +  1; 
end 

if  *0«ml(i)^»mrl/2004te0ttm](i)i«*4inf2/200 
NR2-NR241; 
end 
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if  •0uml(i)i**mr2/200&i0uml(i)j=4amr3/200 
NR3  =  NR3+1; 
end 

if  s0uml(i)i*ainr3/200 
NR4  =  NR4+|; 
end 

if  s0uml(i)i=*amrymiel/200 
NRymlel  =  NRymlel  +  l; 
end 

if  s0uml(i)t*»mrymlel/200&a0uml(i)i=:a».mrymle2/200 
NRymle2=NRymle2+l; 
end 

if  «0uml(i)i#*mrymle2/200&80iiml(i)i=5*mrymle3/200 
NRymle3=NRymle3+l; 
end 

if  «0uml(i)is»mrymle3/200 
NRymle4  =  NRymle4-fl; 
end 

if  s0uml(i)i=»amrymoml/200 
NRy  moml  =  NRymom  1  +  1; 
end 

if  »Ouml (i)4»*mry moml /200^*0uml(i)is=«nmrymom2/ 200 
NRymom2=NRymom2+l; 
end 

if  «0uml(i)^8«mrymom2/200&;s0uml(i)i=s3mrymom3/200 
NRymom3=NRymom3+l; 
end 

if  ft0uml(i)i*amrymom3/200 
NRymom4=NRymom4  +  l; 
end 

if  «Ouml(i)i=«amrl  1/200 
NRll  =  NRil  +  l; 
end 

if  80uml(i)i*amrll/2004£«0uml(i)j=s*mrl2/200 
NRl2=NRl2+l; 
end 

if  s0uml(i)<.samrl2/200&s0u ml (i)j=samrl 3/200 
NR13=NR13+1; 
end 

if  s0oml(i)i*amrl3/200 
NR14=NR14  +  1; 
end 

if  «0aml(i)i=samll/200 
NL1=NL1  +  1; 
end 

if  sOuml(i )i*amll/200&s0iiml(i  ji=*aral2/200 
NL2=NL2+1; 
end 

if  «0uml(i)i*aml2/200&*0uml(i);=«anil3/200 
NL3=NL3+1; 
end 

if  &0uml(i)i&aml3/200 
NL4  =  NL4  +  1; 
end 

if  «Ouml(i)i=:«amnl /200 
NNl=NNl  +  l; 
end 

if  *0uml(i)i*amnl/200&*0aml(i)i=samii2/200 


NN2=NN2+1; 

end 

if  iOuml(i)i«4mB2/200&iOuml(i)is:tama3/200 
NN3=NN3+J; 
end 

if  s0uml(i)i*4mn3/200 
NN4=NN4+1; 
end 
end 

TSB1  =((NBl-0.25*mag)'2)/(0.25*ni4g); 

TSB2=((NB2-0.25*m4g)‘2)/(0.25*»ag); 

TSB3=((NB3-0.25*m4g)*2)/(0.25*m4g); 

TSB4=((NB4-0.25*mag)*2)/(0.25*mag); 

TSB=TSB1+TSB2+TSB3+TSB4 

TSWl=((NWl-0.2S*mag)*2)/(0.2S*ro4g); 

TSW2=((NW2-025*mag)*2)/(0.25*m4g); 

TS  W3=((N  W3-0.25*m*g)  *  2)/(0.25*mag); 

TSW4=((NW4-0.25*mag)'2)/(0.25*mag); 

TSW=TSW1+TSW2+TSW3+TSW4 

TSRl=((NRl-0.25*mag)‘2)/(0.25*mag); 

TSR2=((NR2-0.25*m*g)*2)/(0.25*m4g); 

TSR3=((NR3-0.25*m»g)‘2)/(0.25*mag); 

TSR4=((NR4-0.25"mag)‘2)/(0.25*anag); 

TSR=TSRl-f  TSR2+TSR3+TSR4 

TSRymlel=((NRymlel-0.25*mag)'2)/(0.25**nag); 

TSRymle2=((NRymle2-0.2S*mag)‘2)/(0.2&*mag); 

TSRymle3=((NRymle3-0.25*mag)'2)/(0.2«i*in4g); 

TSRymle4=((NRymle4-0.25*ra4g)*2)/(0.2i>*inag); 

TSRymle=TSRymlel+TSRymle2+TSRymle3+TSRymle4 

TSRymoml=((NRymoml-0.25*mag)*2)/(0.25*mag); 

TSRymom2s=((NRymora2-0  25*m4g)'2)/(0.25*mag); 

TSRymora3=((NRymom3-0.25*m4g)*2)/(0.25",m4g); 

TSRymom4=((NRymom4-0.25*m4g)%2)/(0.25*mag); 

TSRymom=TSRymoml+TSRym4m2+TSRymom3+TSRymom4 

TSRll=((NRll-0.25*mag)-2)/(0.25*m4g); 

TSR12=((NR12-0.25*mag)*2)/(0.25»mag); 

TSR13=((NRl3-0.25*mag)*2)/(0.25*mag); 

TSRU=((NR14-0.25*mag)*2)/(0.25*m4g); 

TSR1=TSR11+TSR12+TSR13+TSR14 

TSLl=((NLl-0.25*m4g)*2)/(0.25*mag); 

TSL2=((NL2-0.25*mag)*2)/(0.25*mag); 

TSL3=((NL3-0.25*mag)'2)/(0.25*inag); 

TSL4=((NL4-0.25*m4g)*2)/(0.25*m4g); 

TSL=TSL1+TSL2+TSL3+TSL4 

TSNI=((NNl-0.2&*m4g)'2)/(0.25*mag); 

TSN2=((NN2-0.25*mag)‘2)/(0.25*mag); 

TSN3=((NN3-0.25*mag)'2)/(0.25"mag); 

TSN4=((NN4-0.25*m4g)‘2)/(0.25"m4g); 

TSN=TSN1+TSN2+TSN3+TSH4 

end 

if  m*gt— 25 
for  i= 1:200 

bcl(i)=4bi(BPl(i)-0.2); 

bc2(i)«ab*(BPl(i)-0.4); 

bc3(i)=aba<BPl(»)-0.6); 

bc4(i)=abi(BPl(i)-0.8); 

wcl(i)=ab*(  WPl(i)-0.2); 

wc2(i)=ab*(WPl(i)-0.4); 


wc3(i)=ab*(WPl(i)-0.<l); 

wc4(i)=4b«(WPl(i)-0.*); 

rcl(i)=ab*(RPl(i)-0.2); 
rc2(i)=ab*(RPl(i)-0.4); 
rc3(i)=ab*(RPl(i)-0.8); 
rc4(i)=ab«(RPl(i)-0.*); 
ryralecl(i)=ab*(RymlePl(i)-0.2); 
rymlee2(»)=abt(RymlePl(i)-Q.4); 
rymlec3(i)=Aba(RymlePl(i)-O.0); 
rymlec4(i)=ab«(RymlePl(i)-0.d); 
rymomcl(i)=abi(  RymomPl(i)-0.2) 
rymomc2(i)=*b»(RyznomPl(i)-0.4) 
rymomc3(i)=ab»(RyraomPl(i)-0.6) 
rymomc4(i)=ab*(RymoinPl(i)-0.8) 
rlcl(i)=ab#(RlPl(i)-0.2); 
rlc2(i)=ab»(RlPl(i)-0.4); 
rlc3(i)=ab»(RlPl(i)-0.6); 
rlc4(i)=»b$(RlPl(i)-0.8); 
lcl(i)  =  abs(LPl(i)-0.2); 
lc2(i)a=ab#(LP  l(i)-0.4); 
lc3(i)=*b»(LPl(i)-0.«); 
lc4(i)  =  abs(LPl(i)-0.8); 
ncl(i)=ab*(NPl(i)-0.2); 
nc2(i)=ab«(NPl(i)-0.4); 
nc3(i)=ab*(NPl(j)-0.6); 
nc4(i)=ab»(NPl(i)-0.8); 
end 


(bcl,*ambl}=:min(bcl); 
fwcl  ,samw  l]=smin(  wcl  ); 

(rcl,samrl]  =  min(rcl); 

[rymlec  1  ,*»mryml«l]smin(ry  ralecl ); 
[rymomcl  ,»*mrymoml)  =  min(rymomcl ); 
[rlcl.samrl  l]  =  min(rlcl); 
(lcl,*amll]=min(lcl); 

(ncl(samnl]  =  min(ncl); 

[bc2,«amb2]  =  min(bc2); 
(wc2rsamw2]=:min(wc2); 

[rc2,#amr2]  =  min(rc2); 
[rymlec2,a*mrymle2]=min(ryin!ec2); 
[rymomc2(*amrymom2]  =  min(rymomc2); 
(rlc2,*amrl2]=min(rlc2); 
(Ic2,saml2]=min(lc2); 

(nc2,*amn2]  =  min(nc2); 
(bc3,«amb3]=min(bc3); 

(wc3,samw3]  =  min(wc3); 
(rc3>Bamr3]=inin(rc3); 
{rymlec3,«amrymle3]  =  min(rymlec3); 
(rymomc3,»amrymom3]=min{rymomc3); 
[r  Jc3,#*mrl3]  =  min(rlc3); 
[Ic3,«aml3]=min(lc3); 

(nc3,aamn3j  =  min(nc3); 

(bc4,*amb4]smtn(bc4); 

(wc4,*amw4)=min(wc4); 

[rc4,«amr4]=min(rc4); 

[rymlec4,aamrymle4]  =  min(rynilcc4); 

(rymomc41i»mrymom4]cmin(rymomc4); 

(rlc4,#amrl4)  =  mtn(rlc4); 

(Ic4,jaml4]  =  min(lc4); 


(nc4,a*mn4]  =  min{  nc4); 

NBl  =0; 

NB2=0; 

NB3=0; 

NB4=0; 

KB5=0; 

NWl=0; 

NW2=0; 

NW3=0; 

NW<=0; 

NW5=0; 

NR1=0; 

NR2=0; 

NR3  =  0; 

KR4=0, 

NR6=0; 

NRymIel=0; 

NRymle2=0; 

NRymle3=0; 

NRymle4=0; 

NRymle5=0; 

NRymoml=0; 

NRymom2=0; 

KRymom3=0; 

NRymom4=0; 

NRymom5=0; 

NRll  =  0; 

NR12=0; 

NR13=0; 

NR14=0; 

NR15=0; 

NL1=0; 

NL2=0, 

NL3=0; 

NL4=0; 

NL3=0; 

NN1=0; 

NN2=0; 

NN3=0; 

NN4=0; 

NN5=0i 

*0uml  =  *orl(»0um); 
for  i  =  l:mag 

if  «0aml(i)i=sainb]/200 
NBl=NBl  +  l; 
end 

if  s0uml(i)i*ambl/200&s0uml(i);:=:8amb2/200 
NB2=NB2+1; 
end 

if  s0ural(i);«atnb2/200&*0uml(i)j=«amb3/200 
NB3=NB3+1; 
end 

if  «0uml(i)t*amb3/200&*0uml(i);s«amb4/200 
NB4  =  NB4+J; 
end 

if  s0uml(i)^.«amb2/200 
NB5*NB5+1; 
end 


B-ll 


if  s0uml(i)i=*amwl/200 

NWlsNWl+1; 

end 

if  s0\iral(i)i,*amwl/200&*0uml(i)i=*amw2/200 
N  W2=sN  W2+ 1; 
end 

if  «0umj(i)i,samw2/200i£«0uml(i)j=»amw3/200 
NW3=NW3+1; 
end 

if  30uml(i)i*amw3/200&*0uml(i)i5=*amw4/200 
NW4  =  NW4  +  1; 
end 

if  s0uml(i)ii*mw4/200 
NW5sNW5  +  1; 
end 

if  s0uml(i)j=#amrl/200 
NRlsNRl  +  3; 
end 

if  sOuml(i)isamrl/2OO&0Ouml(i);=#ainr2/2OO 
NR2sNR2+l; 
end 

if  s0uml(i)i»amr2/200&*0uml(i)i=*amr3/200 
NR3=NR3+1; 
end 

if  30uml(i)is*mr3/200&*0uml(i)i=:**mr4/200 
NR4=NR4  +  1; 
end 

if  s0uml(i)i*amr4/200 
NR5=NR5-fl; 
end 

if  iOuml(i)»=«amrymle  1/200 
NRyml«l=NRymlel  +  l; 
end 

if  50uml(i)t*amryralel/200&fl0uml(i)i=samrymle2/200 
NRymle2=NRymie2+l; 
end 

if  «0uml(i)i«&mrymle2/200<cft0uml(i)is*»mrymle3/200 
NRymle3=NRymle3+l; 
end 

if  sOuml(i)i*amrymle3/2OO&*Ouml(i)i=0amrymle4/2OO 
NRymIe4  =  NRymle4+l; 
end 

if  s0uml(i)i,«*mrymle4/200 
NRymle5=NRymle5+l; 
end 

if  «0uml(i)i=«amrymoml/200 
NRymoml  =  NRymomi+li 
end 

if  *0uml(i)i**mrymoml/200&«0uml(i);s»annrymom2/200 
NRymom2=NRymom2+l; 
end 

if  0Ouml(i)isamrymom2/2OO&«Ouml(i)i=«a(nrymom3/2OO 
NRymom3=NRymom3+l ; 
end 

if  0Ouml(i)i.0amrymom3/2OO&0Ouml(i)is:0amrymom4/2OO 
NRymom4=NRymom4  +  l; 
end 

if  «0uml(i)i«amrymom4/200 


N  Rymom5=NRy  mom5+ 1 ; 
end 

if  sOuml(i)i=t&mrll /200 
NRU  =  NR11  +  1, 
end 

if  sOuml(i)j,0*mrll/2OO&sOuinl(i)j:=«*mrl2/2OO 
NR12=NR12+1; 
end 

if  s0uml(i);«4mrl2/20O&»0uml(i)i=*aTnrl3/2O0 
NR13=NR13+1; 
end 

if  «Ounnl(i)^04mrl3/2OO<£sOuml(i)i=0amrl4 / 200 
NR14=NR14+1; 
end 

if  0Ouml(i)isaiml4/2OO 
NR15=NR1S+1; 
end 

if  s0uml(i)i=«amll/200 
NLls=NLl  +  l; 
end 

if  *Ouml(i);.«amll/200&5Guml(i)i=saml2/200 
NL2=NL2+1; 
end 

if  «0uml(i);sam)2/200&$0uml(i);=:s&inl3/200 
NL=NL3+1; 
end 

if  sOuml(i)i,saml3/2OO&0Ouml(i)i=*aml4/2OO 
NL4=NL4-f  1; 
end 

if  flOuml(i)^0annl4/2OO 
NL5=NLS+1; 
end 

if  s0timl(i)i=s»mn]/200 
NNl=NNl  +  I; 
end 

if  *0uml(i);.samnl/200&s0uml(i)i=*amn2/200 
NN2=NN2  +  1; 
end 

if  sOurol(i)£.0amn2/2OO&sOuml(i)j=«amn3/2OO 
NN3=NN3+1; 
end 

if  s0uml(i);.samn3/2G0££«Guml(i)i=ftafnn4/2OO 
NN4=NN4  +  1; 
end 

if  <0um!(i)^,samn4/200 
NN5=NN5  +  1; 
end 

end 

TSBl=((NBl-0.2*mag)‘2)/(0.2*mag); 

TSB2=((NB2-0.2*mag)‘  2)/(0.2*mag); 

TSB3=((NB3-0.2*mag)2)/(0.2*mag); 

TSB4=((NB4-0  2*mag)*2)/(0.2*mag); 

TSB5=((NB3-0.2»mag)*2)/(0.2*mag); 

TSB=TSBl+TSB2+TSB3+TSB4+TSB5 

TSWl=((NWl-0.2*mag)*2)/(0.2*mag); 

TS  W2=s((NW2-0.2*mag)‘2)/(0.2*mag); 

TSW3=((NW3-0.2*mag)‘2)/(0.2*mag); 

TSW4=((NW4-0.2*mag)-2)/(0.2*mag); 
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TSWJ=((NW5-0.2*m»|)'5)/(0.2”iM|); 

TSW=TSW1+TSW2+TSW3+TSW4+TSWS 

TSRl=((NRl-0.2*mag)‘2)/(0.2-mag); 

TSR2=((NR2-0.2-mag)*2)/(0.2-inag); 

TSR3=((NR3-0.2’raig)‘2)/(0.2*in«.g); 

TSR4=((NR4-0.2“m»g)-2)/(0.2-m»g); 

TSR5=((NR5-0  2-mag)-2)/(0.2“mag); 

TSR=TSRl  +  TSR2+TSR3+TSR4+TSR5 
TSRymlel=((NRymlel-0.2-mag)-2)/(0.2*mag); 
TSRymle2=((N  Ry  mle2-0.2*mag) '  2)/(0.2*mag); 
TSRymle3=((NRymle3-0.2*mag)*2)/(0.2-mag); 
TSRymle4=((NRymle4-0.2*mag)‘2)/(0.2*mag)i 
TSRymle5=((NRymle5-0.2-mag)-2)/(0.2-mag); 
TSRymIe=TSRymlel+TSRymle2+TSRymle3 
+TSRymle44  T3Rymle5 

TSRymoml=((NRy  mom  1-0.2-mag)  *2)/(  0.2-mag); 
TSRymom2:=((NRymom2-0.2*mag)  '2)/(0.2*mag); 
TSRymom3=((NRymom3-0.2-mag)'2)/(0. 2-mag); 
TSRymom4=((NRymom4-0.2-mag)'2)/(0.2-mag); 
T3Rymom5=((  NRymom5-0.2-mag)'2)/(0.2-mag); 
TSRymom=TSRymoml+TSRymom2+TSRymom3 


+  TSRymom4  +  TSRymom5 

TSRIl=((NRn-0.2*m»g)-2)/(0.2-m»g); 

TSR12=((NR12-0.2*m»g)-2)/(0.2'mag); 

TSR13=((NR13-0.2*m»g)-2)/(0.2,m»g); 

TSR14=((NR14-0.2’m»g)-2)/(0.2“rn»g); 

TSR15=((NR15-0.2-mag)'2)/(0.2"mag); 

TSR1=TSR11+TSR12+TSR13+TSR14  +  TSR15 

TSLl=((NLl-0.2*raag)-2)/(0.2'mag); 

TSL2=((NL2-0.2*mag)-2)/(0. 2-mag); 

TSL3=((NL3-0.2*mag)-2)/(0.2-mag); 

TSL4=((KL4-0.2,m»g)-2)/(0.2*m»g); 

TSLS=((NL5-0.2*m»g)-2)/(0.2*m»g); 

TSL=TSL1  +  TSL2+TSL3+TSL4+TSL5 

TSNl=((NNl-0.2'mag)-2)/(0.2*mag); 

TSN2=((NN2-0.2*mag)‘2)/(0  2' mag); 

TSN3=((NN3-0.2*mag)-2)/(0. 2-mag); 

TSN4=((NN4-0.2*mag)'2)/(0.2'mag); 

TSN5=((NNS-0.2*mag)'2)/(0.2'mag); 

TSN=TSN1+TSN2+TSN34TSN4+TSN5 

end 
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